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c)  Case-9,  A2  =  0.05  and  As  =  0.001  216 


6.45  Contours  of  spanwise  vorticity  a;,  in  the  xy-plane  for  t  =  4096At. 
z  =  0.  Twenty  contour  intervals  between  the  minimum  and 
maximum  values.  Solid  lines  denote  positive  vorticity  and  dashed 


lines  denote  negative  vorticHy. 

a)  Case-7,  A2  =  0.00  and  A3  =  0.001; . 217 

b)  Case-8,  A2  =  0.01  and  A3  =  0.001; . 218 

c)  Case-9,  A2  =  0.05  and  As  =  0.001  219 
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6.46  Contours  of  spanwise  vorticity  u;,  in  the  xy-plane  for  t  =  4096At 

and  z  =  0.  Case-8:  A2  =  .01,  Aj  =  .001,  0  =  .51,  and  7  =  .5. 

Thirty  contour  intervals  between  the  minimum  and  maximum 
values.  Solid  lines  denote  positive  vorticity  and 

dashed  lines  denote  negative  vorticity . 220 

6.47  Contours  of  spanwise  vorticity  u>t  in  the  zz-plane  for  t  =  4096At. 

Case-8:  =  .01,  A3  =  .001,  0  =  .51,  and  7  =  .5. 

Thirty  contour  intervals  between  the  minimum  and  maximum 
values.  Solid  lines  denote  positive  vorticity  and  dashed  lines 
denote  negative  vorticity. 


a)  y  =  2  (y/Ay  =  72), .  221 

b)  2/  =  -2  (y/Ay  =  56)  .  222 


6.48  Contours  of  streamwise  vorticity  in  the  zz-plane  for  i  =  4096A<. 
Case-8:  A2  =  .01,  A3  =  .001,  0  =  .51,  and  7  =  .5. 

Thirty  contour  intervals  between  the  minimum  and  maximum 
values.  Solid  lines  denote  positive  vorticity  zmd  dashed  lines 
denote  negative  vorticity. 


a)  y  =  2  (y/Ay  =  72), .  223 

b)  y  =  -2  (y/Ay  =  56)  .  224 


6.49  Contours  of  spanwise  vorticity  u;,  in  the  yz-plane  for  i  =  4096At. 
Ceise-8:  A^  =  .01,  .43  =  .001,  0  =  .51,  and  7  =  .5. 

Thirty  contour  intervals  between  the  minimum  and  maximum 
values.  Solid  lines  denote  positive  vorticity  and  dashed  lines 


denote  negative  vorticity. 

a)  z  =  .172  (z/Az  =  110), . 225 

b)  I  =  .160  (z/Az  =  101) . 226 


6.50  Contours  of  streamwise  vorticity  u;,  in  the  yz-plane  for  t  =  4096 At. 
Czise-8:  A2  =  .01,  A3  =  .001,  0  =  .51,  and  7  =  .5. 

Thirty  contour  intervals  between  the  minimum  and  maximum 
values.  Solid  lines  denote  positive  vorticity  and  dashed  lines 


denote  negative  vorticity. 

a)  z  =  .172  (z/Az  =  110), . 227 

b)  z  =  .160  (z/Az  =  101) . 228 

6.51  Vorticity  amplitude  distributions  for  F  =  1.  Case-7:  A2  =  0, 

A3  =  .001,  0  =  .51,  and  7  =  .5. 

a)  spanwise  vorticity,  fl®; . 229 

b)  spanwise  vorticity,  flj; . 230 

c)  streamwise  vorticity,  . 231 
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ABSTRACT 

Laminar-turbulent  transition  of  high-defidt  fiat  plate  weikes  is  investigated 
by  direct  numerical  simulations  using  the  complete  Naver-Stokes  equations.  The 
simiilations  are  based  on  a  spatial  model  so  that  both  the  base  fiow  and  the 
disturbance  flow  can  develop  in  the  downstream  direction.  The  Navier-Stokes 
equations  are  used  in  a  vorticity-velodty  form  and  are  solved  using  a  combination 
of  finite-difference  and  spectral  approximations.  Fourier  series  are  used  in  the 
spanwise  direction.  Second-order  finite-differences  are  used  to  approximate  the 
spatial  derivatives  in  the  streamwise  and  transverse  directions.  For  the  tempo¬ 
ral  discretization,  a  combination  of  ADI,  Crank-Nicolson,  and  Adams-Bashforth 
methods  is  employed.  The  discretized  velocity  equations  are  solved  using  fast 
Helmholtz  solvers.  Code  validation  is  accomplished  by  comparison  of  the  numeri¬ 
cal  results  to  both  linear  stability  theory  and  to  experiments. 

Calculations  of  two-  and/or  three-dimensional  sinuous  mode  disturbances  in 
the  wake  of  a  fiat  plate  are  undertaken.  For  calculations  of  two-dimensional  dis¬ 
turbances,  the  wake  is  forced  at  an  amplitude  level  so  that  nonlinear  disturbance 
development  may  be  observed.  In  addition,  the  forcing  amplitude  is  varied  in 
order  to  determine  its  effect  on  the  disturbance  behavior.  To  investigate  the 
onset  of  three-dimensionality,  the  wake  is  forced  with  a  small- amplitude  three- 
dimensional  disturbance  and  a  larger  amplitude  two-dimensional  disturbance.  The 
two-dimensional  forcing  amplitude  is  varied  in  order  to  determine  its  influence  on 
the  three-dimensional  fiow  field. 

Two-dimensional  disturbances  are  observed  to  grow  exponentially  at  small 
amplitude  levels.  At  higher  amplitude  levels,  nonlinear  effects  become  important 
and  the  disturbances  saturate.  The  saturation  of  the  fundamental  disturbance 
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appears  to  be  related  to  the  stability  characteristics  of  the  mean  flow.  Larger 
forcing  amplitudes  restilt  in  the  earlier  onset  of  nonlinear  effects  amd  saturation. 
At  large  amplitude  levels,  a  Karman  vortex  street  pattern  develops. 

When  the  wake  is  forced  with  both  two-  and  three-dimensional  disturbances, 
strong  interactions  between  these  disturbances  is  observed.  The  satxiration  of  the 
two-dimensional  disturbance  causes  the  three-dimensional  disturbance  to  saturate. 
However,  this  is  followed  by  a  resumption  of  strong  three-dimensional  growth 
that  may  be  due  to  a  secondary  instability  mechanism.  Larger  two-dimensional 
forcing  amplitudes  accelerate  the  saturation  of  the  two-dimensional  eind  three- 
dimensional  disturbances  as  well  as  accelerate  the  resumption  of  strong  three- 
dimensional  growth.  These  interactions  also  result  in  complicated  distributions  of 
vorticity  and  in  a  significant  increase  in  the  wake  width. 
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CHAPTER  1 

INTRODUCTION 

According  to  Sato  and  Kuiiki  (1961),  in  the  laminar-tiirbulent  transition  of  a 
flat  plate  wake  three  distinct  regions  can  be  identified:  a  linear  region,  a  nonlinear 
region,  and  a  three-dimensional  region.  In  the  linear  region,  small  disturbances 
in  the  background  flow  trigger  two-dimensional  wave-like  fluctuations.  Sato  and 
Kuriki’s  measurements  indicate  that  the  frequency  of  the  dominant  fluctuation  in 
the  wake  corresponds  to  the  frequency  of  the  most  amplified  disturbance  predicted 
by  linear  stability  theory.  Once  present  in  the  wake,  these  fluctuations  grow  ex¬ 
ponentially  in  the  downstream  direction  and  are  well  described  by  linear  stability 
theory. 

Due  to  the  large  ampliflcation  rates  associated  with  wakes,  the  amplitudes  of 
the  disturbances  rapidly  become  quite  large  and  nonlinear  interactions  between  the 
disturbances  become  important.  The  behavior  of  the  disturbances  deviate  from 
linear  stability  theory  predictions  but  the  wake  remains  two-dimensional  (Sato 
and  Kuriki,  1961).  In  this  region,  the  distiirbances  grow  at  rates  significantly  less 
than  those  predicted  by  linear  stability  theory  and  harmonics  of  the  fundament  d 
disturbances  are  generated.  Also  due  to  nonlinear  effects,  the  mean  flow  is  altered. 
As  verified  by  several  researchers  (Zabusky  and  Deem,  1971;  Aref  and  Siggia,  1981; 
and  Meiburg  and  Lasheras,  1988),  this  sequence  of  events  leads  to  the  well  known 
Kirm^  vortex  streets  that  are  observed  in  flat  plate  wakes. 

Finally,  three  dimensionality  becomes  important.  The  work  of  Robinson  and 
Saflman  (1982)  indicates  that  the  two-dimensional  Kirman  vortex  street  is  un¬ 
stable  with  respect  to  three-dimensional  disturbances.  As  a  consequence  of  this 
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instability,  tbree*dimensional  lambda  vortex  patterns  are  formed  in  the  region  be¬ 
tween  consecutive  Karman  vortices  of  opposite  sign  (Meiburg  and  Lasheras,  1988). 
1.1  The  Linear  Region  of  Transition 

Linear  stability  theory  describes  the  behavior  of  small  amplitude  disturbances 
in  parallel  shear  flows.  In  nonparallel  shear  flows,  linear  stability  theory  is  not 
strictly  applicable  but  can  still  be  a  useful  model  of  small  amplitude  disturbance 
development.  The  streamfunction  ^  of  a  disturbance  described  by  linear  stability 
theory  has  the  form 


(1,1) 

where  the  variables  x  =  {x,y,z)  and  t  are  the  spatial  position  vector  and  time 
respectively.  The  amplitude  <f>{y)  is  an  eigenfunction  of  the  Orr-Sommerfeld  equa¬ 
tion 

(l7-c)(#"-(a^+7’)^)-£l'V+ j^(r'-2(a'+7')*"+(o*+7“)V)  =  0  (1.2) 

where  primes  denotes  differentiation  with  respect  to  y.  In  equation  (1.2),  U  = 
U{y)  is  the  streamwise  velocity  distribution  of  the  undisturbed  flow,  a  is  the 
streamwise  wave  nximber  of  the  disturbance,  7  is  the  spanwise  wave  number  of 
the  disturbance,  and  is  the  temporal  frequency  of  the  disturbance.  The  phase 
velocity  of  the  disturbance  is  c  =  ^/a  and  Re  is  the  Reynolds  number  of  the 
undisturbed  flow.  In  general,  a  and  are  complex  and  7  is  real.  However, 
interest  is  usually  confined  to  two  cases,  namely,  to  the  case  of  spatially  growing 
disturbances  for  which  a  =  (ar,a{)  is  complex  and  P  is  real,  and  to  the  case  of 
temporally  growing  disturbances  for  which  /?  =  (/?r«/?t)  is  complex  and  a  is  real. 
A  complete  derivation  of  equation  (1.2)  is  given  by  White  (1974). 
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Wakes  are  inviscidly  unstable  due  to  the  inflection  points  in  their  velocity 
profile.  As  is  characteristic  of  flows  with  inflectional  profiles,  the  disturbances  ob¬ 
served  in  wakes  have  high  aimplification  rates.  Also,  because  the  streamwise  veloc¬ 
ity  distribution  of  wakes  has  two  inflection  points,  two  different  instability  modes 
exist.  The  sinuous  mode  (mode-1)  is  highly  amplified  with  an  anti-symmetric 
streamwise  velocity  distribution.  The  varicose  mode  (mode-2)  is  less  amplified  (at 
most  frequencies)  than  the  sinuous  mode  and  has  a  symmetric  streamwise  velocity 
distribution.  Because  of  the  higher  amplification  rates  associated  with  the  sinuous 
mode,  it  is  this  mode  that  is  usually  observed  in  experiments  (Sato  and  Kuriki, 
1961;  Sato,  1970;  Miksad  et  al.,  1982). 

The  validity  of  linear  stability  theory  in  describing  the  initial  stages  of  transi¬ 
tion  in  near  wakes  has  been  verified  experimentally.  Sato  and  Kuriki  (1961)  studied 
the  transition  of  a  flat  plate  wake  subject  to  acoustical  excitation.  The  results  of 
these  experiments  were  compared  to  linear  stability  calculations  of  temporally- 
amplifying,  sinuous  mode  disturbances.  The  temporal  amplification  rates  /?{,  ob¬ 
tained  from  linear  stability  theory,  were  converted  into  spatial  amplification  rates 
^  using  a  phase  velocity  transformation.  The  theoretical  eigenfunctions  com¬ 
pared  very  well  with  those  obtained  experimentally.  However,  the  experimental 
amplification  rates  did  not  compare  as  well  with  those  predicted  by  linear  stabil¬ 
ity  theory  because  the  phase  velocity  transformation  is  not  accurate  enough  for 
dispersive  waves. 

Mattingly  and  Criminale  (1972)  experimentally  investigated  the  behavior  of 
small  amplitude  disturbances  in  the  near  wake  of  a  thin  airfoil  (NACA  0003)  and 
compaired  their  resiilts  to  both  spatial  and  temporal  linear  stability  theory.  For 
comparison  of  temporal  amplification  rates  to  experimentally  obtained  amplifi¬ 
cation  rates,  the  Gaster  transformation  (Gaster,  1962)  was  employed.  The  ex- 
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perimentally  obtained  amplitude  distributions  and  amplification  rates  compared 
well  to  spatial  linear  stability  theory  predictions.  The  assumption  of  temporally- 
amplifying  disturbances,  and  using  the  Gaster  transformation,  did  not  produce  as 
good  an  agreement  with  the  experimental  resxilts.  The  lack  of  agreement  can  be 
expledned  by  the  fact  that  the  Gaster  transformation  is  based  on  the  assiimption 
of  small  zunplification  rates  |aj  <<  1,  which  is  not  satisfied  for  high-deficit  wakes. 
From  their  results,  Mattingly  and  Ctiminale  concluded  that  spatial  stability  the¬ 
ory  is  more  accurate  than  temporal  stability  theory  in  predicting  the  behavior  of 
small  amplitude  disturbances  in  the  near  wake  of  a  thin  airfoil. 

Mattingly  and  Criminale  also  found  that  the  stability  characteristics  of  a 
spatially  varying  wake  could  be  computed  by  the  use  of  a  quasi-uniform  (quasi¬ 
parallel)  assumption,  i.e.  the  wake  is  eissumed  to  be  locally  paurallcl.  The  stability 
characteristics  of  the  local  mean  velocity  profile  are  then  computed  using  the  Orr- 
Sommerfeld  equation.  Then,  the  streamwise  wave  number  a  is  obtained  as  a 
function  of  streamwise  distance. 

Miksad  et  al.  (1982)  have  shown  that  linear  stability  theory  is  also  valid  when 
there  is  more  than  one  dominant  disturbance  in  the  wake.  They  experimentally 
studied  the  early  development  of  two  small  amplitude  sinuous  mode  disturbances, 
of  different  frequencies,  in  the  wake  of  a  thin  airfoil.  Their  results  indicated  that, 
at  small  amplitudes,  both  disturbances  grew  exponentially  at  rates  predicted  by 
linear  stability  theory. 

1.2  The  Nonlinear  Region  of  Transition 

Due  to  the  large  amplification  rates  associated  with  high-deficit  wakes,  the 
linear  re^on  of  transition  in  these  flows  is  usually  quite  small.  With  the  appearance 
of  large  amplitude  disturbances,  nonlinear  effects  become  important  and  linear 
stability  theory  is  no  longer  applicable.  In  contrast  to  the  linear  region  where 


22 


the  disttirbaaces  grew  qiiite  rapidly,  the  disturbances  evolve  more  gradually  in 
the  nonlinear  region.  This  is  also  in  contrast  with  the  development  in  boundary 
layer  flows  where  rapid  transition  to  turbulence  occurs  once  nonlinear  effects  and 
three-dimensionality  become  important. 

Sato  and  Kuriki  (1961)  found  that  the  early  (upstre2un)  part  of  the  nonlinear 
region  is  similar  in  some  respects  to  the  linear  region  of  transition.  Although  the 
amplification  rates  of  the  disturbances  are  significantly  lower  than  those  observed 
in  the  linear  region,  their  experimental  results  indicate  that  the  amplitude  dis¬ 
tribution  of  the  fundamental  disturbance  does  not  change  significantly  until  the 
nonlinearity  has  persisted  for  some  distance  downstream.  However,  harmonics 
of  the  fundamental  disturbance  are  generated  and  the  second  harmonic  becomes 
the  dominant  disturbance  component  near  the  wake  centerline  (Sato  and  Kuriki, 
1961).  Nonlinear  effects  also  change  the  mean  flow,  leading  to  rapid  increases  in 
the  mean  centerline  velocity  and  the  mean  wake  half- width  (Sato,  1970).  The 
resulting  mean  streamwise  velocity  distribution  is  much  fuller  than  that  of  the 
undisturbed  flow,  indicating  a  greater  degree  of  fluid  mixing. 

As  the  disturbances  travel  farther  downstream  in  the  nonlinear  region,  the 
amplitude  distribution  of  the  fundamental  disturbance  deviates  from  the  shape 
predicted  by  linear  theory.  The  measurements  of  Sato  (1970)  indicate  that  the 
peaks  in  the  streamwise  velocity  amplitude  distribution  of  the  fundamental  dis- 
turbemce  shift  towards  the  outer  edge  of  the  wake  and  also  decrease  in  magnitude 
with  increasing  downstream  distance.  The  peak  value  in  the  streamwise  velocity 
amplitude  distribution  of  the  second  harmonic  component,  which  initially  grew 
rapidly,  also  begins  to  decrease.  Sato  speculates  that  the  energy  lost  by  the  fluctu¬ 
ations  might  be  transferred  to  the  mean  flow.  Sato’s  measurements  also  show  that 
the  mean  half- width  and  the  mean  centerline  Velocity,  which  initially  increased  in 
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the  nonlinear  region,  begin  to  decrease  when  the  fundamental  streamwise  velocity 
fluctuation  begins  to  decrease.  Sato  showed  that  the  production  of  fluctuation 
energy 


—U'V 


dy 


(1.3) 


where  u'  and  v'  are  the  time  fluctuating  components  of  the  streamwise  and  trans¬ 
verse  velocities,  u'v*  is  the  Reynolds  stress  associated  with  the  disturbances,  and  u 
is  the  mean  streamwise  velocity,  has  the  same  sign  as  the  mean  centerline  velocity 
gradient.  Therefore,  when  the  mean  centerline  velocity  increases  the  fluctuation 
energy  should  increase,  and  when  the  mean  centerline  velocity  decreaises  the  fluctu¬ 
ation  energy  should  decrease.  Also  in  this  region,  Sato  and  Kuriki  (1961)  observed 
velocity  “over-shoots”  near  the  outer  edge  of  the  wake  so  that  the  mean  streamwise 
velocity  attained  values  greater  than  its  freestream  value. 

Motivated  by  the  appearance  of  vortices  in  wakes,  Sato  and  Kuriki  proposed 
a  “double  row  vortex  model”  as  a  way  of  explaining  the  nonlinesir  behavior  of 
wake  disturbances.  In  this  model,  the  vortices  alternate  in  sign  and  are  arranged 
in  a  staggered  fashion  on  each  side  of  the  wake  centerline.  The  velocity  induced 
by  this  vortex  pattern  does  indeed  capture  some  of  the  featxires  observed  in  wakes 
such  as  the  second  harmonic  near  the  wake  centerline,  the  velocity  over-shoots 
near  the  outer  edge  of  the  wake,  and  the  decrease  of  the  mean  centerline  velocity 
far  downstream. 

Several  diflerent  numerical  investigations  support  the  double  row  vortex  model. 
Zabusky  and  Deem  (1971)  performed  Navier-Stokes  simulations  of  the  temporal 
evolution  of  an  unstable  disturbEuice  in  a  parallel  wake.  Their  results  indicate 
that  the  saturation  of  the  fundamental  disturbance  in  plane  wakes  leztds  to  the 
development  of  a  vortex  street  pattern.  Phenomena  sinoilar  to  those  observed  by 
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Sato  and  Kuriki,  such  as  the  velocity  over-shoots,  were  observed  in  the  numerical 
results.  Aref  and  Siggia  (1981)  and  Meiburg  and  Lasheras  (1988)  studied  the  tem¬ 
poral  development  of  vortex  rows  using  invisdd  vortex  dynamics.  They  modeled 
a  plzme  wake  as  two  vortex  sheets  of  opposite  sign.  When  the  vortex  sheets  were 
perturbed,  they  evolved  into  a  two-dimensional  vortex  street. 

Several  researchers  including  Sato  and  Onda  (1970),  Sato  (1970),  Motohashi 
(1979)  and  Miksad  et  al.  (1982)  have  experimentally  investigated  the  behavior  of 
wakes  that  were  disturbed  at  two  different  frequencies.  By  perturbing  the  wake 
with  disturbances  at  two  different  frequencies,  these  researchers  hoped  to  gain 
insight  into  the  complex  nonlinear  interactions  that  take  place  in  naturzd  transi¬ 
tion.  Sato  and  Saito  (1975)  performed  experiments  to  study  the  interaction  of  a 
discrete  frequency  disturbance  with  the  broad  band  background  fluctuations  that 
were  naturally  present  in  the  flow.  In  their  experiments,  Sato  (1970)  and  Sato  and 
Onda  (1970)  foimd  that  the  nonlinear  interactions  between  the  two  disturbances 
generated  the  expected  higher  hrirmonics  as  well  as  a  difference  frequency.  This 
difference  frequency  is  similar  to  the  low  frequency  fluctuations  observed  in  natu¬ 
ral  transition.  In  natural  transition,  the  low  frequency  component  is  not  a  single 
frequency  but  a  narrow  band  of  frequencies.  The  conjecture  is  that  these  low  fre¬ 
quencies  are  generated  by  interactions  between  the  higher  frequency  fluctuations 
that  result  from  the  linear  instability. 

Miksad  et  al.  (1982)  fotmd  that  interactions  between  the  difference  frequency 
and  the  two  fundamental  frequencies  produced  sidebands  in  the  spectrum.  These 
sidebands  modulate  the  amplitude  of  the  fundamental  fluctuations.  Amplitude 
modulation  in  combination  with  the  dispersion  relation  of  the  fluctuations  pro¬ 
duces  phase  modulation  of  the  fluctuation-;.  Miksad  et  al.  indicate  that  this  mod- 
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ulation  has  an  important  role  in  spectral  broadening,  that  is  in  the  redistribution 
of  energy  from  a  small  number  of  frequencies  to  a  larger  number  of  frequencies. 
1.3  The  Three-Dimensional  Reeion  of  Transition 

Investigations  of  launinar-turbulent  transition  in  plane  waJces  have  dealt  mainly 
with  two-dimensional  instability  mechanisms.  Notable  exceptions  are  the  numeri¬ 
cal  work  of  Meiburg  and  Lasheras  (1987  and  1988),  the  experimental  investigations 
of  Breidenthal  (1980),  and  the  theoretical  work  of  Robinson  and  Saffman  (1982). 

Robinson  and  Saffman  (1982)  theoretically  investigated  the  stability  of  vortex 
arrays  to  small  amplitude  two-  and  three-dimensional  disturbances.  They  found 
that  the  staggered  pattern  of  the  Karman  vortex  street  is  most  unstable  with 
respect  to  three-dimensional,  spanwise-periodic  disturbances.  This  implies  that 
the  breakdown  of  the  vortex  streets  observed  in  flat  plate  wakes  should  be  of  a 
three-dimensional  nature. 

Breidenthal  (1980)  experimentally  investigated  the  three-dimensional  behav¬ 
ior  of  a  turbulent  flat  plate  wake.  The  plate  trailing  edge  varied  periodically  in 
the  spanwise  direction  in  order  to  induce  three-dimensional  disturbances  in  the 
wake.  Under  the  influence  of  this  three-dimensional  perturbation,  the  wake  was 
observed  to  form  a  spanwise  periodic  array  of  interconnecting  vortex  loops. 

Meiburg  and  Lasheras  (1987  and  1988)  numerically  investigated  the  interac¬ 
tion  of  a  large  amplitude,  three-dimensional  disturbance  with  a  smaller  amplitude, 
two-dimensional  disturbance.  They  found  that  the  two-dimensional  disturbance 
led  to  the  formation  of  a  vortex  street.  The  three-dimensional  disturbance,  through 
an  interaction  with  the  two-dimensional  vortex  street,  formed  lambda  vortices  sim¬ 
ilar  to  those  found  in  boundary  layers.  Further  interaction  between  the  lambda 
vortices  and  the  two-dimensional  structures  led  to  the  development  of  closed  vortex 
loops  similar  to  those  observed  by  Breidenthal. 
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CHAPTER  2 

PROBLEM  STATEMENT 

In  this  work,  laminsu- turbulent  tr&nsition  of  an  incompressible  high-deficit 
fiat  plate  waJce  is  investigated  by  direct  numerical  simulations  using  the  com¬ 
plete  Navier-Stokes  equations.  Research  efforts  are  initially  focused  on  the  early 
stages  of  transition  where  two-dimensional  disturbances  are  dominant.  Of  pri¬ 
mary  interest  is  the  role  of  nonlinear  effects  in  the  development  of  disturbances. 
Subsequently,  the  focus  is  on  the  early  stages  of  three-dimensional  breakdown 
where  both  two-  and  three-dimensional  disturbances  play  an  important  role.  The 
onset  of  this  three-dimensionality  is  investigated  by  simulating  the  interaction  of 
two-dimensional  and  three-dimensional  disturbances.  Investigations  of  these  in¬ 
teractions  might  explain  how  three-dimensional  disturbamces,  which  according  to 
linear  stability  theory  are  more  stable  than  two-dimensional  disturbances,  eventu¬ 
ally  dominate  the  wake  development. 

The  numerical  simulations  are  designed  to  model  the  physicaJ  experiments  of 
Sato  and  Kuriki  (1961)  and  Sato  (1970).  The  results  of  these  simulations  may  then 
be  compared  to  the  results  of  the  physical  experiments.  The  wake  that  is  studied 
is  shown  in  Figure  2.1.  The  wake  is  generated  by  a  thin  flat  plate  aligned  parallel 
to  a  uniform  stream.  The  boundary  layers  on  the  plate  are  assumed  to  be  laminar 
and  steady  over  the  entire  length  of  the  plate,  and  hence  form  a  laminar  wake 
when  they  merge  at  the  plate  trailing  edge.  The  wake  evolves  rapidly  downstream 
of  the  plate. 

The  spatial  domain  in  which  the  numerical  simulations  take  place  lies  down¬ 
stream  of  the  flat  plate  trailing  edge.  The  domain  is  placed  near  the  plate  so  that 
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it  will  contain  the  high-deficit  region  of  the  wake.  However,  the  domain  does  not 
include  the  trailing  edge  of  the  plate.  The  plate  trailing  edge  is  of  primary  interest 
in  studies  of  receptivity  and  in  studies  of  absolute  instabilities.  Although,  these 
are  both  important  aspects  of  wake  transition,  they  are  not  the  main  focus  of  this 
research. 

The  calculations  aire  undertaken  using  a  numerical  method  that  is  designed 
to  solve  both  the  two-dimensional  and  three-dimensional  Navier-Stokes  equations. 
Because  of  the  spanwise  periodicity  of  disturbances  that  was  observed  in  exper¬ 
iments,  the  spanwise  variation  of  the  wake  is  represented  by  finite  Fourier  se¬ 
ries.  Second-order  finite-differences  are  used  to  approximate  the  spatial  deriva¬ 
tives  in  the  streamwise  and  transverse  directions.  Time  integration  of  the  Navier- 
Stokes  equations  is  accomplished  using  a  hybrid  scheme  that  is  second-order  time- 
accurate. 

The  results  of  many  experimental  investigations  (Sato  and  Kuriki,  1961;  Sato, 
1970;  Mattingly  and  Criminale,  1972)  indicate  that  transition  in  wakes  is  the  result 
of  a  spatial  instability  rather  than  a  temporal  instability.  Furthermore,  Mattingly 
and  Criminale  (1972)  showed  that  the  use  of  the  Garter  transformation  (1962) 
to  convert  temporal  stability  results  into  spatiad  stability  results  is  invalid  for 
high-deficit  wakes.  Therefore,  in  this  work  spatial  amplification  and  temporal 
periodicity  of  wake  disturbances  is  assumed.  To  correctly  model  disturbances  of 
this  type,  the  boundary  conditions  employed  in  these  simulations  must  be  of  the 
inflow-outflow  type.  Spatially  growing  disturbances  are  induced  in  the  wake  by 
time-periodic  excitation  at  the  inflow  boundary  of  the  spatial  domain. 

The  remaining  chapters  describe  in  more  detail  other  important  aspects  of 
this  work.  In  Chapter  3,  the  governing  equations  of  these  simulations  are  dis¬ 
cussed.  The  length,  time,  and  velocity  scales  that  are  used  to  nondimensionalize 
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the  governing  equations  are  also  discussed.  The  boundary  and  initial  conditions 
that  are  required  to  complete  the  mathematical  formulation  of  the  physical  prob¬ 
lem  are  outlined  in  Chapter  4.  The  method  used  to  obtain  numerical  solutions  of 
the  governing  equations  is  discussed  in  Chapter  5.  Chapter  6  contains  a  discus¬ 
sion  of  the  results  of  the  numerical  simulations.  Finally,  in  Chapter  7  significsint 
conclusions  are  presented. 
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CHAPTER  3 

GOVERNING  EQUATIONS 


The  Navier-Stokes  equations  govern  momentum  transfer  in  a  viscous  fluid. 
These  equations  are  derived  by  applying  Newton’s  second  law  of  motion  to  a  flmd 
particle.  The  assumptions  required  to  obtain  these  equations  are  that  the  fluid  is 
a  continuum  and  that  its  stress-strain  relationship  is  Newtonian. 

3.1  The  Navier-Stokes  Equations 

The  numerical  method  is  based  on  the  Navier-Stokes  equations  for  incom¬ 
pressible  flow.  In  vector  form  they  can  be  written  as 

~ -I- u  •  Vu  =  — ivp -f  (3.1) 

at  p 

The  continuity  equation  for  incompressible  flow  is 


V.u  =  0  . 


(3,2) 


The  independent  variables  in  equations  (3.1)  and  (3.2)  Me  time  t  and  the 
spatial  position  vector  x  =  (i,y,z).  The  dependent  variables  are  the  Eulerian 
velocity  vector  field  u  =  (u,v,vj)  and  the  fluid  pressure  p.  These  variables  depend 
on  X  and  t.  The  three  coordinate  directions  x,  y,  ^md  z  and  the  corresponding 
fluid  velocity  components  u,  u,  and  w  are  shown  in  Figure  2.1.  The  fluid  density 
p  and  the  kinematic  viscosity  of  the  fluid  F  are  assumed  to  be  constant.  Variables 
with  over-bars  denote  dimensional  quantities,  while  bold-face  variables  represent 
vector  quantities. 

The  operator  V  is  a  vector  which  in  Cartesian  coordinates  is  of  the  form 


(I.  A  I.') 

V^x’  dy'  dz) 


(3.3) 
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The  Laplacian  operator  V  is 

=2  _  ^ 
^  ~  efx^'^  af'^ 

and  the  operator  u  •  V  is 


(3.4a) 


u  •  V  = 


-Q  -d 

tA— -  +  V—  +  ly— 

ax  ay  az 


(3.46) 


Equations  (3.1)  together  with  (3.2)  lead  to  four  scalzur  equations  which  in 
principle  can  be  numerically  solved  for  the  four  dependent  variables  (u,v,w)  and  p 
as  functions  of  (z,  y,  2)  and  1.  However,  one  difficulty  with  solving  these  equations 
is  formulating  boundary  conditions  for  the  pressure.  This  difficulty  is  avoided 
by  solving  an  alternate  formulation  of  the  Navier-Stokes  equations  in  which  the 
pressure  does  not  appear.  This  formulation  of  the  Navier-Stokes  equations  is 
obtained  by  taking  the  curl  of  equation  (3.1)  and  simplifying  the  result  with  the 
appropriate  vector  identities.  The  resulting  vector  equation  (Batchelor,  1967)  is 

—  —  ItV  u>  =  w  •  Vu  —  u  •  V5>  (3.5) 

at 


The  Eulerian  vorticity  vector  field  «•>  =  (u;j,u;j,,u>i)  is  defined  by  the  relation 

=  —V  X  u  (3.6) 


and  the  operator  or  •  V  is 


(3.7) 


It  is  noted  that  the  vorticity  vector  as  defined  by  equation  (3.6)  is  the  negative  of 
the  usual  definition  of  vorticity,  <7  =  V  x  u.  However,  throughout  this  work  the 
term  vorticity  will  be  used  to  denote  the  flow  variable  defined  by  equation  (3.6). 


31 


Furthennore,  because  of  the  negative  sign  in  equation  (3.6),  the  direction  of  fluid 
rotation  associated  with  the  vorticity  vector  ZS  =  —V  x  u  is  given  by  a  left-hand 
rule  instead  of  the  usual  right-hand  rule. 

Equation  (3.5)  is  the  Navier-Stokes  equation  in  vorticity-velocity  formulation, 
or  alternatively  the  vorticity  transport  equation.  Coupled  to  equation  (3.5)  is  the 
Poisson  equation 

V^u  =  V  X  Si  (3.8) 


which  is  derived  from  the  definition  of  vorticity  (Si  =  —V  x  u)  using  the  continuity 
equation  (3.2).  The  vorticity  Si  is  also  a  function  of  the  independent  variables  x 
and  1. 

Equations  (3.5)  and  (3.8)  are  a  system  of  six  scalar  equations  which  can  be 
solved  subject  to  appropriate  boundary  and  initial  conditions.  As  desired,  equa¬ 
tions  (3.5)  and  (3.8)  also  do  not  contain  the  fluid  pressure  p  ais  an  unknown  vari¬ 
able.  However,  in  contrzist  to  equations  (3.1)  and  (3.2),  now  six  scalar  equations 
must  be  solved  instead  of  four. 

3.2  Nondimensionalization 

Equations  (3.5)  through  (3.8)  are  nondimensionalized  using  appropriate  length, 
time,  and  velocity  scales.  The  independent  variables  are  nondimensionalized  using 
the  following  formulas: 


®  y  —  j  j  ^  ^oo 

X  =  •=  ,  y  =  =—  ,  z  —  z'y  ,  and  t  =  — = — 

i'^eu  I 


(3.9) 


where  the  variables  without  overbars  are  dimensionless.  The  streamwise  length 
scale  I  is  the  length  of  the  flat  plate.  In  the  transverse  direction  the  length  scale 
is  which  is  the  momentum  thickness  of  the  wake  at  the  trailing  edge  of  the  flat 
plate.  The  length  scale  that  is  used  in  the  spanwise  direction  is  1  /7  where  7  is  the 
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spanwise  wavenumber  of  the  three-dimensional  disturbance  that  results  from  the 
excitation  at  the  inflow  bovmdary.  The  time  scale  is  the  ratio  of  the  streamwise 
length  scale  i  to  the  freestream  velocity  17  oo* 

The  dependent  variables  u  =  (u^Vyw)  amd  S  =  are  nondimen- 

sionalized  as  follows: 


and 


u 

V 

w 

8 

It:) 

II 

Uoor2 

^~Uoorz  ’ 

Wydig 

WzSte 

■  Uoo 

,  Wz  =  = - 

C^ocra 

(3.10a) 


(3.106) 


The  parameter  r2  is  the  ratio  of  the  transverse  length  scale  6tg  to  the  streamwise 
length  scale  £:  r2  =  Similarly,  the  parameter  r3  is  the  ratio  of  the  spanwise 

length  scale  I/7  to  the  streamwise  length  scade  £:  rs  =  1/(7^). 

With  this  nondimensionalization,  equations  (3.5)  and  (3.8)  become 


-r - - V  t4>  =  f 

dt  Re 


(3.11a) 


and 


V^u  =  g  (3.116) 

which  are  solved  for  the  six  unknown  variables  u  =  (u,v,w)  and  =  (wjjWyjWi) 
as  functions  of  the  dimensionless  spatial  position  vector  x  =  (x,y,z)  and  dimen¬ 
sionless  time  t.  The  right  hand  side  of  equation  (3.11a),  the  vector  f  =  (/z,  /y,  /z), 
is  the  dimensionless  version  of  u>  •  Vu  —  u  •  Vw  which  appears  on  the  right  hand 
side  of  equation  (3.5).  The  three  components  of  f  are 


du  du  du  duix  d(A)x 

-  “IT  - -  “'■57 


dv 


dv 


dv  duv 


dij)x 


dUt 


(3.12a) 

(3.126) 


dw 


dw 


dw 


du>z  dwz 


dwz 


and  /,  =u>,—  +u;„— -u— -V— -xo— 


dx 


dy 


dz 


dx 


dy 


dz 


.  (3.12c) 
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The  right  hand  side  equation  (3.11b),  the  vector  g  =  is  the  dimension¬ 

less  version  of  V  x  ai  which  appears  in  equation  (3.8).  The  three  components  of  g 
are 


Ti  du>x  1  dfjfy 

r|  dy  Ti  dz  ^ 

1  dug  Ts  dug 

rjrl  dz  rl  dx  ’ 

,  1  duty  1  duz 

and  Qz  = - T - 2  “aT  * 

Tz  ax  Tzri  ay 

The  dimensionless  Laplacian  has  the  form 

2  32  1  02  1  52 

“  dx^  r|  ^2  +  ^2  Q^2 


(3.13a) 

(3.136) 

(3.13c) 

(3.14) 


and  the  parameter  Re  is  the  Reynolds  number  Re  = 

Equations  (3.11)  serve  as  the  basic  equations  for  the  simulations  described  in 
this  work.  Use  of  such  a  vorticity- velocity  formulation  for  transition  simulations 
was  suggested  by  Fasel  (1976)  and  used  successfully  for  both  two-dimensional 
(1976)  and  three-dimensional  simulations  of  transition  in  boundary  layers.  The 
same  formulation  was  also  successfully  employed  by  Pruett  (1986)  for  two-  and 
three-dimensional  transition  studies  in  free  shear  layers. 
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CHAPTER  4 

BOUNDARY  AND  INITIAL  CONDITIONS 

The  governing  equations  (3.11)  represent  a  set  of  partial  differential  equa¬ 
tions.  This  equation  system  is  parabolic  with  respect  to  the  time  variable  t  and 
is  elliptic  with  respect  to  the  spatial  variables  x,  y,  and  z.  Solution  of  equations 
(3.11)  requires  specification  of  suitable  spatial  and  temporal  domains.  Further¬ 
more,  the  unknown  variables  ta  and  u  must  be  specified  at  some  initial  time. 
Finally,  boundary  conditions  must  be  specified  along  the  entire  boundary  of  the 
spatial  domain. 

4.1  Spatial  Domain 

The  governing  equations  are  solved  in  a  three-dimensional  rectangular  region 
in  space  that  is  downstream  of  the  fiat  plate.  The  domain  extends  downstream 
from  near  the  trailing  edge  of  the  plate  and  is  assumed  to  extend  in  the  streamwise 
direction  for  several  disturbance  wavelengths.  The  transverse  thickness  of  the 
domain  is  assumed  to  be  much  wider  than  the  vorticity  disturbances  that  are 
present  in  the  wake.  The  domain  extends  into  the  irrotational  region  of  the  flow 
field  above  and  below  the  wake.  The  spanwise  extent  of  the  domain  is  equal  to 
one  nondimensional  spanwise  disturbance  wavelength  A,  =  27r. 

The  spatial  domain  and  coordinate  system  are  shown  in  Figure  2.1.  The  x  di¬ 
rection  corresponds  to  the  streamwise  direction  which  is  parallel  to  the  freestream 
velocity  U oo’  The  transverse  direction  is  denoted  by  the  independent  variable 
y.  The  third  coordinate  direction  z  corresponds  to  the  spanwise  direction  and  is 
parallel  to  the  trailing  edge  of  the  flat  plate.  For  convenience,  the  origin  of  the 
coordinate  system  is  placed  at  the  trailing  edge  of  the  flat  plate  .  The  spatial 
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domain  extends  in  the  streamwise  direction  from  xq  to  in  the  transverse  di¬ 
rection  from  I/O  to  i/M>  and  in  the  span  wise  direction  &om  z  =  0  to  z  =  27r.  The 
spatial  domain  is  then  given  by 


Xo  ^  X  <  Xff  f 

(4.1a) 

yo  <y<yM  , 

(4.16) 

and  0  <  z  <  2?r  . 

(4.1c) 

4.2  Temporal  Domain 

The  time  integration  of  the  governing  equations  is  divided  into  two  separate 
procedures  or  steps.  In  the  first  step,  the  governing  equations  are  solved  to  obtain 
the  steady  laminar  wake  that  appears  behind  the  flat  plate.  This  is  done  by 
integrating  the  governing  equations  (3.11)  with  respect  to  time,  subject  to  time- 
independent  boundary  conditions  (no  forcing),  until  a  steady  state  solution  is 
obtained.  This  calculation  begins  at  an  initial  time  t  =  <  0  for  which  initial 

values  of  the  dependent  variables  u(i,y,z,<£,j)  and  are  specified. 

The  calculation  ends  at  time  <  =  0,  at  which  time  the  steady  solution  is  obtained. 
This  steady  wake  is  referred  to  as  the  b^lse  flow. 

In  the  second  step  of  the  calculation,  the  response  of  the  wake  to  time- 
dependent  excitations  that  are  introduced  at  the  inflow  boundary  is  investigated. 
The  governing  equations  are  integrated  in  time  from  <  =  0  to  <  =  The  initial 
values  of  the  dependent  variables  at  time  t  =  0  are  given  by  the  base  flow. 

4.3  Boundary  and  Initial  Conditions  for  the  Undisturbed  Flow 

The  undisturbed  wake  is  assumed  to  be  two-dimensional,  so  that 


(4.2) 
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The  remaining  flow  variables  w(z,y,  z,<),  and  u>,(z,y,  z,t)  are  solu¬ 

tions  of  the  governing  equations  (3.11).  Restricted  to  two-dimensional  flows,  the 
governing  equations  become 


duiz  dtjJz  dutf  _  1  /  1  d^uiz 

di  dx  dy  RcV  dx^  t\  dy^  )  ’ 

di^u  1  df^u  Tz  duz 

dx^  r|  dy^  rl  dy  ’ 

d^v  1  d^v  _  Vi  d(j3z 

^  dx^  r\  dy^  r|  dx 


(4.3a) 

(4.3b; 

(4.3c) 


Because  the  base  flow  is  two-dimensional,  the  spatial  domain  reduces  from  a  three- 
dimensional  region  to  a  two-dimensional  plane  that  is  perpendicular  to  the  z  axis. 
Boundary  conditions  for  the  unknown  variables  u,  v,  and  a;,  need  to  be  specified 
zJong  the  boundaries  of  this  plane.  These  boundary  conditions  are  independent  of 
time  and  are  chosen  such  that  the  solution  of  equations  (4.3)  represents  a  laminar, 
high-deflcit,  flat  plate  wake. 

Inflow  Boundary  (z  =  zq) 

The  inflow  boundary  of  the  spatial  domain  is  located  downstream  of  the  flat 
plate  trailing  edge.  Therefore,  the  velocity  and  vorticity  distributions  at  this 
boundary  should  correspond  to  those  of  a  spatially-developing,  high-deficit  wake. 
For  most  of  the  calcxilations  in  this  work,  a  Gaussian  distribution  is  used  to  rep¬ 
resent  the  streamwise  velocity  at  the  inflow  boundary.  The  Gaussian  distribution 
compares  reasonably  well  to  experimentally  obtained  streamwise  velocity  profiles 
(Sato  and  Kuriki,  1961).  In  a  few  other  calculations,  a  hyperbolic  secant  function 
is  used  to  represent  the  inflow  streamwise  velocity.  The  hyperbolic  secant  dis¬ 
tribution  also  compares  reasonable  well  with  experimentally  obtained  streamwise 
velocity  profiles  (Mattingly  and  Criminale,  1972) 
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With  the  use  of  either  a  Gaussian  distribution  or  a  hyperbolic  secant  dis¬ 
tribution,  the  streamwise  velocity  is  then  known  as  a  function  of  y  at  the  inflow 
boundary.  However,  since  the  x  dependence  of  the  streamwise  velocity  is  unknown, 
being  given  by  the  solution  of  the  Navier-Stokes  equations,  the  transverse  velocity 
cannot  be  related  to  the  streamwise  velocity  using  the  continuity  equation  (3.2). 
Instead,  v  must  be  given  an  arbitrary  value.  In  this  work,  the  Reynolds  number  of 
the  wake  Re  =  is  large.  Therefore,  except  for  very  near  the  flat  plate  trailing 
edge,  the  transverse  velocity  is  much  smaller  thzin  the  streamwise  velocity  and  can 
be  set  to  zero  at  the  inflow  boundary  with  reasonable  accuracy. 

For  the  base  flow  calculation,  the  inflow  boundary  condition  is 


u{xQ,y,2,t)  =  u^{y)  , 

v(xo,y,z,t)  =  0  , 

and  a;,(zo,y,z,t)  =  wf(y) 


(4.4a) 

(4.46) 

(4.4c) 


where  the  superscript  (  )^  refers  to  the  fact  that  the  functions  u^{y)  and  u;^(y)  are 
the  streamwise  velocity  and  spanwise  vorticity  distributions  at  the  inflow  bound¬ 
ary.  If  the  streamwise  velocity  is  specified  to  have  a  Gaussiam  distribution,  then 
tf^(y)  is  given  by 

u\y)  =  1  -  (1  -  (4.5a) 


where 

ln(2) 


Alternatively,  u^{y)  can  be  given  by 


(4.55) 


» 


u\y)  =  1  -  (1  -  Ue)scch}{<ry) 


(4.6a) 
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where 


c 


arccosh('\/2) 

52 


The  spanwise  vorticity  is  given  by 


dy 


(4.66) 


(4.7) 


In  equations  (4.5)  and  (4.6),  Uc  is  the  centerline  velocity  of  the  wake  at  the  inflow 
boundary  and  6  is  the  wake  half-width  at  the  inflow  boundary.  The  parameters 
Uc  and  6  are  chosen  so  that  the  resulting  base  flow  closely  models  the  high-deflcit 
wakes  observed  in  experiments. 

Outflow  Boundary  (x  =  xn) 

For  spatially-developing  wakes,  a  proper  outflow  boundary  condition  is  not 
easily  found  since  the  flow  condition  at  this  location  would  be  obtained  by  the 
solution  of  the  governing  equations  and  is  therefore  not  known  o  priori.  However,  it 
is  possible  to  determine  an  outflow  boundary  condition  without  requiring  advance 
knowledge  of  the  solution.  Suppose  the  two-dimensional  version  of  equations  (3.5) 
and  (3.8)  are  nondimensionalized  using  the  boundary  layer  scaling 


=  =  ,  y  =  =i2e^  ,  u  =  ,  V  =  ,  and  w,  =  — .■  .  (4.8) 


The  parameters  t,  and  Re  are  defined  in  section  3.2.  The  resulting  dimen¬ 
sionless  equations,  for  steady  flow,  are 


du>. 

d<j)i 

1  d^ug  d^fjjg 

(4.9a) 

dy 

■  Re  dx^  dy^  ’ 

1 

d^u 

dWg 

(4.96) 

Re  dx^ 

+ 

’’W  ’ 

1  d^v 

d^v 

dUf 

(4.9e) 

Redx^ 

+ 

dx 

and 
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For  high  Reynolds  number  flows  JRe  >>  1,  as  in  these  investigations,  the  stream- 
wise  diffusion  term  is  the  smallest  term  in  each  of  the  three  equations  (4.9a) 
through  (4.9c).  Therefore,  the  streamwise  diffusion  terms  are  set  equal  to  zero, 
resulting  in  the  outflow  boundary  conditions 


d^u 

(4.10a) 

d^v 

^(*iv,y,2,t)  =  0  , 

(4.106) 

-^^(*iv,y,i,t)  =  0  . 

(4.10c) 

These  boundary  conditions  were  first  proposed  by  Fasel  (1976). 

Freestream  Boundaries  (y  =  j/o)  and  (y  =  yw) 

As  stated  previously,  the  base  flow  is  assumed  to  be  a  high  Reynolds  number 
flow.  Therefore,  the  rotational  part  of  the  wake  is  confined  to  a  region  very  near 
to  its  centerline  and  the  streamwise  velocity  rapidly  approachs  its  freestream  value 
as  y  -+  ±oo.  The  boundary  conditions  that  satisfy  these  assumptions  are,  for  the 
lower  freestream  boundary, 


«(®iyo,^,0  =  i  . 

(4.11a) 

dv 

^(®,yo,*,0  =  o  > 

(4.116) 

and  w*(x,yo,«,  t)  =  0 

(4.11c) 

Similarly,  for  the  upper  freestream  boundary: 

=  1  , 

and  w,(x,yM,*,0  =  0  • 


(4.12a) 

(4.12fc) 

(4.12c) 
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Equations  (4.11b)  and  (4.12b)  follow  directly  from  the  conditions  (4.11a)  and 
(4.12a)  using  the  continuity  equation. 

Auxiliary  Condition  (y  =  0) 

The  boundary  conditions  described  up  to  this  point  are  sufficient  to  ensure 
the  uniqueness  of  the  solution  of  the  vorticity  equation  (4.3a)  and  the  Poisson 
equation  for  the  streamwise  velocity,  equation  (4.3b).  However,  these  boundary 
conditions  are  not  sufficient  to  ensure  the  uniqueness  of  solutions  of  the  transverse 
velocity  equation  (4.3c).  This  was  shown  by  Pruett  (1986)  who  notes  that  if  a 
function  V (x,  y)  is  a  solution  of  equation  (4.3c),  subject  to  the  boundary  conditions 
described  previously  in  this  section,  then  the  function 

+  c(x  -  zo)  (4.13) 

where  c  is  a  constant,  is  also  a  solution  of  equation  (4.3c).  Therefore,  an  infinite 
number  of  solutions  satisfy  equation  (4.3c)  and  the  boundary  conditions  described 
in  this  section. 

Therefore,  an  additional  condition  has  to  be  specified  that  will  ensure  that 
the  solution  of  equation  (4.3c)  is  unique.  Experimental  evidence  indicates  that 
the  steady  streamwise  velocity  is  symmetric  with  respect  to  the  wake  centerline 
(y  =  0).  Therefore,  the  transverse  velocity  v  is  antisymmetric  with  respect  to  the 
wake  centerline  (y  =  0)  and  is  equal  to  zero  at  the  centerline.  The  transverse 
velocity  v  is  therefore  required  to  satisfy  the  condition 

r(z,0,z,<)  =  0  (4.14) 


at  the  wake  centerline. 
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Initial  Conditions  (t  = 

This  calculation  starts  at  time  t  =  for  which  initial  valuer  of  the  unknown 
variables  u,  v,  and  are  specified.  In  this  work,  the  initial  conditions  are  chosen 
such  that  the  flow  profile  is  independent  of  x  and  is  identical  to  that  of  the  inflow 
boundary.  These  conditions  are 


=  «^(y)  , 

(4.15a) 

f(2:,y,2,fL,)  =  0  , 

(4.156) 

<^zix,y,z,iL^)  =  u[{y)  . 

(4.15c) 

4.4  Boundary  and  Initid  Conditions  for  the  Disturbed  Flow 

For  calculations  of  the  disturbed  flow,  which  may  be  three-dimensional,  bound¬ 
ary  conditions  have  to  be  specified  for  all  boundaries  of  the  three-dimensional 
domain. 

Inflow  Boundary  (x  =  xq) 

In  a  physicjil  experiment,  disturbances  in  wakes  may  originate  from  a  variety 
of  sources.  Wake  disturbances  may  arise  from  freestream  turbulence,  roughness 
on  the  flat  plate,  sound,  etc.  Disturbances  can  also  be  artificially  introduced  into 
wakes  through  loudspeakers,  flaps  at  the  plate  trailing  edge,  or  through  heater 
strips  on  the  plate.  In  this  work,  the  disturbances  are  introduced  at  the  inflow 
boundary.  For  this  purpose,  solutions  of  the  Orr-Sommerfeld  equation  are  used. 

Thus,  the  velocity  and  vorticity  components  at  the  inflow  boundary  are  a 
combination  of  those  of  the  steady  base  flow  discussed  in  the  previous  section  and 
of  time  varying  components  that  then  excite  the  disturbances  in  the  wake.  The 
inflow  boundary  conditions  are  of  the  form 


u(xo , y,  z,  t)  =  «ss(®0, y)  +  P^'^ixo, y,  0  +  t)cos(z)  ,  (4.16a) 
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»(xo,y,2iO  =  0 

+  P^‘'(xo,y,t)  +  P*‘'(xo,y,t)cos(z)  , 

(4.166) 

u>(xo,y,z,<)  =  0 

+  Pd‘'(xo,y,0co8(z)  , 

(4.16c) 

w*(xo,y,z,t)  =  0 

+  P“(*o,y,0cos(^)  . 

(4.16d) 

Wy(®o,y,x,0  =  0 

+  ^“(®o»yi<)cos(z)  , 

(4.16e) 

and  a;,(io,y,«,0  =  ‘^*5s(®o,y)  +  P“(®o,y,<)  +  P^f(*o,y,*)cos(z)  .(4.16/) 

The  functions  uss(®0iy)  and  ‘*'*ss(®o»y)  are  the  streamwise  velocity  and 
spanwise  vorticity  components  of  the  base  flow  at  the  inflow  bouc^a^y.  The  func¬ 
tions  P^‘^(io,yiO»  and  P^'j{xo,y,t)  are  two-dimensional  disturbances 

obtained  from  the  Orr-Sommerfeld  equation.  These  functions  are 

Pu'^{xQ,y,  t)  =r{t)A2U^/{y)cos{al'^xo  -  /?<  +  ^u'^Cy))  ,  (4.17a) 

-P»‘*(*o,y,0  =r(t)>l2t>4'(y)cos(aj^ro  -/?<  +  <f>l'^{y))  ,  (4.176) 

and  Pi‘‘(xo,y,t)  =r(t)yl2u;2^(y)cos(o“io -/3t-t-(A^^(y))  .  (4.17c) 

The  functions  P’‘'(xo,y,0>  Pv‘^(®o,y,<),  P®'^(xo,y,t),  P’^(xo,y,t),  Pw^(*o,y,0. 
andP^;'(io,y,t)  are  three-dimensional  disturbances  obtained  from  the  three-dimen¬ 
sional  Orr-Sommerfeld  equation.  These  functions  are 

P«'^(®o,  y,  0  =r(<)>l3Mx‘*(y)cos(a’'^io  -  0t  +  <^“(y))  ,  (4.18a) 

P^‘*(®o, y,  0  =r(/)>l3t;^‘'(y)cos(a®'*xo  -Pt  +  <f>V^{y))  ,  (4.186) 

Pw‘*(®o,  y,  0  =r{t)Azw^/{y)cos{al^xo  -  /?<  -f-  (i>^J{y))  ,  (4.18c) 

y,  0  =r(<)i4su;’^(y)cos(a’‘'®o  -  /?<  +  (y))  ,  (4.18d) 

Ptti^o,y,i)  =r(<)yl3u;J^(y)cos(a*''io  - /?<  +  <f>l^^{y))  ,  (4.18c) 

and  P’f (xo,  y,  t)  =r(t)i43u;’^(y)cos(a’‘'xo  -  (y))  •  (4.18/) 
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For  the  three-dimensional  disturbances,  a  superpositidn  of  two  oblique  waves  of 
equal  ajid  opposite  spanwise  waveniimbers  is  used. 

In  the  preceding  paragraph,  the  parameter  was  the  reed  part  of  the  eigen¬ 
value  of  the  two-dimensional  Orr-Sommerfeld  equation  corresponding  to  the  fre¬ 
quency  /3.  The  parameter  was  the  real  part  of  the  eigenvalue  of  the  three- 
dimensional  Orr-Sommerfeld  equation  corresponding  to  the  frequency  ^  and  the 
nondimensional  spanwise  wavenumber  7  =  1.  The  functions 

the  amplitudes  and  phases  of  the  eigen¬ 
function  of  the  two-dimensional  Orr-Sommerfeld  equation  corresponding  to  the 
frequency  The  functions  {ul^{y),<i>i‘^{y)),  (y)» (y)- 

{^zlt{y)y^ti{y))^  {^y^Aiy)y‘i>li(y))y  (y)>'^L1(y))  are  the  amplitudes  and 

phases  of  the  eigenfunction  of  the  three-dimensional  Orr-Sommerfeld  equation 
corresponding  to  the  frequency  0  and  the  nondimensional  spanwise  wavenumber 
7  =  1.  Both  the  two-dimensional  and  the  three-dimensional  eigensolutions  cor¬ 
respond  to  the  streamwise  location  zo  at  which  the  steady  streamwise  velocity 
distribution  is  v.ss{^(iiy)-  The  amplitude  of  the  disturbances  at  the  inflow  bound¬ 
ary  is  controlled  by  the  parameters  A2  and  >13.  The  solution  procedure  used  to 
obtain  eigensolutions  of  the  Orr-Sommerfeld  equation  is  discussed  in  Appendix  A. 

The  function  r(t)  that  appears  in  equations  (4.17)  and  (4.18)  is  chosen  to  min¬ 
imize  any  transient  disturbance  that  might  zudse  when  the  excitation  is  initiated. 
The  function  r(t)  has  the  values 

Outflow  Boundztfv  (z  =  zjv) 

For  numerical  simulations  of  spatially-amplifying  disturbances  in  shear  flows, 
the  specification  of  appropriate  outflow  boundary  conditions  represents  a  major 
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difEculty.  Niimerical  simulations  of  spatially-periodic,  temporally  amplifying  dis¬ 
turbances  avoid  this  problem  because  then  periodic  inflow-outflow  boundary  con¬ 
ditions  can  be  employed. 

For  the  numerical  simulation  of  two-dimensional,  spatiaJly-amplifying,  small- 
amplitude  disturbances  in  a  flat  plate  boundary  layer,  Fasel  (1976)  employed  the 
condition 


£V  _ 
5x2 


u 


(4.20) 


at  the  outflow  boundary.  In  equation  (4.20),  u'  =  u  —  ttss  is  the  streamwise  per¬ 
turbation  velocity  and  is  the  real  part  of  the  disturbance  wavenumber.  Similar 
conditions  were  used  for  the  other  flow  veiriables.  In  Fasel’s  numerical  simulations, 
the  boundary  condition  (4.20)  allowed  small-amplitude  disturbances  to  pass  unre¬ 
flected  through  the  outflow  boundary.  However,  equation  (4.20)  is  formally  exact 
only  for  small  amplitude  neutrally  stable  disturbances.  Therefore,  while  equa¬ 
tion  (4.20)  is  a  suitable  approximation  for  the  slowly  amplifying  disturbances  that 
are  observed  in  flat  plate  boundary  layers,  it  is  definitely  not  appropriate  when 
disturbances  are  highly  amplified,  as  in  high-deficit  wakes. 

Pruett  (1986)  attempted  to  use  outflow  boundary  conditions  that  were  sim¬ 
ilar  to  equation  (4.20)  for  the  simulation  of  spatially-amplifying  disturbances  in 
free  shear  layers.  However,  as  Pruett  notes,  amplification  rates  are  quite  large  in 
free  shear  layers,  and  therefore  the  previous  condition  is  not  valid.  Pruett  foxmd 
that  boundary  conditions  similar  to  equation  (4.20)  allowed  the  fundamental  dis¬ 
turbance  component  to  pass  through  the  boimdary  unreflected,  but  reflected  any 
nonlinearly  generated  harmonics  of  the  fundamentaJ  fluctuation. 

As  an  alternative  to  boundary  conditions  such  as  equation  (4.20),  Pruett 
employed  a  moving  outflow  boundary  for  his  calculations.  The  moving  outflow 
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boundairy  is  implemented  in  the  following  manner.  First,  the  base  flow  is  com¬ 
puted  in  the  entire  spatial  domain  for  which  xq  <  x  <  xj\r.  Then,  the  disturbed 
flow  is  calculated  in  a  smaller  subset  of  the  spatial  domain  bounded  by  z  =  xq  and 
X  =  XMB(t)  <  where  x  =  XAfs(t)  is  the  instantaneous  location  of  the  moving 
boimdary.  The  moving  boundary  is  initially  located  at  a  specified  distance  down¬ 
stream  of  the  inflow  boundary.  As  the  calculation  progresses  and  the  disturbances 
propagate  downstream,  the  moving  boundary  propagates  downstream  ahead  of 
the  disturbances  so  that  the  flow  remains  undisturbed  at  this  boundary.  As  a 
consequence,  Dirichlet  boundary  conditions  that  enforce  an  undisturbed  flow  can 
be  employed  at  the  location  of  the  moving  boundary.  When  the  moving  boundary 
reaches  the  fixed  outflow  boundary  at  x  =  xjv,  the  c2dculation  is  stopped. 

For  the  present  simulations,  boundary  conditions  similar  to  those  used  by 
Pruett  are  employed.  However,  instead  of  using  a  moving  outflow  boundary,  the 
undisturbed  flow  condition  is  enforced  at  the  outflow  boundary  located  at  i  =  Xjv- 
The  disturbed  flow  is  then  calculated  in  the  entire  spatial  domain.  Furthermore, 
the  domain  is  specified  to  be  large  enough  so  that  the  disturbances  do  not  reach 
the  outflow  boundary  before  the  calculation  is  stopped.  The  resulting  boundary 
conditions  are 


v(xN,t/,z,i)  =  uss(xN,y)  , 

(4.21a) 

v(xN,yyZ,t)  =  vssi^iN,!/)  , 

(4.216) 

w(x;^,y,z,i)  =  0  , 

(4.21c) 

w*(xiv,y,z,t)  =  0  , 

(4.21d) 

i^y(xN,yyZ,t)  =  0  , 

(4.21c) 

wz(®iv,y,x,0  =  ‘*'*ss(®N,y)  • 

(4.21/) 

As  before,  the  subscript  (  )55  denotes  the  steady  state  wake  or  base  flow. 
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Freestream  Boundaries  (y  =  yo)  and  (y  =  y^f) 

The  freestream  boundary  conditions  are  based  on  the  assumption  that  distur¬ 
bances  at  these  locations  decay  exponentially  at  rates  predicted  by  linear  stability 
theory.  These  boundary  conditions  are  now  derived. 

An  equation  for  the  streamwise  perturbation  velocity  u'  =  u  —  uss  c^n  be 
obtained  by  subtracting  the  base  flow  streamwise  velocity  equation  (4.3b)  from 
the  streamwise  component  of  equation  (3.11b).  The  resulting  equation  is 

d^u'  1  d^u'  1  d^u‘  Ti  duz'  1  duy 

where  the  spanwise  perturbation  vorticity  u>z  —  tt>xss  is  denoted  by  a;^'.  Assum¬ 
ing  that  the  velocity  and  vorticity  fluctuations  at  the  freestream  boundaries  are 
described  by  linear  stability  theory,  then 

u'(i,y,z,t)  =  Real(u(y)e‘(“*+*-^‘))  (4.23) 

and  analogous  formulas  are  valid  for  the  other  velocity  and  vorticity  components. 
Then,  from  Squire’s  equation 


<^y"  —  +  1  +  iaRe  (U  —  c))  uiy  =  —iReU'v  (4.24) 


it  can  be  shown  that  for  Re  »  1,  Jy  ~  0  for  |y)  >>  1.  Furthermore,  from 
the  continuity  equation  and  the  definition  of  and  w,,  it  can  also  be  shown 
that  for  velocity  fluctuations  that  decay  exponentially  when  |y|  >>  1,  then  w*  ~ 
0  and  tJ,  ~  0  for  |y|  >>  1.  Combining  equation  (4.22)  with  equation  (4.23), 
and  assuming  that  the  vorticity  fluctuations  are  negligible  for  large  |y|,  then  the 
equation 


(y)  =  rl{al  -|-  f^)ti(y) 


(4.25) 


47 


is  valid  if  the  disturbance  is  assumed  to  neutrally  stable  (oj  =  0).  In  reality 
however,  wake  disturbances  are  usually  highly  amplified  and  therefore  equation 
(4.25)  is  not  strictly  applicable.  However,  for  the  calculations  presented  in  this 
work,  the  assumption  of  neutrally  stable  disturbances  at  the  freestrezim  boundaries 
did  not  prove  detrimental.  Equation  (4.25)  has  the  solution 


u{y)  =  Ae  y  ^  +  Be  V  . 


(4.26) 


Assuming  that  u{y)  is  finite  for  all  y,  then  for  y  — »  — oo,  B  =  0  and  u{y)  satisfies 
the  difi'erential  equation 


^=,Ja;  +  i4(v)  . 


(4.27) 


Combining  equation  (4.27)  with  equation  (4.23)  results  in  the  following  boundary 
condition  for  the  streamwise  perturbation  velocity  u'(i,y,  z,f)  at  y  =  yo: 


I 


(4.28) 


Similarly,  the  boundary  condition  at  y  =  yM  is 


I  1. 

=  -’■ijal  +  ■^u'{i,yM, z,t)  . 


(4.29) 


In  a  simil2U’  manner,  the  following  boundary  conditions  for  the  spanwise  ve¬ 
locity  w  are  obtained: 


dw ,  .  /  •»  1  /  \ 

—  (i,yo,z,<)  =  r2^jal  +  ^u;(x,yo,  z,f) 


and 


(4.30a) 
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Boundary  conditions  for  the  transverse  velocity  v  are  obtained  from  the  con¬ 
tinuity  equation.  These  boundary  conditions  are 


and 


dv ,  .  du  dw ,  . 

j/o, 


(4.31a) 


(4.316) 


The  vorticity  at  both  freestream  boundaries  is  assumed  to  be  zero.  At  the 
lower  boundary  (y  =  yo) 


w*(x,yo,z,<)  =  0  , 

Wj(x,yo,z,<)  =  0  , 

w,(x,yo,2,<)  =  0 


(4.32a) 

(4.326) 

(4.32c) 


and  at  the  upper  boundary  (y  =  y^ ) 


w*(x,yA/,z,<)  =  0  »  (4.33a) 

=  0  ,  (4.336) 

w*(®,yM,«»0  =  0  •  (4.33c) 

One  further  consideration  is  the  form  of  the  coefficient  r2  in  equa¬ 

tions  (4.28),  (4.29),  and  (4.30).  These  boundary  conditions  were  derived  based 
on  the  sissumption  of  three-dimensional  disturbances  for  which  Or  —  and  for 
which  the  coefficient  rj  yjal  -1-  ^  becomes  rj  -I-  For.  two-dimensional 

disturbances,  the  coefficient  -f  ^  becomes  simply  When  both  two- 

dimensional  and  three-dimensional  disturbances  are  present  in  the  walce  it  is  not 


49 


obvious  whether  r2  +  71  should  be  set  equal  to  aj‘*r2  or  This 

difficulty  is  resolved  by  decomposing  the  boimdary  conditions  (4.28),  (4.29),  and 
(4.30)  into  two-dimensional  and  three-dimensional  components.  Then,  for  the  two- 
dimensional  component  of  equations  (4.28),  (4.29),  and  (4.30),  r2y'aj  -f-  ^  is  set 
equal  to  a$‘*r2.  For  the  three-dimensional  component  of  equations  (4.28),  (4.29), 
and  (4.30),  r2  ^  is  set  equal  to  4-  This  decomposition  will  be 

discussed  further  in  the  next  chapter. 

Lateral  Boundaries  (z  =  0)  and  (z  =  27r) 

As  discussed  in  section  4.1,  the  spanwise  width  of  the  spatizd  domain  is  equal 
to  one  spanwise  wavelength  =  27r.  Therefore,  periodic  lateral  boundary  condi¬ 
tions  are  employed.  These  boundary  conditions  for  the  u  component  are 

u(i,y,0,t)  =  u(i,i/,27r,t)  ,  (4.34) 

with  analogous  conditions  for  all  other  variables. 

Initial  Conditions  {i  —  0) 

This  calculation  starts  at  time  <  =  0.  At  this  time  an  undisturbed  flow 
is  cissumed  and  therefore  the  flow  variables  are  set  equal  to  the  vzdues  of  the 
previously  computed  base  flow.  Thus 


u(i,y,z,0)  =uss(*,y)  , 

(4.35a) 

v(i,y,z,0)  =  x;s5(*,y)  , 

(4.356) 

u;(i,y,z,0)  =  0  , 

(4.35c) 

w,(®,y,a:,0)  =  0  , 

(4.35d) 

=  0  > 

(4.35c) 

w,(x,y,z,0)  =  w,55(i,y)  . 

(4.35/) 
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CHAPTER  5 

NUMERICAL  METHOD 

The  governing  equations  (3.11)  together  with  the  boundary  and  initial  condi¬ 
tions  described  in  the  previous  chapter  are  solved  using  a  combination  of  spectral 
and  finite-difference  methods.  Fourier  series  are  used  to  represent  the  spanwise 
variation  of  the  dependent  variables.  Spatial  derivatives  in  the  streamwise  and 
transverse  direction  are  approximated  using  second-order  finite-differences.  The 
vorticity  equation  (3.11a)  is  numerically  integrated  in  time  using  a  combination  of 
an  Alternating-Direction-Implicit  (ADI)  method,  the  Crank-Nicolson  method,  and 
the  second-order  Adams-Bashforth  method.  The  nonlinear  terms  are  computed 
spectrally  at  each  time  step.  The  discretized  Helmholtz  equations  are  solved  using 
fast  Helmholtz  Solvers.  The  numerical  method  discussed  here  is  based  on  the  work 
of  Pruett  (1986). 

5.1  Spanwise  Spectral  Approximation 

As  discussed  in  Chapter  4,  the  disturbance  flow  is  assumed  to  be  periodic 
in  the  spanwise  direction  with  wavelength  A,  =  27r.  Because  of  this  spanwise 
periodicity,  the  streeunwise  velocity  can  be  represented  by  a  truncated  Fourier 
series  of  the  form 

4 

u(x,y,z,<)~  ^  C/^*(x,y,<)c’**  (5.1) 

with  analogous  Fourier  representations  for  the  other  flow  variables.  The  Fourier 
coefficients  V*,  fl*,  fij,  and  fl*  are  functions  of  the  two  spatial  variables 
X  and  y  and  the  time  variable  t.  The  Fourier  coefficients  can  be  rewritten  as  the 
vectors 


u‘  =  (i^*,y*,iy*) 


I 


(5.2a) 
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and  =  (5.2a) 

which  are  the  Fourier  representations  of  the  velocity  vector  u  =  (u,  v,  w)  and  the 
vorticity  vector  ta  =  Furthermore,  because  the  physical  variables 

u  =  (ti,  v,to)  and  <if  =  are  real  numbers,  the  Fourier  coefficients  satisfy 

the  relations 


U‘(x,2/,t)  =  U  *(i,y,t) 

**  lb 

and  n*(x,y,<)  =  fi  (x,y,t) 


(5.3a) 

(5.36) 


where  "  denotes  the  conjugate  of  a  complex  number.  Equations  (5.3)  are  valid  for 
adl  values  of  k. 

In  the  governing  equations  (3.11),  the  physical  variables  are  replaced  by  their 
Fourier  representations  given  by  equation  (5.1).  Orthogonzdity  of  the  exponential 
functions  in  the  spanwise  interval  0  <  z  <  27r  decouples  the  three-dimensional 
governing  equations  into  y  +  1  sets  of  two-dimensional  equations  of  the  form 


at 

at 

an!; 

at 


1 

Re 

1 

Re 

1_ 

Re 


a^n^  1  a^n^ 

^  ax^  r|  ay^ 

'a^nj;  j_a^ 

^  ax^  r|  ay^ 

^  r\  ay^ 

a^u’‘  1  a^u’’ 

ax^  r|  dy^ 

a^v^  1 

**”  r|  ay^ 

a^w>‘  1  a^w'^ 

ax^  r*  ay^ 


(5.4a) 

(5.46) 

(5.4c) 

{5.4<i) 

(5.4e) 

(5.4/) 


Equations  (5.4)  are  solved  for  0  <  Jb  < 


^  to  obtain  the  complete  set  of  Fourier 


coefficients.  The  Fourier  coefficients  for  which  k  <  0  are  obtained  from  equations 


f 
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(5.3).  Once  the  Fourier  coefHcients  are  obtained  £rom  equations  (5.4),  the  physical 
variables  u  and  ut  can  be  reconstructed  using  equations  (5.1). 

The  functions  F*,  F*,  and  F*  are  Fourier  representations  of  the  nonlinear 
terms  /*,  /,,  and  /*.  These  fimctions  are  related  by  the  expressions 


(5.5a) 

‘4" 

1 

II 

.« 

F*(x,y,<)e‘*'  , 

(5.56) 

*=-4: 

K 

T 

and  f;i(x,y,2,t)  Fj'(x,y,/)e‘''*  . 

(5.5c) 

L_  K 
k - j- 

The  functions  /*,  /j,,  and  are  defined  by  equations  (3.12). 

The  functions  Gy,and  G\  are  Fourier  representations  of  Qy,  and  p,. 
These  functions  are  related  by  the  expressions 


9x(x,y,z,t)oi  ^  G‘(x,y,t)e*‘*  , 

(5.6a) 

4 

gy{x,y,z,t)oi  ^  Gj(x,y,t)c’*'  , 

(5.66) 

*=-4: 

4 

and  p,(x,y,z,<)~  ^  G^{x,y,t)e^^‘  . 

(5.6c) 

fc=-4 

The  functions  y*,  and  y,  are  defined  by  equations  (3.13). 

GJ,  and  G*  have  the  simple  form 

The  functions  G*, 

fc  r,  an*  i*  * 
^‘■r|ay  ’ 

{5.7a) 

*■*  r^fc 

»"r3r|"*  rl  dx  ’ 

(5.75) 

r^T2  '2 

rj  ax  rjrj  ay 


(5.7c) 
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Equations  (5.4)  are  solved  for  the  Fourier  coefticients,  subject  to  the  Fourier 
repr‘‘3entations  of  the  boundary  and  initial  conditions  from  Chapter  4.  For  the 
undisturbed  flow  calctilation,  Fourier  representations  of  the  boundary  and  initial 
conditions  &om  Section  4.3  axe  employed.  For  the  disturbed  flow  calculation, 
Fourier  representations  of  the  boundary  and  initial  conditions  from  Section  4.4 
are  employed. 

5.1.1  Fourier  Representation  of  the  Boundary  and  Initial  Conditions  for  the 
Undisturbed  Flow 

The  Fourier  representation  of  the  boundary  and  initial  conditions  for  the 
undisturbed  flow  are  discussed  in  this  section.  The  undisturbed  flow  is  two- 
dimensionad,  so  that 

C/*  =  U''  =  Pt"*  =  fi*  =  fij  =  n*  =  0  for  ib  >  0  .  (5.8) 

Furthermore,  for  i  =  0 


=  0  •  (5-9) 


The  remaining  Fourier  coefficients  U^{x,y,t),  V’®(i,y,t),  and  n°(i,y,  <)  satisfy 
the  equations 


an;  i  i 

di  Re\  az2  r|  ay2  y  *  ’ 

a^u”  1  a^u°  _  rs  an; 

ax^  Xj  dy^  t\  dy  ’ 

a^v°  1  a^u°  _  ra  an; 

dx^  Xj  dy^  Xj  dx 

which  are  the  Fourier  analogs  of  equations  (4.3)  in  Chapter  4.  The 
F°(x,y,  <)  in  equation  (5.10a)  is 


(5.10a) 

(5.10b) 

(5.10c) 


nonlinear  term 


^yo 


dy  ) 


(5.11) 
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Solution  of  equations  (5.10)  requires  the  Fourier  representation  of  the  boundary 
and  initial  conditions  from  Section  4.3. 

Inflow  Boundary  (x  =  iq) 

The  Fourier  representation  of  the  inflow  boundary  conditions  (4.4)  are 


(5.12a) 

(5.126) 

and  Q^(xo,y,i)  =  u;^(y)  . 

(5.12c) 

Outflow  Boundary  (x  =  xiv) 

The  Fourier  representation  of  the  outflow  boundary  conditions  (4.10) 

are 

(5.13a) 

^^2  =  0  » 

(5.136) 

and  -^~(xAr,y,t)  =  0  . 

(5.13c) 

Freestream  Boundaries  (v  =  unl  and  (v  =  vm) 

The  Fourier  representation  of  the  lower  freestream  boundary  conditions  are 

(5.14a) 

— (*,yo,0  =  0  > 

(5.146) 

and  fi®(z,yo,<)  =  0  . 

(5.14c) 

Similarly,  the  Fourier  representation  of  the  upper  freestream  boimdary  conditions 

are 

f^®(®,yM,f)  =  1  , 

(5.15a) 

dV° 

-^(®,yM,f)  =  o  , 

(5.156) 

and  fl®(x,yM,<)  =  0  • 

(5.15c) 
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Auxili2irv  Condition  fv  =  0) 

The  auxiliary  condition  is 

V®(z,0,t)  =  0  . 

(5.16) 

InitiaJ  Conditions  (i  =  ir. ) 

The  initial  conditions  for 

the  base  flow  calculation  are 

=  u^(y)  , 

(5.17a) 

=  0  , 

(5.176) 

and 

=  wj(y)  . 

(5.17c) 

5.1.2  Fourier  Representation  of  the  Boundary  and  Initial  Conditions  for  the 
Disturbed  Flow 

The  Fourier  representation  of  the  boundary  and  initial  conditions  that  aire 
used  for  the  calculation  of  the  disturbed  wake  are  discussed  in  this  section.  In  gen¬ 
eral,  the  disturbed  wake  is  three-dimensional.  This  requires  solving  the  governing 
equations  for  the  three  velocity  components  zuid  the  three  vorticity  components. 
Boundary  and  initial  conditions  are  required  for  all  of  these  components. 

Inflow  Boundary  (i  =  zq) 

The  Fourier  representation  of  the  inflow  boundary  conditions  are  shown  below. 
For  k  =  0,  the  boundary  conditions  are 


^°(xa,y,t)  =  «ss(®o,y)  +  Pl’^{xo,y,t)  , 

(5.I80) 

V°{xo,y,t)  =  0  + P^'^{xo,y,i)  , 

(5.186) 

W\xo,y,t)  =  0  , 

(5.18c) 

J^*(®o,y,t)  =  0  , 

(5.18d) 

f^5(*o,y,<)  =  0  , 

(5.18c) 

and  n°(xo,y,t)  =  u>,55(zo,y)  +  P3f(io,y,t)  • 

(5.18/) 
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For  k  =  1,  the  inflow  boundary  conditions  are 

with  similar  conditions  for  the  other  flow  variables.  For  k  >  1,  the  Fourier  com¬ 
ponents  of  the  velocity  amd  vorticity  at  the  inflow  boundary  are 

?7*(io,y,<)  =  0  etc.  (5.20) 

Outflow  Boundary  (i  =  zjv) 

At  the  outflow  boundary,  the  wake  is  assumed  to  be  undisturbed.  For  k  =  0, 
the  boundary  conditions  (4.21)  become 


=  f^ssi^Nty)  , 

(5.21a) 

^°(®iv,y,f)  =  uss(*N,y)  , 

(5.216) 

^°(®;v,y.f)  =  0  , 

(5.21c) 

=  0  > 

(5.21(i) 

^l{xN,y,i)  =  o  , 

(5.21e) 

=  i^2ss{^N,y) 

(5.21/) 

and  for  fc  >  0,  become 

U'‘{x/^,y,i)  =  0  etc.  (5.22) 


Freestream  Boundaries  (y  =  yo)  and  (y  =  y^f) 

In  Chapter  4,  it  was  shown  that  a  disturbance  described  by  linear  stability 
theory  behaves  like 

'{x,y,z,t) 


■{x,y,z,t)  =  ±r2A  al  +  -ju 


(5.23) 
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as  y  — f  ±oo.  For  two-dimensional  disturbances,  equation  (5.23)  simplifies  to 

du* 

=  ±r2aru'(x,y,z,i)  .  (5.24) 

However,  as  noted  in  Chapter  4,  equations  (5.23)  and  (5.24)  are  not  easily  em¬ 
ployed  in  their  present  form.  The  freestream  boundary  conditions  are  dealt  with 
more  easily  in  terms  of  the  Fourier  transformed  variables. 

If  equation  (4.28)  is  Fourier  decomposed  with  respect  to  the  spanwise  direction 
z,  then  satisfies  boundary  conditions  similar  to  equation  (5.23)  and  17®  satisfies 
boundary  conditions  similar  to  equation  (5.24).  One  other  question  concerns  the 
value  of  Or  in  these  equations.  For  i  =  0,  Or  is  set  equal  to  which  is  associated 
with  the  two-dimensional  excitation  at  the  inflow  boundary  (see  equations  4.17). 
For  fc  =  1,  Or  is  set  equal  to  which  is  associated  with  the  three-dimensional 
excitation  at  the  inflow  boundary  (see  equations  4.18).  For  fc  >  1,  the  wake  is 
assumed  to  be  undisturbed  at  the  freestream  boundaries. 

The  Fourier  representation  of  the  lower  freestream  boundary  conditions  (4.28), 
(4.30a),  (4.31a),  and  (4.32)  are  shown  below.  For  fc  =  0,  the  boimdary  conditions 
are 


ar/'o 

Qy  (®,yo,t) -»-2aJ'^C7'®(®,yo,<)  , 

(5.25a) 

av® ,  ,  du\ 

(5.256) 

W°  =  0  , 

(5.25c) 

fi*(®,yo,t)  =  o  , 

(5.25d) 

(5.25c) 

and  n®(a;,yo,t)  =  0  . 

(5.25/) 
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As  mentioned  previously,  primes  denote  the  perturbation  component  of  a  variable. 
For  jb  =  1,  the  boundary  conditions  are 


dV^  dU^ 

dW^ 
dy 

fii(a:,yo,0  =  0  » 


(a:,yo,0  =  *'2yaJ^N-^w^^(®»yo,0  , 


^^i(a:>yo,<)  =  0  » 

and  Ql{x,yo,t)  =  Q  . 


For  k  >  1,  the  boundary  conditions  are 


(5.26a) 

(5.266) 

(5.26c) 

(5.26d) 

(5.26e) 

(5.26/) 


U'^(x,y(i,i)  =  0  etc. 

The  upper  freestream  boundary  conditions  are  similar  to  those 
freestream  boundary.  For  k  =  0,  these  conditions  are 

QU'^ 

-^(®,yM,0  =  -»-2a“t/'°(z,yAf,<)  , 

dV\  .  dU\ 

-^{x,yMyt)  =  — ’ 

W°  =  0  , 

n*(®,yAf,0  =  o  . 

n5(®,yM,0  =  o  > 

and  f2°(x,yAf,0  =  0  * 

For  fc  =  1,  the  boundary  conditions  are 


(5.27) 
at  the  lower 


(5.28a) 

(5.286) 

(5.28c) 

(5.28d) 

(5.28c) 

(5.28/) 


dU^ 

-^(*,yM,0 


U^ix,yM,t) 


(5.29a) 
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-^(a5,yM,<)  =  — -iW^{x,yM,t)  ,  (5.296) 


=  -r2^al^^  +  ^W^{x,yM,i)  ,  (5.29c) 

n*(*,yM,0  =  0  ,  (5.29«i) 

nl{x,yM,i)  =  0  ,  (5.29e) 

and  Q,l{x,yM,t)  =  0  .  (5.29/) 

For  A:  >  1,  the  upper  freestream  boundary  conditions  are 

U'‘{x,yM,t)  =  0  etc.  (5.30) 


Equations  (5.25a)  and  (5.28a)  enforce  transverse  exponential  decay  of  the  two- 
dimensional  velocity  disturbaxH:es.  Equations  (5.26a),  (5.26c),  (5.29a),  and  (5.29c) 
enforce  transverse  exponential  decay  of  the  three-dimensional  velocity  disturbances 
that  correspond  to  the  first  spanwise  mode  fc  =  1. 

Initial  Conditions  (t  =  0) 

As  initial  conditions  for  the  calciilation  of  the  disturbed  wake,  the  velocity 
and  vorticity  distributions  of  the  two-dimensional  bzise  flow  are  prescribed  in  the 
entire  integration  domain.  Thiis,  for  ib  =  0  the  initial  conditions  are 


I7‘’(z,y,0)  =  1*55(1,  y)  , 

(5.31a) 

^”(®»y,0)  =  V55(®,y)  > 

(5.315) 

W\x,y,Q)  =  0  , 

(5.31c) 

n®(*>y.o)  =  o  , 

(5.31d) 

fiJ(*,y.o)  =  o  , 

(5.31e) 

J^°(*»y»o)  =  w*ss(*>y) 

(5.31/) 
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zind  for  k  >  0,  the  initial  conditions  zire 

U^{x,y,0)  =  0  etc.  (5.32) 

5.2  Computational  Domain 

As  a  result  of  the  Fourier  decomposition  of  the  dependent  variables,  the  spatial 
domain  reduces  to  a  two-dimensional  plane  in  space  (z,y)  that  is  perpendicular 
to  the  z-ajds.  This  domain  is  discretized  into  a  finite  number  of  uniformly  spaced 
grid  points.  These  grid  points  have  coordinates  (in,ym)  given  by  the  relations 

=  lo  +  nAx  for  n  =  0, . . . ,  JV  (5.33o) 

and  ym  =  yo  +  mAy  for  m  =  0,  ...,M  .  (5.336) 

The  number  of  grid  increments  in  the  streamwise  and  transverse  directions  are 
denoted  by  N  amd  M  respectively.  The  constants  Ax  and  Ay  are  the  grid  sizes 
in  the  streamwise  and  transverse  directions  respectively.  The  discretized  spatial 
domain  is  shown  in  Figure  5.1. 

The  temporal  domain  is  divided  into  discrete  uniformly  spaced  time  levels  ti, 
such  that 

ti  =  lAt  for  1  =  Xi,...,0,  ...,^2  (5.34) 

where  At  is  the  time  step  between  each  time  level.  The  base  flow  calculation  is 
performed  in  the  time  interval 


tL,  <ii<0 


(5.35) 


which  corresponds  to 


Li<l<0  . 


(5.36) 
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The  disturbed  flow  is  calculated  in  the  time  interval 


0  <  <  ti, 


(5.37) 


which  corresponds  to 

0  <  1  <  I2  .  (5.38) 

Approximate  numerical  solutions  of  the  governing  equations  (5.4)  are  obtained 
at  the  spatial  and  temporal  locations  (xntymt^i)*  approximate  solutions  and 
their  Correspondence  to  the  spatial  and  temporal  grid  points  are  denoted  by 


~  ^  (®n»yin»^l)  ctc. 


(5.39) 


5.3  Discretization  of  Spatial  Derivatives 

All  spatial  derivatives  in  equations  (5.4)  are  approximated  using  the  second- 
order  finite-difference  formulas  shown  below.  Consider  a  complex  function  4>{(,ri) 
with  known  values  at  the  equally  spaced  points,  (fo>^o)>  (^i>^o)>  •  •  •  >  (^/-i)^j)> 
{(itlj)-  Then  the  first  partial  derivative  ^^(|^^can  be  approximated  by 


and  the  second  partial  derivative  ^  can  be  approximated  by 

82^  -  ^^2  ) 


(5.40) 


(5.41) 


where  is  the  grid  increment  in  the  (  direction. 

5.4  Discretized  Vorticitv  Equations 

A  combination  of  different  methods  are  used  to  solve  the  vortidty  equations 
(5.4a)  through  (5.4c).  The  methods  have  been  chosen  because  of  their  numerical 
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stability,  their  accuracy,  and  their  computational  efficiency.  The  Altemating- 
Direction-Implicit  (ADI)  method  is  used  to  discretize  the  streamwise  and  trans¬ 
verse  diffusion  terms.  This  method  wais  selected  because  it  is  second-order  time- 
accurate  and  is  unconditionally  stable  when  used  for  solving  linear  diffusion  equa¬ 
tions.  Use  of  the  Crank-Nicolson  method  for  the  spanwisc  diffusion  term  re¬ 
tains  both  second-order  time-accuracy  and  good  stability  characteristics.  Fi¬ 
nally,  the  second-order  Adams- Bashforth  method  is  used  for  the  nonlinear  terms. 
This  method  possesses  favorable  stability  characteristics  and  is  second-order  time- 
accurate.  The  spatial  derivatives  are  discretized  using  equations  (5.40)  and  (5.41). 
The  resulting  discretized  vorticity  equations  can  be  solved  using  noniterative  meth¬ 
ods. 

This  choice  of  numerical  techniques  results  in  a  two-step  method  for  inte¬ 
grating  the  vorticity  equations.  In  the  first  integration  step,  the  streamwise  dif¬ 
fusion  term  is  explicit  and  the  transverse  diffusion  term  is  implicit.  In  this  step, 
the  vorticity,  which  is  known  at  tj  =  /At,  is  obtained  at  the  intermediate  time 
tj^i  =  (/  +  j)At.  For  the  first  integration  step,  the  discretized  streamwise  vortic¬ 
ity  equation  is 


ADI— Explicit 


At/2  ~  Re  \  Ax^ 


AOI— Implicit 


Crank  — NicoUon 


+ 


D*  , ,  , ,  1  —  2D*  , .  1  +  D*  ,  , ,  1  l2  D*  , .  I  -I-  D 


r|Ay^ 

Adam*  — Baahforth 

/<  . .  s 

_  f  s,n,m,{  — 

“  V  O  /  ’ 


2 


(5.42a) 


63 


the  discretized  transverse  vortidty  equation  is 

AOI— Explicit 


^V,n,Tn,l+^  _ ^  ^v,n-l,m,l 


A</2 


Re\ 


Ax" 


ADI— Implicit 


Crank  — Ni  colson 


+ 


^Ia?  rl^  2 


Adams  — Baskfortk 


3J*’*  _  J!’* 

_  ^  V.n.m.l 

-  V  2  ^  ’ 

and  the  discretized  spanwise  vorticity  equation  is 


(5.426) 


n'‘ 


ADI— Explicit 


_ ±f: _ 

At/2  Re  V  Ai^ 


ADI— Implicit 


Crank— Nicolson 


+ 


rjAy^  rV  2  ’) 


Adams  — Bashforth 
=  (  2  ' 


(5.42c) 

In  the  second  integration  step,  the  streamwise  diffusion  term  is  implicit  and 
the  transverse  diffusion  term  is  explicit.  In  this  step,  the  vorticity  is  obtained  at 
t  =  (/  4-  l)At.  For  this  step,  the  discretized  streamwise  vorticity  equation  is 

ADI— Implicit 


At/2  Re  \  Ax^ 


ADI— Explicit 


Crank— Nicolson 


+ 


rjAy^ 

Adams— Bashforth 
_  I?* 

—  (  2  '  ’ 


■)) 


(5.43o) 
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the  discretized  transverse  vorticity  equation  is 


ADI— Implicit 


At/2 


-U 


A®' 


ADI— Explicit 


Craak— NicoUon 


+ 


.  1  -  2fl* 


^y,n,m— 1,1+^  ^y,n,m,l+'j  \\ 

W  2  V 


Adam* — Baahfortk 


—  F* 

—  V  /  » 


(5.436) 


and  the  discretized  spanwise  vorticity  equation  is 


ADI— Implicit 


^z,n,m,t4-l  ^z,n,m,l+-^  1  /  1  ^^x,n,in,l+l  ^«,n— 

-  Re' 

ADI- Explicit 


_ 

At/2  Re  V  A®2 


Craak  —Nicolson 


+ - - -j{ - 5 - 


r|Ar 

Adams— Bashforth 


SF*' 

^  z,n,fTi,44* 'j  *»^i^**^ 


(5.43c) 


Equations  (5.42)  and  (5.43)  are  valid  in  the  interior  of  the  spatial  domain  at 
the  points  (in>ym)  where 


1  <  n  <  ^■  -  1  (5.44a) 

and  1  <  m  <  Af  —  1  .  (5.446) 

Equations  (5.42)  and  (5.43)  are  also  valid  for  times  U  where 

t  =  Zrl,  .  .  .  ,  0,  .  .  .  ,  X2  .  (5.45) 
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At  the  domain  boundary,  the  boundary  conditions  from  Sections  5.1.1  and  5.1.2 
are  used  to  obtain  the  vorticity  components. 

5.4.1  Discretized  Boundary  and  Initial  Conditions  for  the  Undisturbed  Vorticity 
For  the  base  flow  calculation,  the  vorticity  boundary  and  initial  conditions 
from  Section  5.1.1  are  discretized.  These  boundary  conditions  axe  discussed  here. 
At  the  inflow  boundary,  the  spanwise  vorticity  is 

=  ‘*>iiym)  •  (5.46) 


At  the  outflow  boundary,  equation  (5.13c)  combined  with  equation  (5.42c) 
results  in 

_ 


At/2  Re 

opO  _  pO 

_  _  )  . 


Equation  (5.13c)  combined  with  equation  (5.43c)  results  in 


(5.47) 


n! 


Aij2  V  rlAy^  J 


o  jpQ  _  pO 

=  (  2  ^  ' 


(5.48) 


At  the  freesteam  boundaries,  the  spanwise  vorticity  is 


Kn,0,l  =  0 


and  =  0  • 


The  initial  condition  for  the  base  flow  vorticity  is 


(5.49a) 

(5.496) 


n®,n,m,L,  = 


(5.50) 


t 


66 


5.4.2  Discretized  Boundary  and  Initial  Conditions  for  the  Disturbed  Vorticity 
The  boundary  and  initial  conditions  from  Section  5.1.2  are  used  for  the  cal¬ 
culation  of  the  disturbed  vorticity.  The  discretized  versions  of  these  boundary 
conditions  sure  discussed  here. 

At  the  inflow  boundary,  the  vorticity  boundary  conditions  for  which  A:  =  0 
are 


and 


®  » 

^2,0,tn,i  —  ^xSsi^Otym)  +  -Pw,  (®0>  J/m»  tj) 


For  fc  =  1,  the  boundary  conditions  are 

and  o  j  =  • 

For  ib  >  1,  the  vorticity  components  at  the  inflow  boundary  are 


n 


k 


=  0 


etc. 


(5.51a) 

(5.51i) 

(5.51c) 


(5.52a) 

(5.526) 

(5.52c) 


(5.53) 


At  the  outflow  boundary,  the  vorticity  boundary  conditions  for  which  k  =  0 
are 


(5.54a) 

=  0  » 

(5.546) 

(5.54c) 
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and  for  fc  >  0,  are 

^z,N,m,i  =  0  etc-  (5-55) 

At  the  freestream  boundaries,  the  wake  is  assumed  to  be  irrotational.  For 
fc>0,  the  lower  freestream  boundary  conditions  are 


(5.56) 

For  k  >  0,  the  upper  freestream  boundary  conditions  are 

, 

=  0  etc. 

(5.57) 

Finally,  the  discretized  initial  conditions  for  the  disturbed  vorticity  calciilation 

are  presented.  For  fc  =  0,  the  initiail  conditions  are 

(5.58a) 

n®  —  n 

(5.586) 

^z,n,m,0  —  ^zSsi^ntVm) 

(5.58c) 

and  for  k  >  0,  are 

=  0  etc.  (5.59) 

5.4.3  Equation  Systems  for  the  CiJculation  of  the  Vorticity 

The  discretized  vorticity  equations,  (5.42)  and  (5.43),  combined  with  the  dis¬ 
cretized  vorticity  boundary  conditions  form  systems  of  linear  algebraic  equations 
that  are  solved  for  the  unknown  vorticity.  To  obtain  the  vorticity  at  the  interme- 
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diate  time  tj+j,  the  tridiagonal  system 


(hi  ai 
di  h\ 


f  \ 


\ 

/ 


di 


hi 

di 


hj 


\  liJ+'j  } 


(5.60) 


\ 


IS 


+  ^ni^z^n+i^M-1,1  +  ^z,n-l,M-l,l)i^  ~  ^N.n)  +  ■ffn.M-l,/'' 
solved.  For  the  base  flow  cailculation,  equation  (5.60)  is  solved  for 


jfc  =  0  ,  1  <n  <  N  ,  and  Li  <  I  <0  . 

For  the  disturbed  flow  calculation,  equation  (5.60)  is  solved  for 

jr 

0  <  k  <  —  ,  l<n<JV  —  1  ,  and  0  <l  <  L2 

2 


The  constzmts  in  equation  (5.60)  su’C 


At 


2i2e(r2Ay)^ 


hi  =  1-1- 


Ai 


+ 


At  /^Y 
4i2e  \r3  J 


Cl  = 


ReirzAy)^ 

i!n=N; 


and  d 


I  0  if  n  =  iV. 


(5.61) 

(5.62) 

(5.63d) 

(5.635) 

(5.63c) 

(5.63d) 
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At. 


The  function  H!!  _  ,  is 

^ 

and  ^N.nj  the  Kronecker  delta,  is 

‘"■”  =  {0;  n 


(5.64) 


=  JV; 

j^N. 


(5.65) 


To  obtain  the  unknown  vorticity  at  time  ti+iy  the  tridiagonal  system 


^  62  ^2 

02  62  0,2 


V 


02  62  02 

o;  hi) 


Vn* 


(5.66) 


^  1  1  +  d2(fi^  1  02nj,,,„^j^j  ^ 

^2^2,2,m,J+i  2,m+l,I+|  ^z,2,m-l ,/+ j  )  ■*" ‘^2,7n,J+ 1 


‘^2(f^*^;V',m+l,i+A+  ^z,JV',m-l.J+p  ■*“  ^N',m,l+\~  “2^*,A?,m,i  +  l  ^ 
is  solved.  For  the  base  flow  csJculation,  equation  (5.66)  is  solved  for 


N'  =  N  ,  Jb  =  0  ,  l<m<M-l  ,  and  Li  <  I  <  0  .  (5.67) 

For  the  disturbed  flow  calculation,  equation  (5.66)  is  solved  for 

IT 

N'  =  N-l  ,  0  <  ifc  <  —  ,  l<Tn<M-l  , 

2 


and  0  <  I  <  L2  •  (5.68) 
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The  function  l  is 

The  constants  in  equation  (5.66)  are 

At 

“^““2i2cA®2  ‘ 

^  ReAx^  ^  4Re  \ri  ) 

At  At 

^  4Re  Vrs/  Re(r2Ayy 

di  2Re{r2AyY 

The  constcint  a!;  is 

,  _  J  0,  for  the  base  flow  calculation; 
^  ~  \  ®2  5  otherwise. 


(5.69) 


(5.70a) 

(5.706) 

(5.70c) 

(5.70<i) 


(5.71) 


The  constant  63  is 


f  1  +  for  ffl®  base  flow  calculation; 

i  62,  otherwise. 


(5.72) 


Equations  (5.60)  and  (5.66)  are  solved  for  the  spanwise  vorticity.  Equation 
systems  for  the  streamwise  and  transverse  vorticity  components  are  identical  to 
equations  (5.60)  and  (5.66)  and  are  not  shown. 

For  the  base  flow  calculation,  equations  (5.60)  and  (5.66)  are  solved  beginning 
with  the  initial  data  at  time  t  =  t^^.  For  the  disturbed  flow  calculation,  equations 
(5.60)  and  (5.66)  are  solved  beginning  with  the  initial  data  at  time  t  —  0.  For  the 
initial  integration  step  of  both  the  base  flow  and  the  disturbed  flow  calculations, 
the  nonlinear  terms  are  evaluated  using  the  first-order  Euler  method  instead  of 
the  second-order  Adams-Bashforth  method. 
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5.5  Discretized  Velocity  Equations 

The  Helmholtz  equations,  (5.4d)  through  (5.4f),  are  discretized  using  the 
second-order  finite-difference  formulas  derived  in  Section  5.3.  The  discretized 
Helmholtz  equations  are,  for  the  streamwise  velocity  component 
,  ,  -  2f/*  ,  +  I/*"  ,  , 


+ 


Ai’ 


(^)  =  ,  (5.73a) 


for  the  transverse  velocity  component 
yfc  _  ‘2V’*'  -t-  V* 

'n— l,m,/ 


+ 


Ai^ 

V*  —  4.  V'*' 

rlAy^ 


f~ik 


(5.736) 


and  for  the  spanwise  velocity  component 


+ 


Ax' 

K 


w. 


n,m,l 


=  G 


z,n,m,l 


(5.73c) 


The  right  hand  sides  of  equations  (5.73)  are 


2Ay  J  rj  ’ 

^3  f  ^x,n+l,m,l  ^x,n— 


ik 

rsrl 


2Ax 


(5.74a) 

(5.746) 
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and  =  — 


—  n* 

V.n+l,m,l  **y.n-l,m,< 


2A2 


1  f  ^a,n,m+l,i  1,/ 

2Ay 


rjr| 


(5.74c) 


The  Helmholtz  equations  are  solved  for  each  integration  step,  subject  to  the  ap¬ 
propriate  boundary  conditions  from  Sections  5.1.1  and  5.1.2. 

Equations  (5.73)  are  valid  in  the  interior  of  the  spatial  domain.  At  the  domain 
boundary,  equations  (5.73)  combined  with  the  discretized  boundary  conditions  are 
used  to  obtain  the  velocity  components. 

5.5.1  Discretized  Boundary  and  Initial  Conditions  for  the  Undisturbed  Velocity 
For  the  base  flow  calculation,  the  velocity  boundary  and  initi2j  conditions 
from  Section  5.1.1  are  discretized  as  discussed  below. 

At  the  inflow  boundary,  the  streamwise  and  transverse  velocity  components 
are 


and  =  0  . 


(5.75a) 

(5.756) 


The  discretized  version  of  the  outflow  boundary  conditions,  (5.13a)  and 
(5.13b),  combined  with  equations  (5.73a)  and  (5.73b)  resiilt  in 


(r2Ay)2 

for  the  streamwise  velocity  and 

(rjAy)* 


—  r*® 


(5.76) 


^0 


(5.77) 


for  the  transverse  velocity. 
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At  the  freestream  boundaries,  the  streamwise  velocity  is 


(5.78a) 

(5.786) 


Equations  (5.14b)  and  (5.15b)  combined  with  equation  (5.73b)  result  in 


.  2C...  -  2V„%.I 


at  the  lower  freestream  boundary  and 


=  0 


I  2K.V-..1  - 


Ai' 

at  the  upper  freestream  boundary. 
The  auxiliary  condition 


r|As= 


V“  „  ,  =  0 


(5.79) 


(5.80) 


(5.81) 


imposes  antisymmetry  on  the  transverse  velocity  field.  Because  of  this  antisym¬ 
metry,  the  transverse  velocity  j  is  solved  only  in  the  lower  half  of  the  spatial 
domain  for  which 


and 


1  <  n  <  N 

0  <  m  <  ^  -  1  . 

-  -  2 


(5.82a) 

(5.826) 


The  transverse  velocity  in  the  upper  half  of  the  spatial  domain  is  found  from  the 
relation 


—  _V® 

+  n,'^— m,J 


for 


-  2 


and  1  <  n  <  iV 


(5.83) 


Finally,  the  initial  conditions  for  the  streamwise  and  transverse  velocities  are 

{5.84o) 

and  =0  •  (5.841.) 
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5.5.2  Discretized  Boundary  and  Initial  Conditions  for  the  Disturbed  Velocity 
Boundary  eind  initial  conditions  from  Section  5.1.2  are  used  for  the  calcxilation 
of  the  disturbed  velocity.  The  discretized  versions  of  these  boundary  and  initial 
conditions  are  discussed  in  this  section. 

At  the  inflow  boundary,  the  velocity  components  for  Jb  =  0  are 


“  ^Ss(®0»yTn)  +  Pu  > 

(5.85a) 

(5.856) 

and  =  0  • 

(5.85c) 

For  =  1,  the  velocity  components  are 

(5.86a) 

(5.866) 

and  . 

(5.86c) 

For  k  >  1,  the  velocity  components  are 

^o,m,i  =  0  etc. 

(5.87) 

At  the  outflow  boundary,  the  velocity  components  for  fc  =  0  are 

=  y’Ss{xs,ym)  , 

(5.88a) 

=  »5s(*N,ym)  , 

(5.886) 

n,nv.l  =  0 

(5.88c) 

and  for  jfc  >  0,  are 

=  0  etc- 


(5.89) 
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At  the  lower  freestream  boundary,  equation  (5.25a)  combined  with  equation 
(5.73a)  results  in  the  following  equation  for  the  fc  =  0  streamwise  velocity: 

Ax’  T\^y^ 

=  (2uss(gn>yi)  -  (2  +  2C'’‘^)u5s(*n,yo))  (5.90) 

where  C’**  =  Equation  (5.25b)  combined  with  equation  (5.73b)  results 


in 


Ax’ 

for  the  fc  =  0  transverse  velocity.  For  fc  =  0,  the  spanwise  velocity  is 


<0./  =  0 


(5.91) 


(5.92) 


at  the  lower  boundary. 

For  fc  =  1,  equation  (5.26a)  combined  with  (5.73a)  results  in  the  following 
equation  for  the  streamwise  velocity: 


+ 


Ax’ 

-  (2  + 

r|Ay’ 


(5.93) 


where  C’**  =  Ayr2 ya*^^-|-(^^.  Similarly,  equation  (5.26c)  combined  with 
(5.73c)  results  in 


Ax’ 

^K,u-i^  +  2Ci^WUi 

rlAy’ 


+ 


(5.94) 
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for  the  jb  =  1  spanwise  velocity.  Equation  (5.26b)  combined  with  equation  (5.73b) 
results  in 

r|Av’  V's/ 

which  is  valid  for  the  k  =  1  transverse  velocity.  For  >  1,  the  velocity  components 
at  the  lower  freestream  boundary  are 


Un,Q,l  =  ^  etc- 


(5.96) 


At  the  upper  freestream  boundary,  the  following  equations  are  valid.  For 
k  =  0,  the  velocity  boundary  conditions  are 


rlAy^ 


and 


=  (2u55(gn>yAf-i)  -{2  +  2Cl^)uss{xn,yM))  »  (5.97a) 


2v;m,, + v.°-. 

Az^  rlAy^ 


Km.i  =  0  •  (5.97c) 


For  A:  =  1,  the  velocity  boundary  conditions  are 

Az^ 

rlAy^ 


98a) 


77 


A*’ 


,(5.3.) 


and 


-  ^Ko.l  +  -1M,1  ,  ,M,l  f  1  ^  ^  T/1 

- zp - — ;iz? - UJ 


=  (;;|;jT(|f(t^i+i.M,.+  l^;-..M..)  +  iAyW'i„,,)  .  (5.98c) 

For  fc  >  1,  the  velocity  components  at  the  upper  freestream  boundary  Me 


(5.99) 


Finally,  initial  conditions  must  be  specified  for  the  velocity  components.  For 
A:  =  0,  the  initial  conditions  are 


^n,m,0  ~  ^5s(®n>  Vm)  > 

(5.IOO0) 

^n,m,0  ~  ®Ss(®n»yTn)  > 

(5.1006) 

O 

II 

o 

i 

(5.100c) 

and  for  A;  >  0,  the  initial  conditions  are 

K,m,o  =  0  etc. 

(5.101) 

5.5.3  Eauation  Svstems  for  the  Calculation  of  the  Velocity 

The  discretized  velocity  equations  (5.73)  combined  with  the  discretized  veloc- 

ity  boundary  conditions  form  systems  of  linear  zdgebraic  equations  that  are  solved 
to  obtain  the  velocity  components. 
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Base  flow  Streamwise  Velocity 

Equation  (5.73a)  combined  with  boundary  conditions  (5.75a),  (5.76),  and 
(5.78)  form  a  system  of  linear  algebraic  equations  that  is  solved  for  the  base  flow 
streamwise  velocity.  To  derive  this  system  of  equations,  let 
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Then,  the  resulting  system  of  equations 


(A'  I 
I  A'  I 


\ 

/  Az^gi  -  hiA  -^o\ 
I  Ai^gj  -  h2A  I 


I  A> 
0 


NxN 


\^N/ 


Nxl 


(5.103) 


\  Ai^giv  -  hivA  /  ;vxi 


is  solved  for  the  base  flow  streamwise  velocity.  The  constant  A  is 


the  matrix  A'  is 


>l'  =  Ayi-2/  , 


(5.104) 


(5.105) 


and  I  is  the  identity  matrix. 

Base  flow  Transverse  Velocity 

Equation  (5.73b)  combined  with  boundary  conditions  (5.75b),  (5.77),  (5.79), 
and  (5.80)  form  a  system  of  linear  algebraic  equations  that  is  solved  for  the  base 
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flow  transverse  velocity.  To  derive  this  system  of  equations,  let 


/-2  1 

1  -2  1 


and  A  — 


Then,  the  system  of  equations 


/A'  / 

I  A'  I 


Sn  = 


,  (5.IO60) 


'  ^  if  Xl 


\ 


•  • 


1  -2  1 
2  -2/ 


(5.1066) 


ifx^ 


\ 


/ 


V 


\  /  Ai^gi  \ 


(5.107) 


\Ax^gN/ Nxl 


I  A’  I 

is  solved  for  the  base  flow  transverse  velocity  in  the  lower  half  of  the  spatial  domain, 
for  which 


I <n<N  , 

M 

and  0  <  m  <  — —  1 
-  2 


(5.108a) 

(5.1086) 


81 


The  base  flow  transverse  velocity  in  the  upper  half  of  the  spatial  domain  is  obtained 
from  the  relation 


n,-y+TH,t  n 


M 

1  <  m  <  — 


(5.109) 


As  previously  , 


jind 


A'  =  KA- 21  . 


(5.110a) 

(5.1106) 


Disturbed  Velocity  Components  for  k<\ 

Equation  (5.73a)  combined  with  the  boundary  conditions  (5.85a),  (5.88a), 
(5.90),  and  (5.97a)  form  a  system  of  equations  that  is  solved  for  the  zeroth  spanwise 
mode  of  the  disturbed  streamwise  velocity  To  derive  this  system  of  equations. 


f>2 


(  Km,  \ 

0 


and  hm  = 


(5.111c) 


(N-l) 

Then,  the  system  of  equations 


f  A'  -  B'  2M 
AI  A'  A1 


\ 


/*o\ 


V 


!  Ai^go  -  ho  +  2AyAf' 
Ai^gi  -  hi 


AI  A*  AI 
2AI  >1' + 


<M+l)xl 


(5.112) 


Ai^gjvf-i  -  hAf-i 
\  Az^gjvf  —  hAf  —  2AyAf"  /  (Af+i)xi 

is  solved  for  the  zeroth  spanwise  mode  (fc  =  0)  of  the  disturbed  streamwise  velocity. 
In  equation  (5.112),  A'  is 

A'  =  (^)  (5.113a) 

and  A  is 


A=r^v 


VrjAy/ 


(5.1136) 
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The  matrix  A'  is 

A' =  A  -  (2A  +  A')/  ,  (5.114) 

the  matrix  B'  is 

B'  =  2AAyc'J  ,  (5.115) 

and  the  matrix  B"  is 

B"  =  2AAyc'I  .  (5.116) 

The  vectors  f'  =  [/^]  and  f"  =  [/"]  and  the  sccdars  c'  and  c"  result  from 
the  freestream  boundary  conditions.  For  the  zeroth  spanwise  mode  (ifc  =  0)  of  the 
streamwise  velocity  C7°, 

ft  2d  /  \  .  ^Ssl^^niVl)  '^Ss(^ni  J/o)  /  r  1  1  \ 

fn  =  -^r  r2Uss(Xn,yo)+  - ^  ,  (5.117a) 

Ay 

and 

c'  =  a^‘^r2  ,  (5.117c) 

c"  =  -a$S  •  (5.117d) 


Equation  systems  for  the  velocity  components  V®,  17^,  ,  and  are  almost 

identical  to  equation  (5.112).  The  only  difference  between  equatic.  i  (5.112)  and 
systems  of  equations  for  the  other  velocity  components  V°,  Z7\  ,  and  are 

the  vectors  f'  =  [/^]  and  f"  =  [/"]  and  the  constants  c'  and  c".  Therefore,  for 
these  other  velocity  components,  only  /^,  /^,  c'  and  c"  will  be  shown. 

For  the  zeroth  spanwise  mode  {k  =  0)  of  the  transverse  velocity  F®, 


f/®  —  TJ^ 

ft  _  *-^n-n,o,i 

“  2Ax 

u°  —  r/® 

fii  _  ^n+l,M,l 

2Az 


(5.118a) 


(5.118i) 
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and 


c'  =  0  ,  c"  =  0 


(5.118c) 


For  the  first  spanwise  mode  (fc  =  1)  of  both  the  streamwise  and  spanwise 
velocities,  and 


/:  =  o  , 

fn=0 


(5.119a) 

(5.1196) 


and 


c"=r,y^a>'''+  (i)  . 


(5.119c) 

(5.119(i) 


Finedly,  for  the  first  spanwise  mode  {k  =  1)  of  the  transverse  velocity  , 


f  _  ^n+l,0,l  ^n-l,0,l  .^1 

2Ax  ’ 

ell  _  ^n+l,Af,/  -xTrl 


(5.120a) 

(5.1206) 


and 


c'  =  0  ,  c"  =  0 


(5.120c) 


Disturbed  Velocity  Components  for  ib  >  1 

Equation  (5.73a)  combined  with  the  boundary  conditions  (5.87),  (5.89),  (5.96), 
and  (5.99)  form  a  system  of  equations  that  is  solved  for  the  disturbed  streamwise 


J 
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velocity  for  fc  >  1.  To  derive  this  system  of  equations,  let 


-2 

1 


(5.1216) 


(N-l)x(/V-l) 
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Then,  the  system  of  equations 

(A!  AI 
A/  A'  AI 


V 

/  \ 


AI  A'  AI 

^  (M-l)x(M-l) 


i  ^ 


\^M  -1  /  (M-l)xl 


(5.122) 


\Ax^gAf-l/  (Af_i)xl 

is  solved  for  the  streamwise  velocity  component  U*.  As  before, 

A'=(^' 

\  ’'3  . 


and  A  =  f 


The  matrix  A'  is 


A'  =  A  -  (2A  +  A')/ 


(5.123a) 

(5.1236) 

(5.124) 


Equation  systems  for  W'’®  and  for  F*  and  TF*  with  k  >  \,  are  identical  to  equation 
(5.122)  and  will  not  be  shown. 

Equations  (5.103),  (5.107),  (5.112),  and  (5.122)  are  systems  of  linear  alge¬ 
braic  equations  that  are  solved  for  the  velocity  components.  For  each  system, 
the  coefficient  matrix  is  large  and  sparse.  However,  because  of  the  structure  of 
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these  matrices,  feist  Helmholtz  solvers  can  be  used  to  solve  these  equation  systems. 
Details  of  these  techniques  are  given  by  Swarztrauber  (1977). 

5.6  Evaluation  of  the  Nonlinear  Terms 

Up  to  this  point  in  the  development  of  the  numerical  method,  the  evaluation 
of  the  nonlinear  terms  and  _  i  at  each  time  step  has  been 

ignored.  However,  in  order  to  calculate  and  from  equa¬ 

tions  (5.42)  and  (5.43),  the  nonlinear  terms  must  be  computed  at  times  ,  ti, 
and  Generally,  computing  the  nonlinear  terms  is  the  most  time  consuming 

aspect  of  using  spectral  methods  to  solve  the  Navier-Stokes  equations.  The  diffi¬ 
culties  arise  from  the  fact  that  the  nonlinear  terms  do  not  have  a  simple  Fourier 
representation. 

The  nonline2ir  terms  can  be  computed  either  pseudo-spectrally  or  spectrally. 
The  pseudo-spectral  method  works  in  the  following  manner.  The  physical  repre¬ 
sentation  of  the  velocity  and  vorticity  are  obtained  from  the  Fourier  coefficients 
using  equations  (5.1).  Then,  the  nonlinear  terms  and  are 

computed  from  the  formulae 


k  k  k  k 

rk  _  fc  ^n— l,Tn,i  ,  k  ^n,Tn-H,i 

J  z,n,m,l  ~  ^x,n,m,l  2Ax 


2Ay 


du 


4.^,*  — I*  u* 

r  ln,Tn,4  n.m, 


07^  —  (J** 


2Ax 


—  V 


k  ‘^z,n,Tn+l,l  fc  ifc  /cioc-X 

ln,m,<  »  (5.125a) 


n,m,l 


2Ay 


fk  _  k 


2Ax 


,  fc  ^jk  _  u* 


fc  .  .fc 


2Ay 


dz 


2Ax 


k  k 

„fc  ~ fc  ^y  ifc  (k, 

®n,m,J  2Ay  g^  ln,m,l  >  (5.1256) 
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f z,n,m,<  ^x,n,m, 


^n,m+l,Z  ^n,m  — 1,Z 


2Ax 


y,n,m,r 


2Ay 


k  k 

,  k  ^^\k  _  *,fc  ‘*^z,n+l,m,Z  ~  ^z,n-l,m,l 


2Az 


ib  ik 

^n,m,/  2Ay  “’n.m.I  In, m,Z  5  (O.lZOCj 


which  are  the  finite-diflerence  versions  of  equations  (3.12). 

The  z  derivatives  in  equations  (5.125)  are  computed  using  the  formula 


|^{x)  =  £  ikS^e 

*=-f 


ikz 


(5.126) 


where 


Jf 

T 


/i(z)=  ^  . 


(5.127) 


L-.  K 
k - -j- 


The  nonlinear  terms  fz,n,m,i  "e  transformed  back  to 

Fourier  space  using  equations  (5.5).  This  results  in  the  nonlinear  terms 
Fy,n,m,o  ^  ^  pseudo-spectral  method,  together 

with  the  use  of  fast  Fourier  transforms,  has  an  aisymptotic  operation  count  of 


0{NMKLoz2K)  . 


(5.128) 


An  alternative  to  the  pseudo-spectral  method  is  the  spectral  method.  For  the 
representative  nonlinear  term 

^  (5.129) 


*dx  ’ 


the  spectral  method  proceeds  as  follows:  The  Fourier  series  representation  of 
the  dependent  variables  u(zn,ym?2fc|t()  and  substituted  into 
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equation  (5.129)  to  get 

*  dx ~ 

^  f;  ^  Y,  •  (5-130) 

The  multiplication  of  the  two  series  in  equation  (5.130)  is  performed  in  order 
to  obtain  the  Fourier  series  representation  of  equation  (5.129).  This  procedure 
is  repeated  for  the  other  terms  that  make  up  the  nonlinear  part  of  the  govern¬ 
ing  equations.  In  this  way,  the  Fourier  representation  of  the  nonlinear  terms  is 
obtained.  The  asymptotic  operation  count  for  the  spectred  method  is 

0{NMK^)  .  (5.131) 

For  large  values  of  K,  the  pseudo-spectral  method  reqmres  fewer  operations 
than  does  the  spectrad  method.  However,  for  small  values  of  K  the  spectral  method 
is  faster.  In  this  work,  which  is  limited  to  smadl  values  of  K,  the  spectral  method 
is  used. 

5.7  Consistency,  Stability,  and  Convergence  of  the  Numerical  Method 

Finally,  consistency,  numerical  stability,  and  convergence  of  the  numericad 
method  aire  discussed.  A  consistent  numerical  method  is  one  whose  truncation 
error  approachs  zero  as  the  spatial  and  temporal  grid  increments  approach  zero. 
A  numerical  method  is  stable  if  the  roimd-off  error  contained  in  the  numerical 
solution  does  not  grow  with  time.  A  convergent  numerical  method  is  one  in  which 
the  numericad  solution  of  the  discretized  partial  differential  equation  approachs 
the  exact  solution  of  the  partial  differential  equation  as  the  grid  sizes  approach 


zero. 
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The  relationship  between  consistency,  stability,  and  convergence  of  a  dis¬ 
cretization  scheme  is  given  by  the  Lajc  Equivalence  Theorem  (Smith,  1985): 

Given  a  properly  posed  linear  initial- vedue  problem  and  a  linear 
hnite-difTerence  approximation  to  it  that  satisfies  the  consis¬ 
tency  condition,  stability  is  the  necessary  and  sufficient  condi¬ 
tion  for  convergence. 

In  this  work,  the  governing  equations  are  not  linear.  However,  it  is  still  useful  to 
examine  the  consistency  and  stability  of  the  numerical  method. 

The  trimcation  error  of  the  numerical  method  described  in  this  chapter 
is  now  discussed.  Let 

•Lai,As,ai(^)  —  0  (5.132) 

represent  the  difference  equations  (5.42)  and  (5.43)  and  let  be  the  solution  of 
those  difference  equations.  Furthermore,  let  Cl  be  the  exact  solution  of  the  partial 
differential  equations  (5.4).  Then,  the  truncation  error  Et  a*, Ay, At  associated 
with  the  discrete  operator  LAi,Ay,At  is 


Et  Ai,Ay,At  —  L ^x,Ay,At{Cl)  •  (5.133) 

It  can  be  shown  that  the  hybrid  ADI,  Crank-Nicolson,  Adams-Beishforth  dis¬ 
cretization  scheme,  coupled  with  second-order  finite-difference  approximations  of 
the  spatial  derivatives,  has  the  local  truncation  error 

-E7’.A.,Ay.At  =  0(AI^Ay^A^2)  .  (5.134) 

This  truncation  error  approachs  zero  as  Az,  Ay,  and  At  approach  zero  and  there¬ 
fore  the  scheme  is  consistent.  Similarly,  the  discretization  scheme  used  to  solve 
the  Helmholtz  equations  is  consistent  and  exhibits  second-order  spatial  accuracy. 
The  numerical  method  has  spectral  accuracy  in  the  spanwise  direction. 
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Evaluating  the  numerical  stability  of  the  discretization  scheme  used  in 
this  work  is  very  difEcvdt  due  to  the  nonlinearity  of  the  governing  equations.  How¬ 
ever,  Pruett  (1986),  who  used  the  same  discretization  scheme  as  discussed  here, 
found  that  for  simulations  of  instability  waves  in  free  shear  layers,  the  condition 

l«|ma»^  +  +  klma«~  <  1  (5.135) 

was  sufficient  to  ensure  the  numerical  stability  of  the  hybrid  ADI,  Crank-Nicolson, 
Adams-Bashforth  discretization  scheme. 

For  the  simulation  of  instability  waves  in  high-deficit  wakes,  it  is  noted 
that,  because  of  numerical  stability  considerations,  simulations  of  three-dimensionEil 
disturbances  required  a  smaller  time  step  At  than  did  simulations  of  two-dimen¬ 
sional  disturbances.  Furthermore,  although  a  stability  criterion  has  not  been  de¬ 
rived,  for  these  simulations  the  numerical  method  was  stable  for  the  spatial  and 
temporal  resolution  that  was  employed. 
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CHAPTER  6 

RESULTS 


As  discussed  in  Chapter  3,  the  laminar-turbulent  transition  of  wakes  is  influ¬ 
enced  by  many  nondimensional  parameters.  A  thorough  study  of  wake  transition 
would  require  investigating  the  role  of  many  of  these  parameters.  However,  due 
to  the  large  amount  of  computer  time  and  memory  that  is  required  to  solve  the 
Navier-Stokes  equations,  a  detailed  study  of  the  effect  of  each  parameter  on  wake 
transition  is  not  possible. 

In  this  work,  investigations  are  limited  to  two  areas.  First,  the  influence  of 
different  levels  of  excitation  on  the  behavior  of  two-dimensional  disturbances  is 
investigated.  The  results  of  these  investigations  are  discussed  in  Section  6.1.  Sec¬ 
ondly,  the  interaction  of  two-  and  three-dimensional  disturbances  is  investigated. 
This  latter  topic  will  be  discussed  in  Section  6.2. 

Before  proceeding,  the  techniques  used  to  emalyze  the  data  that  result  from 
the  numerical  simulations  are  discussed.  The  harmonic  content  of  the  flow  vari¬ 
ables  is  obtained  by  Fourier  time  series  analysis.  With  representing  any 

one  of  the  flow  quantities,  the  Fourier  analyzed  vEuiables  are  obtained  from 


^*(i,y,F)  =  — - ^  cV>*(x,y,tj)e 

.  «  1  2  L'-L[ 

with  F  =  Q,- — ,- — - - 

^per  iper 


-ilftrFtl 


(6.1) 


The  parameters  L\  and  Xj  are  the  first  and  last  time  steps  respectively  of  the 
time  interval  in  which  the  data  is  Fourier  analyzed.  The  parameter  /per  denotes 
the  number  of  oscillation  periods  that  are  analyzed.  If  one  period  of  oscillation  is 
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analyzed  then  Iper  =  and  if  two  periods  of  oscillation  are  analyzed  then  /per  =  2. 
The  parameter  c  ztssumes  different  values  depending  on  the  flow  variable  that  is 
analyzed.  If  V’*  represents  or  fl*  then  c  =  1.  If  represents  F*  or  flj  then 
c  =  r2.  Finally,  if  represents  IF*  or  fl*  then  c  =  rs.  The  parameter  i  is  the 
imaginary  number 

The  transformed  variable  ^*(x,y,F)  is  the  time  harmonic  component  of 
the  flow  variable  ^*(x,  y,  t).  The  F**^  harmonic  oscillates  in  time  with  the  frequency 
Fl3,  where  (3  is  the  frequency  of  the  wake  excitation. 

The  amplitude  of  the  F***  harmonic  ^*(x,y,  F),  for  fc  =  0,  is 


and  for  1:  >  0  is 


ifF>0 


^  l4|^*(x,y,F)|,  ifF>0. 


(6.2a) 


(6.26) 


The  phase  of  the  F^^  harmonic  ^  *(x,y,F)  is 


<^|j,(x,y,F)  =  arctan 


-Im(t6*(z,y,F)) 
Re(^*(x,y,  F)) 


(6.3a) 


if  V"*  represents  17*,  V*,  or  fi*.  The  phase  of  the  F*^  harmonic  ^*(x,y,F)  is 


<^^(x,y,F)  =  arctan 


/  Re(^*(x,y,F))\ 
\Im(t6*(i,y,F))y 


(6.36) 


if  ^*  represents  IF*,  fl*,  or  fij.  Im(  )  and  Re(  )  denote  the  real  and  imaginary 
parts  of  ^*(x,y,  F)  respectively. 

Additional  quantities  that  are  calculated  are  the  wavenumber  a* ,  the  amplifi¬ 
cation  rate  a*,  and  the  phase  velocity  c*  of  the  various  harmonic  components.  For 


94 


linear  stability  theory,  in  which  the  base  flow  is  assumed  to  be  parallel,  these  quan¬ 
tities  eire  uniquely  defined  for  a  given  disturbance  frequency,  spanwise  wavenumber, 
and  flow  Reynolds  number.  Furthermore,  for  a  disturbance  that  is  governed  by 


linear  stability  theory,  ajl,  a^,  and  c*  can  be  computed  from 


*rLST 


“i  L3T  =  y))  . 


^ilST 


‘^yLST 


(6.4a) 

(6.46) 

(6.4c) 


The  values  of  arLSTt  *^iLSTy  ^plst  independent  of  *  and  y  and  are 
independent  of  which  flow  variable  is  represented  by  V"**  However,  if  expressions 
analogous  to  equations  (6.4)  are  used  to  obtain  similar  quantities  for  a  nonparallel 
flow,  then  the  coordinates  x  and  y  as  well  as  the  flow  variable  that  is  represented 
by  have  considerable  influence  on  the  values  of  a*,  a*,  and  c*.  Gaster  (1974) 
discusses  some  of  the  various  ways  to  calculate  a*,  aj*,  and  c*  for  nonparallel  flows. 

For  this  work,  the  disturbance  quantities  a*,  a*,  and  c*  are  calculated  in  the 
following  manner.  The  streamwise  wavenumber  a*  of  the  F^^  h2umonic  ^*(z,  y,  F) 


is  calculated  from 


aJ(x,F)  =  4t(x,y,F) 


(6.5a) 


where  the  transverse  coordinate  y  is  taken  to  be  constant  in  equation  (6.5a).  The 
^lmplification  rate  a*  of  the  F**‘  harmonic  ^*(x,y,  F)  is  calculated  from 


<»!'(a:,-F’)  =  (ln(^^(x,y,„a,(x),F))  . 


(6.56) 


In  formula  (6.5b),  ymaxi^)  is  the  transverse  coordinate,  as  a  function  of  x,  at 
which  the  amplitude  tl;\{x,y,F)  attains  its  maximum  value.  The  phase  velocity 
of  the  F***  harmonic  ■0*(x,y,  F)  is  computed  from 


c;{z.F)  = 


a*(x,F) 


(6.5c) 
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Finally,  the  disturbance  kinetic  energy  is  calculated  &om 

^  Jyo 

The  kinetic  energy  £*(x,F)  is  a  measure  of  the  disturbance  amplitude.  Also, 
because  the  kinetic  energy  is  an  integral  quamtity,  it  does  not  depend  on  the 
transverse  coordinate  y.  Amplification  rates  of  the  harmonic  components  can  be 
calcxilated  from  the  kinetic  energy  using 

aHx,F)  =  -~(\n(E-{x,F)))  .  (6.7) 

Use  of  equation  (6.7)  to  calculate  of  is  advantageous  because,  unlike  equation 
(6.5b),  it  does  not  depend  on  the  transverse  coordinate  y. 

6.1  Investigations  of  Two-Dimensional  Disturbance  Development 

In  this  section,  investigations  of  two-dimensional  disturbance  development  in 
a  high-deficit  fiat  plate  wake  are  discussed.  Several  different  calculations  are  under¬ 
taken.  With  the  first  two  calculations,  termed  Case-1  and  Case-2,  the  suitability  of 
the  numerical  method  for  the  calculation  of  wake  disturbances  is  tested.  In  Case-1, 
the  method  is  tested  by  comparing  the  results  for  small  amplitude  disturbances  to 
linear  stability  theory.  In  Case-2,  the  results  of  the  calculation  aire  compared  to 
the  experimental  measurements  of  Sato  (1970).  In  Case-2,  the  numerical  method 
is  tested  for  larger  amplitude  disturbances  than  in  Case-1.  In  Case-3  and  Case-4, 
the  large  amplitude  behavior  of  wake  disturbances  is  investigated.  Case-4  is  iden¬ 
tical  to  Case-3  except  that  the  amplitude  level  of  the  excitation,  denoted  by  is 
larger  than  in  Case-3.  By  calculating  disturbances  for  different  excitation  levels, 
the  influence  of  the  initial  amplitude  on  the  disturbance  development  is  observed. 
Finally,  in  Case-5  the  effect  of  the  outflow  boundary  conditions  on  the  results  is 
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investigated.  This  is  done  by  repeating  Case-3  in  a  longer  spatial  domain.  In  all 
of  these  calculations,  the  wake  is  excited  with  a  two-dimensional  sinuous  mode 
disturbance  of  frequency  /3.  Finally,  because  the  calculated  disturbances  are  two- 
dimensional,  only  V®,  and  fl®  are  calculated. 

Case-1 

In  this  case,  the  behavior  of  a  small  amplitude  disturbance  in  a  nonparallel  flat 
plate  wake  is  investigated.  Because  the  base  flow  is  nonpaxallel,  there  is  no  proper 
way  to  compare  these  results  to  linear  stability  theory.  In  fact,  Gcister  (1974) 
states  that  the  agreement  between  linear  stability  theory,  for  which  the  flow  is 
assumed  to  be  parallel,  and  experiments  (numerical  or  physical),  for  which  the 
flow  is  nonparallel,  can  never  be  better  than  0(Re~^ ).  Therefore,  comparisons 
of  the  numerical  results  to  linear  stability  theory  are  qualitative  in  nature  and 
can  not  be  used  to  test  the  accuracy  of  the  numerical  method.  However,  this 
calculation  can  be  used  to  obtziin  some  measure  of  the  suitability  of  the  numerical 
method  for  cadculations  of  small  amplitude  disturbances. 

This  calculation  was  performed  in  a  spatizd  domain  bounded  by 

Ig  =  .3  ,  Xjv  =  .7  ,  j/o  = -10  ,  and  t/m  =  10  .  (6.8) 

The  transverse  extent  of  the  spatial  domain  was  specified  so  that  the  vorticity 
distvirbances  would  be  approximately  zero  at  the  freestream  boimdaries.  For  the 
streamwise  extent  of  the  domain,  it  was  required  that  the  domain  be  long  enough 
to  contain  several  wavelengths  of  the  fundamental  disturbance.  The  domain  is 
^hown  in  Figmes  2.1  and  5.1. 

The  base  flow  wais  calculated  subject  to  the  inflow  streamwise  velocity  distri¬ 
bution  u\y)  given  by  equation  (4.6a).  At  the  inflow  boundary,  the  wake  centerline 
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velocity  Ue  and  the  wake  half-width  b  were 

Ue  =  .5  and  6  =  1.3  .  (6.9) 

The  streamwise  velocity  component  of  the  computed  base  flow  is  shown  in  Figure 

6.1. 

As  discussed  in  Section  4.4,  the  wake  is  excited  using  the  time-dependent 
boundary  conditions  (4.16).  For  this  calculation,  the  disturbance  amplitude  was 
Aj  =  .001.  In  terms  of  the  streamwise  velocity  disturbance  u'(xo,y)  at  the  inflow 
boundary,  A2  is  given  by  the  relation  Aj  =  max(|u'(xo»!/)l)  where  the  function 
max(|u'(xojy)|)  denotes  the  maximum  value  of  |u'(xojy)|  with  respect  to  the  y 
direction.  The  disturbance  frequency  was  =  .317.  For  this  frequency,  the 
amplification  rate  is  less  than  the  maximum  amplification  rate  that  is  predicted 
by  linear  stability  theory.  Furthermore,  the  disturbance  that  corresponds  to  this 
frequency  experiences  amplification  throughout  the  spatial  domain.  The  Orr- 
Sommerfeld  eigenfunctions  that  correspond  to  this  frequency  are  shown  in  Figures 
6.2. 

The  influence  of  Ax  and  Ay  on  the  numerical  results  was  investigated  to 
determine  the  spatial  discretization  that  results  in  solutions  that  are  reasonable 
independent  of  the  grid  sizes.  A  detailed  discussion  of  these  investigations  is  given 
in  Appendix  B.  For  this  calcxilation,  the  grid 

iV  =  128  and  M  =  64  (6.10) 

was  sufficient  for  this  purpose.  For  the  time  step  At,  numerical  stability  rather 
them  temporal  accuracy  was  the  most  severe  constraint.  Therefore,  the  influence  of 
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the  time  step  on  the  numerical  results  was  not  investigated.  To  maintain  numerical 
stability,  the  time  step  was  specified  to  be 

At  =  ^mox  =  64  . 

The  par2imeter  Lmax  denotes  the  number  of  time  steps  per  fundamental  distur¬ 
bance  period.  The  response  of  the  wake  to  the  excitation  at  the  inflow  boundary 
was  calculated  for  five  periods  Tp  of  the  fundamental  disturbance  {L2  =  320) 
where  Tp  =  ^. 

For  this  calculation,  the  spatial  domain  contained  approximately  eight  wave¬ 
lengths  of  the  fimdamental  disturbance.  Each  disturbance  wavelength  was  dis¬ 
cretized  by  approximately  sixteen  grid  points.  In  the  transverse  direction,  the 
spatial  domain  contained  approximately  twenty  momentum  thicknesses  (based  on 
the  inflow  velocity  distribution)  with  approximately  3.2  grid  points  per  momen¬ 
tum  thickness.  As  mentioned  previously,  investigations  of  the  influence  of  the  grid 
sizes  on  the  results  of  the  numerical  calculation  (see  Appendix  B)  have  established 
the  suitability  of  the  spatial  discretization. 


(6.11) 
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The  relevajit  parameters  in  this  calculation  are  summairized  below: 


Case-l: 

the  streamwise  location  of  the  inflow  boundary,  xq  =  .3; 

the  streamwise  location  of  the  outflow  boundary,  =  .7; 

the  transverse  location  of  the  lower  freestream  boimdary,  yo  =  —10; 

the  transverse  location  of  the  upper  freestream  boundary,  yi4  =  10; 

the  Reynolds  number  Re  =  Re  =  200000; 

the  amplitude  level  of  the  wake  excitation,  A2  =  .001; 

the  frequency  of  the  wake  excitation,  ^  =  .317; 

the  number  of  streamwise  grid  increments,  N  =  128; 

the  number  of  transverse  grid  increments,  M  —  64; 

the  number  of  time  steps  calculated,  X2  =  320; 

and  the  time  steps  per  fundamental  disturbance  period,  Lma*  =  64. 


Results  of  this  calculation  are  shown  in  Figures  6.3,  6.4,  6.5,  and  6.6.  Am- 
pliflcation  curves  based  on  the  disturbance  kinetic  energy  J5®(x,F)  are  shown  in 
Figure  6.3.  These  amplification  curves  are  compared  with  analogous  curves  from 
linear  stability  theory  that  are  based  on  the  amplification  rate  a^.  However,  the 
amplification  of  the  kinetic  energy  that  results  from  the  Navier-Stokes  calculation 
arises  from  two  diflerent  sources.  First,  the  kinetic  energy  amplifies  due  to  the 
linear  instability.  Secondly,  the  kinetic  energy  chzmges  due  to  the  slow  divergence 
of  the  base  flow.  However,  despite  the  differences  between  linear  stability  theory 
and  the  calculation,  the  kinetic  energy  of  the  fundamental  disturbance  (F  =  1) 
still  compares  closely  to  the  linear  stability  theory  prediction  of  the  fundamen¬ 
tal  disturbance  kinetic  energy.  The  mean  disturbance  component  (F  =  0)  and 
the  second  harmonic  (F  =  2)  are  also  present,  but  are  much  smaller  than  the 
fundamental  disturbance. 
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For  the  streamwise  location  x  =  .35,  the  amplitude  and  phase  distributions  of 
the  fundamental  disturbance  component  (F  =  1)  are  compared  to  the  amplitude 
2uid  phstse  distributions  of  the  Orr-Sommerfeld  eigenfunctions.  These  compar¬ 
isons  are  shown  in  Figures  6.4  and  6.5.  The  amplitudes  of  the  Orr-Sommerfeld 
eigenfunctions  are  multiplied  by  the  constant 


_  ymaz^l) 


(6.12) 


SO  that  at  y  =  ymax,  they  will  be  exactly  equal  to  the  amplitude  distribution 
of  the  calculated  fundamental  disturbance.  The  phase  distributions  of  the 
Orr-Sommerfeld  eigenfunctions  are  shifted  by  the  constant 


^LST  —  —  *35,  1)  ~  l,s'j'{ymax) 


(6.13) 


so  that  at  y  =  ymax^  they  will  be  equal  to  the  phase  distributions  of  the  calculated 
fundamental  disturbance.  The  constant  ymax  was  defined  previously.  The  function 
’^ALsriy)  the  amplitude  of  the  Orr-Sommerfeld  eigenfunction  and  <l>^isj<iy)  is 
the  phase  of  the  Orr-Sommerfeld  eigenfunction.  As  before,  can  represent  any 
one  of  the  flow  variables. 

As  can  be  observed  in  Figures  6.4,  the  amplitude  distributions  of  the  fun¬ 
damental  disturbance  are  very  similar  tc  lUc  amplitude  distributions  of  the  Orr- 
Sommerfeld  eigenfunctions.  In  Figure  6.4c,  vortidty  disturbances  are  observed 
to  be  qmte  small  (0(10"*))  at  the  freestream  boundaries.  This  confirms  that  the 
transverse  domain  was  wide  enough  for  this  calculation. 

The  comparison  of  the  phase  distributions  of  the  fundamental  disturbance 
to  those  of  the  Ori-Sommerfeld  eigenfunctions  is  shown  in  Figures  6.5.  In  these 
figmes,  reasonably  good  agreement  between  the  niunerical  and  theoretical  phase 
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distributions  is  observed.  The  differences  between  the  ntunerical  auid  theoretical 
phase  distributions  that  are  observed  at  the  traoisverse  locations  y  =  ±10  are  due 
to  the  freestream  boundary  conditions.  For  the  streamwise  and  transverse  velocity 
components,  the  observed  differences  are  due  to  the  exponential  decay  boundary 
conditions,  (5.25a)  and  (5.28a).  For  the  spanwise  vortidty,  the  differences  arise 
because  the  theoretical  spanwise  vortidty  satisfies  exponential  decay  freestream 
boxmdary  conditions  while  the  calculated  spanwise  vortidty  is  set  to  zero  at  the 
freestream  boundaries. 

Finally,  in  Figures  6.6  the  streamwise  wavenumber  a°,  the  amplification  rate 
a®,  and  the  phase  velocity  c®  of  the  fundamentad  disturbance  are  compared  to  the 
analogous  quantities  from  linear  stability  theory.  As  mentioned  previously,  the 
agreement  between  linear  stability  theory  and  the  numerical  experiments  can  not 
be  better  than  0{Re~^)  Gaster  (1974).  Furthermore,  the  values  of  a®,  a®,  and  c° 
depend  on  the  flow  variable  and  transverse  location  that  are  used  for  calctilation 
of  these  quantities.  The  streamwise  wavenumber  a®  and  the  phase  velocity  c° 
are  computed  from  the  flow  variables  17®,  V®,  and  fi®  for  the  transverse  location 
y  =  5.  The  amplification  rate  a®  is  computed  from  all  three  flow  variables  as  well 
as  from  the  kinetic  energy  .F®(x,  1).  Equations  (6.5)  and  (6.7)  are  used  to  calculate 
these  values.  In  Figures  6.6,  the  observed  differences  between  the  theoretical  and 
calculated  values  are  approximately  of  the  order  0{Re~i).  Therefore,  based  on  the 
results  of  Figures  6.6,  the  agreement  between  the  calculation  and  linear  stability 
theory  is  reasonable.  Additionally,  comparable  results  were  obtained  when  other 
trainsverse  locations  were  used  to  calculate  Or  and  Cp. 

The  agreement  between  the  calculated  results  and  linear  stability  theory  is 
within  the  limitations  imposed  by  comparisons  of  paralld  linear  stability  the¬ 
ory  to  nonparallel  numerical  simulations.  The  disturbance  amplification  rate  and 
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wavenumber,  Oi  and  Or,  as  well  as  the  disturbance  amplitude  and  phase  distribu¬ 
tions,  all  exhibit  reasonable  agreement  with  linear  stability  theory.  It  is  concluded 
that  the  numerical  method,  with  the  spatial  and  temporal  discretization  that  was 
employed,  is  suitable  for  the  calculation  of  small  amplitude  disturbances  in  flat 
plate  wakes. 

Case-2 

As  a  further  test  of  the  numerical  method,  the  results  of  this  calculation  are 
compared  to  the  experiments  of  Sato  (1970).  In  contrast  to  Case-1,  in  this  calcula¬ 
tion  the  disturbances  attain  large  amplitude  levels.  To  facilitate  the  comparison  of 
the  numerical  results  to  the  experimental  data,  the  parameters  for  this  calculation 
are  selected  so  that  the  calcTilated  wake  models  the  wake  from  the  experiments  of 
Sato. 

This  calculation  was  performed  in  a  spatial  domain  with  boundaries  at 

xo  =  .03  ,  XN  =  -75  ,  yo  =  —16  ,  and  yM  =  16  .  {6-14) 

As  in  Case-1,  the  transverse  extent  of  the  domain  was  designed  so  that  the  vorticity 
disturbcinces  would  be  approximately  zero  at  the  freestream  boundaries.  In  the 
discussion  of  Case-5,  it  will  be  shown  that  the  streamwise  extent  of  the  domain  is 
sufficiently  long  for  the  present  case. 

The  base  flow  was  computed  subject  to  the  Gaussian  inflow  streamwise  veloc¬ 
ity  distribution  given  by  equation  (4.5a).  This  distribution  was  employed  because, 
as  shown  in  Figure  6.7,  it  compares  well  to  experimentally  obtained  streamwise 
velocity  distributions  of  high-deficit  flat  plate  wakes  (Sato  and  Kuriki,  1961).  At 
the  inflow  boundary,  the  wake  centerline  velocity  Ug  and  the  wsdce  half-width  b 
were 

Uc  =  .234826  and  h  =  1.15  .  (6.15) 
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The  values  of  Ue,  b,  and  zo  were  chosen  so  that  the  computed  base  flow  would 
closely  model  the  wake  in  the  experiments  of  Sato  (1970). 

The  streamwise  velocity  C/®  of  the  computed  base  flow  is  shown  in  Figure 
6.8.  As  observed  in  Figure  6.9,  the  streamwise  variation  of  the  centerline  (y  =  0) 
velocity  of  the  calculated  base  flow  compares  closely  to  the  similarity  solution  of 
Goldstein  (1929).  The  similarity  solution  was  calculated  using  the  series  represen¬ 
tation  for  the  inner  wake  given  by  Goldstein  (1929).  The  large  variation  of  the 
centerline  velocity  indicates  that  nonparallel  effects  may  play  an  important  role  in 
the  initial  development  of  disturbances  in  high-deficit  wakes. 

In  his  experimental  work,  Sato  (1970)  found  that  the  frequency  of  the  pre¬ 
dominant  small  amplitude  disturbance  in  flat  plate  wakes  corresponded  almost 
exactly  to  the  frequency  of  maximum  amplification  as  predicted  by  linear  stabil¬ 
ity  theory.  Therefore,  in  this  calculation  the  undisturbed  wake  was  excited  at  its 
most  unstable  frequency.  This  frequency  was  determined  from  a  linear  stability 
analysis  of  the  inflow  streamwise  velocity  distribution.  The  eigenvalues  (stream- 
wise  wavenumber  and  amplification  rate)  of  the  Orr-Sommerfeld  equation  were 
obtauned  for  a  range  of  frequencies  in  order  to  determine  the  frequency  of  greatest 
instability.  The  streamwise  wavenumber  Qr,  amplification  rate  a^,  amd  phcise  ve¬ 
locity  Cp  obtained  from  this  analysis  are  shown  in  Figures  6.10.  As  seen  in  Figure 
6.10b,  the  amplification  rate  — of  the  sinuous  mode  (mode-1)  achieves  its  max¬ 
imum  value  at  the  frequency  /?  =  .51.  The  eigenfunctions  of  the  Orr-Sommerfeld 
equation,  corresponding  to  the  frequency  ^  =  .51,  are  shown  in  Figures  6.11. 

When  this  calctilation  was  undertaken,  it  was  thought  that  Sato  (1970)  had 
not  indicated  the  initial  disturbance  levels  in  his  experiments.  Therefore,  it  was 
necessary  to  estimate  the  amplitude  level  for  the  disturbance  excitation  in  the 
calcidation.  The  amplitude  of  the  wake  excitation  was  chosen  to  be  Aj  =  .000667 
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as  this  was  thought  to  be  a  good  estimate  of  the  disturbance  level  in  Sato’s  (1970) 
experiments.  However,  after  these  calculations  were  completed,  it  was  learned 
that  the  forcing  level  in  Sato’s  experiments  corresponded  to  A2  %  .001. 

As  in  Caise-1,  the  influence  of  the  grid  sizes  Ax  and  Ay  on  the  numericed  results 
was  investigated  in  order  to  determine  the  spatial  discretization  that  resulted  in 
solutions  that  were  reasonable  independent  of  the  grid  sizes.  A  discussion  of  these 
investigations  is  ^ven  in  Appendix  B.  For  this  calculation,  the  grid 

N  =  1024  and  M  =  256  (6.16) 


was  sufficient  for  this  purpose.  Again,  the  time  step  was  determined  by  stability 
considerations  rather  than  considerations  of  temporal  accuracy.  Therefore,  the 
influence  of  the  time  step  on  the  numerical  results  was  not  investigated.  For  this 
calculation,  the  time  step 


At  = 


Tf 

‘•max 


with  Lmax  =  128 


(6.17) 


was  sufficiently  small  to  ensure  numerical  stability.  The  response  of  the  wake 
to  the  excitation  at  the  inflow  boundary  was  calculated  for  seventeen  oscillation 
periods  Tp  of  the  fundamental  disturbance  {L2  =  2176). 

For  this  calculation,  greater  spatial  resolution  than  in  Case-1  was  required 
because  of  the  presence  of  large  amplitude  harmonics  of  the  fundeunental  distur¬ 
bance.  The  spatial  domain  contained  approximately  thirty  wavelengths  of  the 
fundamental  disturbance  (based  on  the  wavelength  at  the  inflow  boundary),  with 
approximately  34  grid  points  per  wavelength.  In  the  transverse  direction,  the 
spatial  domain  contained  approximately  32  momentum  thicknesses  (based  on  the 
inflow  boundary  velocity  distribution)  with  approximately  eight  grid  points  per 
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momentum  thickness.  As  mentioned,  investigations  of  the  influence  of  the  grid 
sizes  on  the  results  of  the  numerical  calculation  (see  Appendix  B)  have  established 
the  suitability  of  the  spatial  discretization. 

The  relevant  parameters  in  this  calculation  are  summarized  below: 

Case-2: 

the  streamwise  location  of  the  inflow  boundary,  zo  =  .03; 

the  streamwise  location  of  the  outflow  boundairy,  Xff  =  .75; 

the  transverse  location  of  the  lower  freestream  boundary,  yo  =  —16; 

the  transverse  location  of  the  upper  freestream  boundary,  j/m  =  16; 

the  Reynolds  number  Re  =  Re  =  200000; 

the  amplitude  level  of  the  wake  excitation,  Aj  =  .000667; 

the  frequency  of  the  wake  excitation,  /3  =  .51; 

the  number  of  streamwise  grid  increments,  N  =  1024; 

the  number  of  transverse  grid  increments,  M  =  256; 

the  number  of  time  steps  calcidated,  L2  =  2176; 

and  the  time  steps  per  fundamental  disturbeince  period,  ima*  =  128. 

Amplitude  distributions  of  the  fundzimental  disturbance  component  of  the 
streamwise  velocity  were  obtained  by  Fourier  time  series  analysis  in  the  time  in¬ 
terval  15Tf  <  t  <  17Tf-  Comparisons  of  these  lunplitude  distributions  to  those 
from  the  experiments  of  Sato  (1970)  are  shown  in  Figures  6.12.  The  experimen¬ 
tal  results  are  plotted  on  an  arbitrary  scale,  but  the  relative  amplitude  level  at 
each  streamwise  location  is  accurately  represented  in  the  figure.  The  streamwise 
locations  z  for  which  the  amplitudes  in  Figure  6.12a  are  plotted  relate  to  the 
streamwise  locations  X  in  Figure  6.12b  according  to 

_  £ 
t 


(6.18) 
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The  flat  plate  length  I  was  deflned  in  Chapter  3  and  was  equal  to  SOOnam  in 
Sato’s  experiments.  There  is  good  qualitative  agreement  between  the  numeri¬ 
cal  and  experimental  amplitude  distributions  at  the  streamwise  locations  x  = 
.067  {X  =  20mm),  x  =  .1  {X  =  30mm),  and  x  =  .133  {X  =  40mm).  At 
X  =  .200  (JC  =  60mm),  the  shapes  of  the  amplitude  distributions  for  both  the  nu¬ 
merical  simulation  and  the  experiment  are  similar,  but  the  amplitude  level  relative 
to  the  level  at  the  previous  location  (x  =  .133)  is  smaller  for  the  numerical  calcu¬ 
lation  than  for  the  experiment.  This  is  due  to  differences  between  the  excitation 
zunplitude  in  the  calculation  and  in  the  experiment.  Additionally,  it  is  noted  that  in 
Figure  6.12b,  the  experimentally  obtained  amplitude  distribution  for  X  =  20mm 
is  larger  than  the  amplitude  distribution  for  X  —  30mm.  As  this  is  somewhat 
inconsistent  with  expected  behavior,  it  is  possible  that  the  curves  are  mislabeled 
so  that  the  curve  labeled  X  =  30mm  actually  corresponds  to  X  =  20mm  and 
vice-versa.  The  comparison  of  the  numerical  results  to  the  experimental  results 
that  was  discussed  in  this  paragraph  was  bsised  on  the  assumption  that  the  curves 
were  mislabeled. 

The  streamwise  variation  of  the  mean  centerline  velocity  and  the  mean  wake 
half-width  for  both  the  calculation  and  the  experiments  of  Sato  (1970)  are  shown 
in  Figures  6.13.  The  horizont2d  scale  0  <  x  <  2.67  in  Figure  6.13a  corresponds  to 
the  horizontal  scale  0  <  X  <  800mm  in  Figure  6.13b.  The  mean  flow  is  the  zero 
frequency  [F  =  0)  component  of  the  disturbed  wake  and  it  therefore  contains  any 
nonlinearly  generated  0*^  hairmonic  of  the  fimdamental  disturbance. 

In  Figures  6.13,  the  calculated  mean  centerline  velocity  compares  reasonably 
well  to  the  experimentally  obtained  mean  centerline  velocity  (open  circles  in  Figure 
6.13b).  Both  the  numerical  and  experimental  centerline  velocities  increase  rapidly. 
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This  rapid  rise  is  due  to  nonlinear  interactions.  As  observed  in  Figures  6.13,  the 
weike  half- width  for  the  calctilation  is  similar  to  the  half- width  for  the  experiment. 

The  results  of  this  calculation  exhibit  good  qualitative  comparison  with  the 
experimental  results  of  Sato  (1970).  Good  agreement  between  the  numerical  and 
experimental  results  was  obtained  for  large  amplitude  levels  where  nonlineau:  inter¬ 
actions  were  important.  Therefore,  it  is  concluded  that  the  numerical  method  is 
smtable  for  the  calculation  of  large  amplitude  disturbances  in  high-deiicit  wakes. 
Case-3 

The  purpose  of  this  calculation  was  to  investigate  the  effect  of  larger  ampli¬ 
tudes  on  the  disturbance  development.  Therefore,  for  this  calculation  all  parame¬ 
ters  except  for  the  excitation  amplitude  were  kept  the  same  as  in  Case-2.  For  this 
calculation  (Case-3),  the  excitation  amplitude  was  A2  =  .001.  The  base  flow  that 
was  used  as  the  initial  condition  was  identical  to  the  base  flow  from  Ceise-2  (see 
Figures  6.8  and  6.9). 

Figures  6.14  display  the  flow  variables  U®,  V®,  and  fl®  at  the  final  time 
of  the  calculation,  i  =  2176A<.  Due  to  the  rapid  growth  of  the  disturbances, 
the  wake  changes  significantly  from  its  undisturbed  state.  Very  near  the  inflow 
boundary,  the  wake  appears  undisturbed  because  the  disturbances  are  still  quite 
small  in  this  region  of  the  spatial  domzun.  Beyond  x/Ax  k  128  (x  k  .12),  the 
disturbances  become  large  enough  so  that  they  are  the  dominant  feature  in  the 
wake.  Observing  the  spanwise  vorticity  D®  in  Figure  6.14c,  a  pattern  develops 
that  resembles  that  of  a  Karman  vortex  street.  This  pattern  develops  as  the 
disturbances  reaf:h  large  amplitude  levels.  As  mentioned  in  Chapter  3,  the  vorticity 
is  defined  as  u)  =  —V  x  u. 

As  seen  in  Figures  6.14,  at  this  time  the  disturbances  are  still  qiiite  far  up¬ 
stream  of  the  outflow  boundary.  Based  on  the  wavelength  of  the  vortex  street 
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observed  in  Figure  6.14c,  the  leading  edge  of  the  disturbance  wave  is  approxi¬ 
mately  seven  wavelengths  upstream  of  the  outflow  boundary.  Therefore,  at  this 
time  the  assumption  that  the  wake  is  undisturbed  near  the  outflow  boundary  is 
clearly  satisfied.  Furthermore,  the  vortidty  disturbances  are  confined  to  a  region 
nezu:  the  wake  centerline  and  appear  to  be  quite  small  at  the  freestrezun  boundaries. 

In  order  to  see  more  details  of  the  vorticity  field,  vorticity  plots  that  corre¬ 
spond  to  smaller  regions  of  the  spatial  domain  are  shown  in  Figures  6.15  and  6.16. 
Figure  6.15a  shows  the  spanwise  vorticity  fl®  in  the  upstreeim  half  of  the  spatial 
domain  (.03  <  x  <  .39,  0  <  z/As  <  512)  and  Figure  6.15b  shows  the  spanwise 
vorticity  for  the  streamwise  interval  .03  <  x  <  .27  (0  <  z/Az  <  340).  In 
both  of  these  figures,  rapid  disturbance  development  is  observed.  In  Figures  6.16, 
the  spanwise  vorticity  fl®  is  shown  for  the  streamwise  intervals  .03  <  z  <  .12 
(0  <  z/Az  <  128),  .12  <  z  <  .21  (128  <  z/Az  <  256),  .21  <  z  <  .30 
(256  <  z/Az  <  384),  and  .30  <  z  <  .39  (384  <  z/Az  <  512).  These  figures 
show  in  great  detail  the  disturbamces  development  that  results  from  the  wake 
excitation.  In  Figure  6.16d,  a  very  distinct  vortex  street  pattern  is  visible. 

Amplification  curves  based  on  the  kinetic  energy  .E°(z,F)  are  shown  in  Fig¬ 
ure  6.17.  The  kinetic  energy  is  calculated  using  equation  (6.6).  The  harmonic 
content  of  the  flow  variables  wais  obtained  by  Fourier  time  series  analysis  of  the 
velocity  components  in  the  time  interval  13Tf  <  ^  <  15rF'.  For  z  <  .11,  the 
streamwise  variation  of  the  kinetic  energy  of  the  fundamental  disturbance  (F  =  1) 
agrees  closely  with  that  obtained  from  lineau:  stability  theory  calculations.  In 
the  linear  stability  theory  calcidations,  the  streamwise  variation  of  the  base  flow 
was  eurcounted  for  by  the  use  of  a  quasi-uniform  assumption.  The  mean  distur¬ 
bance  component  (F  =  0),  the  second  harmonic  (F  =  2),  and  the  third  harmonic 
(F  =  3)  also  grow  rapidly  in  the  region  z  <  .11.  In  particular,  the  second  and 
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third  harmonics  grow  more  rapidly  than  the  fundamental  disturbance.  Sato  aind 
Kuriki  (1961)  and  Sato  (1970)  have  also  observed  large  mean  and  second  harmonic 
components  in  flat  plate  wakes. 

The  saturation  of  the  fundamental  disturbance  occurs  at  x  .11.  There,  the 
amplitude  of  the  streamwise  velocity  component  of  the  fundamental  disturbance 
is  approximately  twenty  percent  of  the  freestreim  streamwise  velocity.  Beyond 
X  Ai  .11,  the  fundamentad  disturbaince  varies  little  in  the  streamwise  direction. 
The  mean  disturbance  component  amd  the  second  hau^monic  also  saturate  and 
reach  a  state  in  which  they  vary  little  in  the  streamwise  direction. 

The  saturation  of  the  fundamentad  disturbance,  which  does  not  occur  for  the 
small  amplitude,  lineair  development  of  disturbances,  may  nevertheless  be  indi¬ 
rectly  explained  by  the  changing  stability  characteristics  of  the  mean  flow.  As  a 
simple  way  of  determining  the  stability  of  the  meam  flow,  the  spatial  amplification 
rate  — aj(x)  of  the  fundamentad  disturbance  is  computed  from  a  linear  stability 
theory  auiadysis  of  the  mean  flow.  Then,  the  streamwise  variation  of  the  kinetic 
energy  of  the  fundamentad  disturbance  is  computed  from 

E\x,l)  =  i;‘»(xo,l)e~'C^“*‘'*  (6.19) 

using  the  amplification  rates  of  linear  stability  theory.  The  streamwise  variation  of 
the  kinetic  energy  of  the  fundamental  disturbamce,  ais  predicted  by  a  linear  stabil¬ 
ity  analysis  of  the  mean  flow,  is  denoted  in  Figure  6.17  by  the  label  ‘LST  MEAN 
FLOW’.  This  simplified  model  also  exhibits  saturation  of  the  fundamentad  distur¬ 
bance.  However,  the  amplitude  of  the  fundamental  disturbance  after  saturation 
is  over-predicted.  These  results  indicate  that  the  saturation  of  the  fundamentad 
disturbance  is  at  least  partly  due  to  the  changing  stability  characteristics  of  the 
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nonlinearly  generated  mean  flow.  Altematiyely,  the  variation  of  the  mean  flow 
can  be  thought  of  as  being  due  to  the  growth  and  eventual  saturation  of  the 
fundamental  disturbance. 

The  presence  of  the  harmonic  components  F  =  0.5,  F  =  1.5,  and  F  =  2.5 
in  Figure  6.17  is  due  to  the  fact  that  the  wake  has  not  yet  achieved  a  truly  time 
periodic  state  in  the  time  interval  ISTj?  <  i  <  ISTp.  Evidence  of  the  nonperiodic 
nature  of  the  wake  is  given  in  Figure  6.18  which  shows  the  temporal  vairiation  of 
the  treinsverse  velocity  V®  in  the  time  interval  IZTp  <  t  <  IZTp-  It  is  apparent 
from  Figure  6.18  that,  in  the  time  interval  shown,  the  transverse  velocity  V®  is 
not  temporally  periodic  at  z  =  .27.  Instead,  the  amplitude  of  V°  that  corresponds 
to  the  streamwise  location  x  =  .27  appeeirs  to  be  increasing  with  time.  Fourier 
analysis  of  this  data  results  in  the  F  =  0.5,  F  —  1.5,  and  F  =  2.5  harmonics 
that  appear  in  Figure  6.17.  Later,  in  the  discussion  of  Case-5  in  which  a  larger 
streamwise  domain  was  employed,  it  will  be  shown  that  for  the  same  streamwise 
locations  the  flow  variables  eventually  become  temporally  periodic  after  more  time 
has  elapsed.  When  periodicity  is  reached,  the  harmonic  components  F  =  0.5, 
F  =  1.5,  and  F  =  2.5  decrease  to  negligible  levels.  However,  as  will  also  be  shown 
in  Case-5,  the  F  =  0.5,  F  =  1.5,  and  F  =  2.5  harmonics  do  not  influence  the  other 
harmonic  components. 

The  streamwise  variation  of  the  mean  centerline  velocity  and  the  mean  wake 
half-width  are  shown  in  Figures  6.19  and  6.20  respectively.  From  these  figures,  it 
is  apparent  that  the  mean  flow  is  quite  different  from  the  base  flow  and  may  have 
very  different  stability  characteristics. 

The  streamwise  wavenumber  a,,  and  the  phase  velocity  Cp  of  the  fundamental 
disturbance  are  shown  in  Figures  6.21.  They  are  compared  with  the  corresponding 
quantities  that  are  obtained  from  a  linear  stability  analysis  of  the  calculated  mean 
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flow.  The  wavenxunber  and  phase  velocity  that  are  obtained  from  this  analysis 
compare  well  to  the  wavenumber  and  phase  velocity  of  the  calculated  fundamental 
disturbance. 

Case-4 

The  results  of  Case-3  have  shown  that  a  larger  excitation  amplitude  has  signif¬ 
icant  effects  on  the  disturbance  development.  For  this  calcidation,  the  amplitude 
was  increased  by  a  factor  of  ten  over  that  in  Case-S,  resulting  in  A2  =  .01.  The 
base  flow  was  identical  to  that  of  Case-2  and  Case-3  (see  Figures  6.8  and  6.9).  The 
response  of  the  wake  to  the  excitation  at  the  inflow  boundary,  was  calculated  for 
fifteen  fundamental  disturbance  periods  (Xj  =  1920).  All  other  parameters  were 
identical  to  C2ise-3. 

Amplification  curves  based  on  the  disturbance  kinetic  energy  of  the  mean 
component  F  =  0,  the  fundamental  disturbance  F  =  1,  and  the  second  harmonic 
F  —  2  are  shown  in  Figures  6.22.  For  comparison,  the  analogous  curves  of  Case-3 
for  which  A2  =  .001  are  also  displayed  in  these  figures.  The  amplification  curves 
for  the  larger  excitation  level,  A2  =  .01,  saturate  further  upstream  and  at  higher 
amplitude  levels  them  those  for  Case-3.  However,  qualitatively  the  amplification 
curves  are  quite  similar  for  the  two  excitation  levels. 

The  streamwise  variation  of  the  mean  centerline  velocity  and  the  mean  wake 
haJf-width  are  shown  in  Figures  6.23  and  6.24  and  are  compared  with  the  corre¬ 
sponding  curves  from  Case-3.  For  both  excitation  levels,  the  centerline  velocity 
increases  rapidly  due  to  nonlinear  interactions.  However,  for  the  higher  excitation 
amplitude,  the  rapid  increase  in  the  centerline  velocity  begins  further  upstream. 
The  wake  half-width,  shown  in  Figure  6.24,  behaves  in  a  manner  similar  to  the 
centerline  velocity.  The  increase  in  the  wake  half-width,  which  is  also  due  to 
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nonlinear  interactions,  begins  further  upstream  when  the  excitation  amplitude  is 
larger. 

The  queditative  behavior  of  the  disturbed  wake,  as  observed  in  Figures  6.22, 
6.23,  and  6.24,  is  similar  regardless  of  the  excitation  amplitude.  The  main  influence 
of  the  larger  excitation  amplitude  is  to  accelerate  the  onset  of  nonlinear  interactions 
and  the  saturation  of  the  various  disturbance  components. 

Case-5 

With  this  calculation,  the  influence  of  the  outflow  boundary  conditions  on  the 
numerical  results  is  investigated.  For  this,  the  calculation  of  Case-3  is  repeated 
with  a  streamwise  domain  that  is  1.5  times  longer  than  the  one  in  Case-3.  All 
other  parameters  for  this  calculation  were  identical  to  those  of  Case-3.  Thus,  the 
influence  of  the  outflow  boundary  conditions  can  be  assessed. 

This  calculation  was  performed  in  a  spatial  domain  bounded  by 

xq  =  .03  ,  xn  =  1.11  ,  yo  =  —16  ,  and  yM  =  16  .  (6.20) 

The  spatial  domain  was  discretized  into  N  streamwise  grid  increments  and  M 
transverse  grid  increments  with 

N  =  1536  and  Af  =  256  .  (6.21) 

The  spatial  resolution  that  resulted  from  this  discretization  was  identical  to  that 
of  Case-3.  The  response  of  the  wake  to  the  excitation  at  the  inflow  boundary  was 
calculated  for  25  periods  Tf  of  the  fundamentaJ  disturbance  (X2  =  3200)  using 
the  same  time  step  as  in  Case-3.  The  base  flow  was  computed  subject  to  the  same 
inflow  boundary  conditions  as  discussed  in  Case-2. 
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The  relevant  paramaeters  in  this  calculation  are  summarized  below: 

Case-5: 


the  streamwise  location  of  the  inflow  boundary, 

*0 

= 

.03; 

the  streamwise  location  of  the  outflow  boundairy, 

= 

1.11; 

the  transverse  location  of  the  lower  freestream  boimdary. 

yo 

= 

-16; 

the  transverse  location  of  the  upper  freestream  boundary. 

VM 

= 

16; 

the  Reynolds  number  Re  = 

Re 

= 

200000; 

the  amplitude  level  of  the  wake  excitation. 

A2 

= 

.001; 

the  frequency  of  the  wake  excitation. 

= 

.51; 

the  number  of  streamwise  grid  increments, 

N 

= 

1536; 

the  number  of  transverse  grid  increments. 

M 

= 

256; 

the  number  of  time  steps  calculated. 

L2 

3200; 

and  the  time  steps  per  fundamental  disturbance  period, 

^max 

128. 

The  streamwise  variation  of  the  base  flow  centerline  velocity  for  this  calcula¬ 
tion  {N  =  1536)  is  shown  in  Figure  6.25  together  with  the  corresponding  curve  of 
Case-3  {N  =  1024).  It  is  obvious  that  the  centerline  velocity  is  practically  identi¬ 
cal  for  both  cases.  This  is  an  indication  that  the  outflow  boundary  had  negligible 
effect  on  the  base  flow  for  Case-3. 

Contours  of  instantaneous  spanwise  vorticity  for  i  =  2l76Af  and  t  =  3200At 
are  shown  in  Figures  6.26.  For  t  =  2176A<,  the  spanwise  vorticity  contours  are  sim¬ 
ilar  to  those  in  Figure  6.14c  (Ca8e-3).  The  disturbances  experience  rapid  stream- 
wise  amplification  downstream  of  the  inflow  boundary,  attain  relatively  large  am¬ 
plitude  levels,  and  eventually  dominate  the  flow  field.  The  wake  develops  a  Karmin 
vortex  street  pattern  similar  to  what  was  observed  in  Figure  6.14c.  Figure  6.26b 
shows  the  spanwise  vorticity  at  f  =  3200Af.  A  vortex  street  pattern  is  also  visible 
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for  this  time.  However,  the  vortex  street  has  propagated  much  farther  downstream 
than  in  Figure  6.26a. 

In  Figures  6.27,  amplification  curves  for  the  disturbance  kinetic  energy 
E^{xyF)  of  the  mean  component  JP  =  0,  the  fundamental  disturbance  F  =  1, 
and  the  second  harmonic  F  =  2  are  compared  to  the  corresponding  curves  from 
Case-3.  The  disturbance  kinetic  energy  for  Case-3  (N  =  1024)  is  computed  using 
equation  (6.6)  and  the  Fourier  amaJyzed  velocity  components  from  the  time  interval 
13rir  <  <  <  ISTf.  Similarly,  the  disturbance  kinetic  energy  for  Case-5  (JV  =  1536) 
is  computed  using  the  Fourier  analyzed  velocity  components  from  the  time  inter¬ 
val  23rj?  <  t  <  25Tf-  From  Figures  6.27,  it  is  obvious  that  the  amplification 
curves  for  the  various  harmonic  components  corresponding  to  Case-5  {N  =  1536) 
are  almost  identical  to  the  corresponding  curves  from  Case-3  {N  =  1024).  The 
only  differences  are  observed  in  the  region  x  >  .20.  These  differences  are  due 
to  the  different  time  intervals  used  to  Fourier  analyze  the  velocity  components 
for  each  case.  The  otherwise  good  agreement  between  the  results  of  Case-3  and 
Case-5  indicates  that  the  influence  of  the  outflow  boundary  conditions  on  C2ise-3 
is  minimal. 

Figure  6.28  shows  the  amplification  curves  for  the  kinetic  energy  E^{x,F) 
that  results  from  the  Fourier  time  series  analysis  of  the  velocity  components  in 
the  time  interval  23tf  <  i  <  25Tf.  It  is  observed  that  the  harmonic  components 
F  =  0.5,  F  =  1.5,  and  F  =  2.5  are  smaller  by  several  decades  when  compared  to 
the  results  of  Case-3  as  displayed  in  Figure  6.17.  This  reduction  is  due  to  the  fact 
that  the  wake  disturbsmees  are  periodic  in  the  time  interval  23Tf  <  t  <  25Tf‘ 
This  periodic  behavior  is  apparent  from  Fig\ire  6.29  which  shows  the  temporal 
behavior  of  the  transverse  velocity  V°  in  the  time  interval  23riP  <  i  <  25Tf’ 
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Clearly,  V’*’  is  approximately  periodic  for  all  the  streamwise  locations  that  are 
shown  in  Figiire  6.29. 

6.2  Investigations  of  Three-Dimensional  Disturbance  Development 

In  this  section,  investigations  of  three-dimensional  disturbance  development 
in  a  high-deficit  flat  plate  are  discussed.  Several  different  calculations  are  imder- 
taken.  With  the  first  calculation  (Case-6),  the  ability  of  the  numerical  method 
to  accurately  simulate  three-dimensional  disturbances  is  demonstrated.  This  cal- 
ctdation  is  the  three-dimensional  analog  of  C2ise-1.  In  Case-7,  the  response  of 
a  wake  when  subject  to  a  three-dimensional  excitation' is  investigated.  In  con¬ 
trast  to  Case-6,  in  this  calculation  the  base  flow  exhibits  considerable  streamwise 
variation.  In  Case-8  and  Case-9,  the  response  of  the  wake  to  a  combination  of 
two-dimensionsJ  and  three-dimensional  excitations  is  investigated. 

Case-6 

For  this  case,  the  development  of  a  small  amplitude,  three-dimensional  dis¬ 
turbance  in  a  fiat  plate  wake  is  calculated.  The  resiilts  of  this  calculation  are 
compared  to  linear  stability  theory  in  order  to  verify  the  ability  of  the  numerical 
method  to  accurately  simulate  three-dimensional  disturbances.  For  this  calcula¬ 
tion,  the  spatial  2md  temporal  domains,  and  the  base  flow  were  identical  to  Case-1. 

The  wake  was  excited  with  a  three-dimensional  sinuous  mode  disturbance  of 
amplitude  As  =  .001,  spanwise  wave  number  7  =  .5,  and  frequency  =  .28.  The 
two-dimensional  component  of  the  excitation  was  set  to  zero,  so  that  Aj  =  0. 
In  terms  of  the  three-dimensional  streamwise  velocity  disturbance  Uy^{xo,y)  at 
the  inflow  boundary,  As  is  given  by  the  relation  A3  =  max(|u3j(xo,y)|).  The 
amplitude  and  phase  of  the  Orr-Sommerfeld  eigenfunctions  that  correspond  to  the 
frequency  0  =  .28  and  spanwise  wavenumber  7  =  .5  are  shown  in  Figures  6.30. 
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The  spatial  discretization  was  identical  to  that  iised  in  Case-1.  Three  span  wise 
modes  {K  =  6)  were  calctilated.  Restrictions  on  the  allowable  time  step,  due 
to  numerical  stability  considerations,  are  more  severe  for  calculations  of  three- 
dimensional  disturbances  them  for  the  calculation  of  two-dimensional  disturbances. 
For  this  calculation,  the  time  step 


At  = 


with  imax  =  256 


(6.22) 


was  sufficiently  small  to  ensure  numerical  stability.  Five  periods  Tp  of  the  funda¬ 
mental  disturbance  component  were  calculated  {L2  =  1280). 

The  relevant  parameters  in  this  calculation  are  summarized  below: 

Case-6: 


the  streamwise  location  of  the  inflow  boundary. 

Xo 

= 

.3 

the  streamwise  location  of  the  outflow  boundary. 

XN 

= 

.7 

the  transverse  location  of  the  lower  freestream  boundary, 

Vo 

= 

-10 

the  transverse  location  of  the  upper  freestream  boundMy,yjvf 

= 

10 

the  Reynolds  number  Re  = 

Re 

= 

200000 

the  amplitude  level  of  the  two-dimensional  excitation, 

A2 

= 

0 

the  amplitude  level  of  the  three-dimensional  excitation. 

As 

= 

.001 

the  frequency  of  the  wake  excitation, 

0 

= 

.28 

the  spanwise  wave  number  of  the  wake  excitation. 

7 

= 

.5 

the  number  of  streamwise  grid  increments. 

N 

= 

128 

the  number  of  transverse  grid  increments. 

M 

= 

64 

the  number  of  spanwise  modes  computed. 

KI2 

= 

3 

the  number  of  time  steps  calculated. 

L2 

1280; 

and  the  time  steps  per  fundamental  disturbance  period,  Lmax  = 


256. 
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The  amplification  curves  for  the  disturbance  kinetic  energy  E^{x,F)  of  the 
first  spanwise  mode  (k  =  1)  are  shown  in  Figure  6.31  and  are  compared  to  the 
corresponding  curves  from  linear  stability  theory.  The  amplification  curve  for  the 
fundamental  disturbance  (F  =  1)  compares  well  to  the  linear  stability  theory 
prediction  of  the  kinetic  energy.  In  the  Navier-Stokes  calculation,  the  second 
haurmonic  (F  ==  2)  is  present  but  is  much  smaller  than  the  fundamental  disturbamce. 

For  the  streamwise  location  z  =  .35,  the  aunplitude  and  phaise  distributions 
of  the  fundamental  disturbance  component  (k  —  1,F  =  1)  are  compaired  to  the 
amplitude  and  phase  distributions  obtained  from  a  spatial  linear  stability  theory 
analysis  of  the  baise  fiow.  This  compairison  is  shown  in  Figures  6.32  and  6.33.  The 
Orr-Sommerfeld  amplitude  and  phaise  distributions  are  nbrmadized  in  the  same 
mamner  as  discussed  in  connection  with  Caise-1.  The  amplitude  distributions  of  the 
fundamental  disturbance  are  virtually  identical  to  the  Orr-Sommerfeld  amplitude 
distributions.  The  phase  distributions  of  the  fxmdamental  disturbance  also  exhibit 
good  agreement  with  the  Orr-Sommerfeld  phase  distributions. 

The  results  of  this  calculation,  ais  represented  by  the  disturbauice  kinetic  en¬ 
ergy  as  well  as  the  disturbance  amplitude  and  phase  distributions,  compau-ed  rea¬ 
sonably  well  to  linear  stability  theory.  It  is  concluded  that  the  numericad  method 
is  suitable  for  cadculations  of  three-dimensional  disturbances. 

Case-7 

With  this  cadculation,  the  development  of  a  three-dimensional  disturbance  in 
a  high-deficit  flat  plate  wadce  is  investigated.  In  contrast  to  Caise-fi,  the  base  fiow 
changes  significantly  in  the  streamwise  direction  amd  was  approximately  the  saune 
as  the  base  flow  that  was  used  for  Caise-2,  Caise-3,  amd  Case-4.  For  this  caise,  the 
frequency  of  the  excitation  wau  chosen  to  be  identical  to  the  forcing  frequency  for 
Caise-2  through  Case-5  so  that  the  results  of  this  calculation  can  be  compared  to 
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those  earlier  cases.  The  results  of  this  calculation  will  serve  as  a  reference  for  other 
calculations  for  which  both  two-dimensional  amd  three-dimensional  disturbances 
were  introduced. 

This  calculation  was  performed  in  a  spatial  domain  that  is  boimded  by 
iq  =  .03  ,  Xfi  =  .36  ,  yo  =  —16  ,  and  ym  =  16  .  (6.23) 


As  in  Case-2,  the  undisturbed  wake  was  computed  subject  to  a  Gaussian 
inflow  streamwise  velocity  distribution.  The  resulting  base  flow  modeled  the  wake 
from  the  experiments  of  Sato  (1970). 

The  wadce  was  excited  with  a  three-dimensional  sinuous  mode  disturbance  of 
amplitude  A3  =  .001,  spanwise  wavenumber  7  =  .5,  and  frequency  yS  =  .51.  The 
amplitude  of  the  two-dimensional  component  of  the  excitation  was  Aj  =  0.  The 
amplitude  and  phase  of  the  Orr-Sommerfeld  eigenfunctions  that  correspond  to  the 
frequency  /?  and  spanwise  wavenumber  7  are  shown  in  Figiires  6.34. 

The  spatial  domadn  was  discretized  into  N  streamwise  grid  intervals  and  M 
transverse  grid  intervals  where 


N  =  256  and  M  =  128 


(6.24) 


Three  spanwise  modes  {K  =  6)  were  computed.  For  this  calculation,  the  time  step  . 


At  = 


Tf 

^max 


with  Lmax  =  512 


(6.25) 


wais  sufHcient  to  ensure  numerical  stability.  The  time-dependent  response  of  the 
wadke  was  calculated  for  eight  fundamental  disturbance  periods  Tp  which  is  equiv¬ 
alent  to  t  =  4096A^ 

For  this  calculation,  the  spatial  domain  contained  approximately  eight  stream- 
wise  wavelengths  of  the  fundamental  disturbance  (based  on  the  wavelength  at  the 
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inflow  boundary)  with  approximately  sixteen  grid  points  per  streamwise  wave¬ 
length.  In  the  trzmsverse  direction,  the  domain  contained  approximately  32  mo- 
mentmn  thicknesses  (based  on  the  inflow  boundary  velocity  distribution)  with 
approximately  four  grid  points  per  momentum  thickness.  Based  on  the  results  of 
the  two-dimensional  calculations,  it  is  felt  that  the  spatial  discretization  discussed 
here  res\ilts  in  solutions  that  are  sufficiently  independent  of  the  grid  sizes. 

The  relevant  parameters  in  this  calculation  are  summarized  below: 

Case-7: 


the  streamwise  location  of  the  inflow  boundary, 

xo 

= 

.03 

the  streamwise  location  of  the  outflow  boundary. 

Xn 

= 

.36 

the  transverse  location  of  the  lower  freestream  boundary,  yo 

= 

-16 

the  transverse  location  of  the  upper  freestream  boundary, yjvf 

= 

16 

the  Reynolds  number  Re  = 

Re 

= 

200000 

the  amplitude  level  of  the  two-dimensional  excitation. 

M 

= 

0 

the  amplitude  level  of  the  three-dimensional  excitation. 

Ai 

= 

.001 

the  frequency  of  the  W2ike  excitation. 

P 

= 

.51 

the  spanwise  wave  number  of  the  wake  excitation, 

7 

.5 

the  number  of  streamwise  grid  increments, 

N 

= 

256 

the  number  of  transverse  grid  increments, 

M 

= 

128 

the  number  of  spanwise  modes  computed, 

K12 

= 

3 

the  number  of  time  steps  calculated. 

Li 

= 

4096 

and  the  time  steps  per  fundamental  disturbance  period. 

Lmax 

512. 

Amplification  curves  of  the  kinetic  energy  F)  are  shown  in  Figure  6.35a. 
The  disturbance  kinetic  energy  is  calculated  using  the  Fourier  analyzed  velocity 
components  that  correspond  to  the  time  interval  STp  <  t  <  STp.  The  variation  of 
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the  kinetic  energy  of  the  fundamental  disturbance  {F  =  1)  agrees  closely  with  that 
obtained  from  linear  stability  theory  cadculations.  As  wais  the  case  for  calculations 
of  two-dimensionad  disturbances,  the  large  amplitude  level  of  the  F  =  0.5  and 
F  =  1.5  harmonics  is  due  to  the  fact  that  the  wadce  has  not  achieved  a  time 
periodic  state  in  the  time  interval  STp  <  t  <  STjp.  At  i  «  .12,  the  fimdmental 
disturbance  appears  to  saturate.  However,  this  saturation  may  be  due  to  the 
nonperiodic  component  of  the  wake  disturbances. 

For  comparison  purposes,  amplification  curves  of  the  disturbance  kinetic  en¬ 
ergy  that  correspond  to  the  time  interval  STp  <  i  <  lOIV  are  shown  in  Figure 
6.35b.  In  this  figure,  the  fundamental  disturbance  does  not  saturate.  In  addition, 
the  other  harmonic  components,  particularly  the  F  =  0.5  and  F  ~  1.5  harmonics, 
have  decreased  relative  to  their  values  in  Figure  6.35a.  This  indicates  that  the 
disturbances  are  not  periodic  in  the  time  interval  6Tp  <  t  <  STp.  The  saturation 
of  the  fundamental  disturbance  that  was  observed  in  Figure  6.35a  is  apparently 
due  to  this  nonperiodicity. 

Figures  6.36  show  contours  of  streamwise  vorticity  w,  in  the  j/z— plane  for 
t  =  4096 At.  The  horizontal  scale  of  the  plots  in  Figures  6.36,  0  <  z/Az  <  32, 
is  equivalent  to  two  spanwise  wavelengths  of  the  fundamental  disturbance.  The 
vorticity  shown  in  Figure  6.36a  corresponds  to  the  inflow  boundary  x  =  .03  and 
is  due  solely  to  the  three-dimensional  excitation  at  this  boundary.  Figure  6.36b 
shows  the  streamwise  vorticity  corresponding  to  x  =  .146.  For  this  location,  the 
vorticity  concentrations  have  rotated  relative  to  the  vorticity  concentrations  at  the 
inflow  boundary.  The  numerical  results  of  Meibiirg  and  Lasheras  (1988),  which  are 
shown  in  Figure  6.37,  display  a  similar  behavior.  Meiburg  and  Lasheras  attributed 
the  rotation  of  the  vorticity  concentrations  to  the  velocity  field  induced  by  these 
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concentrations.  In  Figiires  6.36,  the  solid  contours  denote  positive  streamwise 
vorticity.  Because  the  vorticity  is  defined  as 

la  =  —V  X  u  ,  (6.26) 

positive  vorticity  induces  a  counterclockwise  azimuthal  velocity  component  and 
a  counterclockwise  rotation  of  the  positive  streamwise  vorticity  concentrations. 
For  similar  reasons,  the  negative  streamwise  vorticity  concentrations  rotate  in  the 
clockwise  direction.  The  vorticity  in  Figures  6.36  is  also  distributed  in  a  more  com¬ 
plicated  pattern  than  the  vorticity  observed  in  Figure  6.37  (Meiburg  and  Lasheras, 
1988).  It  is  believed  that  the  simpler  vorticity  distribution  in  Figure  6.37,  as  com¬ 
pared  to  Figure  6.36b,  is  due  to  the  simpler  initial  perturbation  that  was  employed 
by  Meiburg  and  Lasheras  (they  used  inviscid  vortex  dynamics  and  disturbed  the 
wake  by  perturbing  the  vortex  filaments  sinusoidally  in  the  streamwise  and  span- 
wise  directions). 

Case-8  and  Case-9 

With  these  calculations,  an  attempt  is  made  to  investigate  certeun  aspects  of 
the  secondary  instability  of  wakes.  These  studies  are  undertaken  by  calculating  the 
interaction  of  a  large  amplitude  two-dimensional  disturbance  with  a  smaller  ampli¬ 
tude  three-dimensional  disturbance.  The  results  of  these  calculations  are  compared 
to  the  results  of  Case-7  and  Casc-4  in  order  to  observe  how  the  two-dimensional 
disturbance  influences  the  development  of  the  three-dimensional  disturbance  and 
vice-versa. 

The  base  flow  and  most  relevant  parameters  were  identical  to  those  in  Case-7. 
For  both  Case-8  and  Case-9,  the  amplitude,  frequency,  and  spanwise  wavenumber 
of  the  three-dimensional  excitation  had  the  same  values  as  in  Case-7,  so  that 

Ai  =  .001  ,  =  .51  ,  and  7  =  .5  .  (6.27) 
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For  Case-8,  the  zunplitude  of  the  two-dimensional  excitation  had  the  same  value 
as  in  Case-4,  so  that 

A2  =  .01  .  (6.28) 

For  Case-9,  the  amplitude  and  frequency  of  the  two-dimensional  excitation  W8is 

A2  =  .05  .  (6.29) 

Amplification  curves  of  the  kinetic  energy  E^{x,F  —  1)  for  Case-8  are  shown 
in  Figure  6.38.  The  curve  corresponding  to  ib  =  0  represents  the  kinetic  energy 
of  the  two-dimensional  disturbance.  The  curve  corresponding  to  A:  =  1  is  for 
the  kinetic  energy  of  the  three-dimensional  disturbance.  The  curve  labeled  ‘LST 
MEAN  FLOW’  is  for  the  kinetic  energy  E^{x,F  =  1)  of  the  three-dimensional 
disturbance  component,  as  predicted  by  a  linear  stability  theory  analysis  of  the 
mean  flow.  This  last  curve  is  computed  in  the  same  manner  as  discussed  previously 
(Case-3).  The  two-dimensional  disturbance  {k  =  0)  grows  very  rapidly  at  first  zmd 
then  saturates  at  an  amplitude  of  £°(x,F  =  1)  ss  .20.  This  behavior  is  very  simileir 
to  that  observed  for  the  two-dimensional  disturbance  development  in  Case-4  where 
three-dimensional  disturbances  were  not  present. 

The  three-dimensionaJ  disturbance  (/b  =  1)  initially  grows  in  a  manner  similar 
to  the  three-dimensional  disturbance  from  Case-7.  However,  the  saturation  of  the 
two-dimensional  disturbance  causes  a  temporary  reduction  in  the  amplification 
rate  of  the  three-dimensional  disturbemce.  After  a  brief  period  of  lower  ampli¬ 
fication,  the  three-dimensional  disturbance  resumes  stronger  growth.  However, 
the  amplification  rate  is  smaller  than  the  amplification  rate  that  was  observed 
before  the  saturation  of  the  two-dimensional  disturbance.  The  reduction  in  the 
three-dimensional  growth  rate,  that  is  probably  caused  by  the  saturation  of  the 
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two-dimensional  disturbance,  may  be  due  to  the  changing  stability  characteristics 
of  the  mean  flow.  This  proposition  is  supported  by  the  close  comparison  in  Fig¬ 
ure  6.38  of  the  three-dimensional  disturbance  and  the  curve  labeled  ‘LST  MEAN 
FLOW’.  The  resumption  of  rapid  three-dimensional  disturbance  growth,  follow¬ 
ing  the  saturation  of  the  two-dimensional  disturbance,  could  be  explciined  by  a 
secondary  instability  mechanism. 

Amplification  curves  of  the  kinetic  energy  E^{x,F  =  1)  for  Case-9  are  shown 
in  Figure  6.39.  Due  to  the  now  larger  two-dimensional  excitation  amplitude,  the 
two-dimensional  disturbance  saturates  further  upstream  than  in  Case-8.  As  in 
the  previous  case,  the  saturation  of  the  two-'^imensional  disturbance  causes  a 
strong  reduction  in  the  three-dimepslonai  amplification  rate.  After  a  short  period 
of  reduced  three-dimensional  growth,  the  three-dimensional  disturbance  continues 
its  amplification.  Qualitatively,  this  behavior  is  the  same  as  in  Case-8.  However,  as 
a  result  of  the  larger  two-dimensional  excitation  amplitude,  the  reduction  in  the 
three-dimensional  amplification  rate  and  the  subsequent  resumption  of  stronger 
growth  occur  further  upstream  than  in  Case-8. 

In  Figure  6.39,  the  resumption  of  strong  three-dimensional  growth  is  short¬ 
lived:  the  three-dimensional  disturbance  tends  to  saturate  for  a  second  time.  How¬ 
ever,  this  second  saturation  is  in  contrast  to  the  observed  behavior  in  Figure  6.38 
(Case-8).  To  further  check  the  validity  of  the  results  of  Case-9,  this  calculation 
is  repeated  in  a  spatial  domain  that  is  twice  as  long  {xn  —  *69,  N  =  512)  as 
the  spatial  domain  of  the  original  calculation  but  which  in  all  other  respects  is 
identical  to  the  original  calculation.  The  amplification  curves  of  the  disturbance 
kinetic  energy  E^{x,F)  for  both  the  original  calculation  (xjv  =  -36,  N  =  256)  and 
the  long  domain  calculation  (xjy  =  .69,  iV  =  512)  are  shown  in  Figures  6.40.  For 
the  long  domain  calculation  (Figure  6.40b),  the  three-dimensional  fundamental 


124 


distiirbance  (F  =  1)  initially  grows  more  slowly  after  its  initial  saturation  than 
was  the  case  for  the  original  calcxilation  (Figure  6.40a).  However,  the  fundamental 
disturbeince  does  not  saturate  for  a  second  time  as  it  did  in  the  original  calcula¬ 
tion.  Therefore,  it  appears  that  the  second  saturation  of  the  three-dimensional 
fundamental  disturbance  that  is  observed  in  Figures  6.39  and  6.40a  is  due  to  the 
influence  of  the  outflow  boundary. 

The  three-dimensional  disturbance  behavior  observed  in  Figures  6.38  and  6.39 
is  similar  to  that  observed  by  Metcalfe  et  al.  (1987)  in  their  investigations  of  sec¬ 
ondary  instability  in  free  shear  layers  (see  Figure  6.41).  They  numerically  calcu¬ 
lated  the  temporal  development  of  a  large  amplitude  two-dimensional  disturbance 
as  it  interacted  with  a  smaller  amplitude  three-dimensional  disturbance.  They 
found  that  the  saturation  of  the  two-dimensional  disturbance  temporarily  inhib¬ 
ited  the  amplification  of  the  three-dimensional  disturbance.  After  a  brief  period  of 
reduced  growth,  the  three-dimensional  disturbance  continued  growing  and  even¬ 
tually  surpassed  the  two-dimensional  disturbance. 

For  Case-8  and  Case-9,  the  presence  of  the  two-dimension2il  disturbance  ap¬ 
pears  to  alter  the  distribution  of  vorticity  as  compared  to  Case-7.  The  streamwise 
vorticity  in  the  yz— plane  for  t  =  4096At  is  shown  in  Figures  6.42.  Figure  6.42a 
shows  the  streamwise  vorticity  for  Case-8  which  corresponds  to  the  streamwise 
location  z  =  .146.  Figure  6.42b  shows  the  streamwise  vorticity  for  Case-9  for  the 
same  streamwise  location.  Due  to  the  presence  of  the  two-dimensional  disturbance, 
the  streamwise  vorticity  is  distributed  quite  differently  than  in  Figure  6.36b  (Case- 
7).  In  Figure  6.36b,  the  streamwise  vorticity  is  distributed  in  a  pattern  of  discrete 
streamwise  vortices.  For  Case-8  and  Case-9,  in  which  both  two-dimensional  and 
three-dimensional  disturbances  are  present,  the  streamwise  vorticity  is  distributed 
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in  a  large  number  of  small  vortices.  Furthermore,  in  contrast  to  Case-7,  these 
vortices  do  not  appear  to  rotate. 

Additional  displays  of  the  interactions  of  the  two-  and  three-dimensional  dis¬ 
turbances  are  shown  in  Figures  6.43  through  6.45.  In  Figures  6.43,  the  spanwise 
vorticity  w*  in  the  yz— plane  for  t  =  4096At  is  shown  for  C^tfe-7,  Case-8,  and 
Case-9.  The  interaction  of  the  two-  and  three-dimensional  disturb2uices  (Case-8 
and  Case-9)  significantly  alters  the  vorticity  distribution  as  compared  to  Case-7 
for  which  only  a  three-dimensional  disturbance  was  introduced. 

In  Figures  6.44,  the  streamwise  vorticity  w,  in  the  rj/— plane  for  t  =  4096At 
is  shown  for  Case-7,  Case-8,  and  Case-9.  The  two-dimensional  disturbance,  which 
is  present  in  Case-8  2ind  Case-9,  causes  a  much  broader  transverse  distribution  of 
vorticity  as  compared  to  Case-7  for  which  the  two-dimensional  disturbance  was 
absent.  In  Figures  6.45,  similar  plots  of  the  spanwise  vorticity  u;*  are  shown.  For 
Case-7  (Figure  6.45a)  the  spanwise  vorticity  disturbance  is  much  smaller  thar  for 
Case-8  amd  Case-9  (Figure  6.45b  and  6.45c).  For  Case-8  and  Caise-9,  the  vorticity 
develops  in  a  manner  similar  to  Caise-3  (see  Figure  6.14c). 

Figures  6.46,  6,47,  6.48,  6.49,  and  6.50  (Case-8)  give  detailed  views  of  the 
vorticity  fields  in  the  braid  regions  between  the  large  concentrations  of  spanwise 
vorticity.  It  is  in  these  regions  that  Meiburg  amd  Lasherais  (1988)  observed  the  for¬ 
mation  of  lambda  vortices.  However,  from  the  results  of  these  cadculations,  no  con¬ 
clusive  evidence  of  these  vortices  was  observed.  This  might  be  due  to  the  fact  that 
the  forcing  amplitudes  that  were  used  for  Catse-8  were  different  from  those  used  for 
the  simulations  of  Meiburg  and  Lasheras  (1988).  For  Case-8,  the  waike  was  forced 
with  a  lau’ge  two-dimensional  disturbance  and  a  small  three-dimensional  distur- 
bamce  while  Meiburg  and  Lasherais  employed  a  large  three-dimensional  disturbance 
and  a  smaller  two-dimensional  disturbance.  Furthermore,  there  is  a  considerable 
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difference  between  the  Reynolds  number  for  this  work,  Reo  =  =  594  and 

the  Reynolds  number  for  Meiburg  and  Lasheras’s  calculations,  Res  =  74.  Finally, 
it  was  not  possible  to  duplicate  the  amplitude  levels  used  by  Meiburg  and  Lasheras 
because  a  larger  three-dimensional  forcing  amplitude  would  have  required  the  use 
of  more  Fourier  modes  in  these  calculations.  This  was  not  feasible  due  to  the 
limitations  of  the  available  computer  resources. 

In  Figures  6.51  through  6.53,  the  amplitude  distributions  of  the  spanwise 
vorticity  for  (fc  =  0,i^  =  1)  and  (fc  =  1,F  =  1)  and  the  aunplitude  distributions  of 
the  streamwise  vorticity  for  (A:  =  1,F  =  1)  are  shown.  Case-7,  Case-8,  and  Case- 
9  are  represented  in  these  figures.  Due  to  the  influence  of  the  two-dimensional 
disturbance  in  Case-8  and  Case-9,  the  amplitude  distributions  of  the  streamwise 
and  spanwise  vorticity  are  significantly  different  from  those  for  Case-7  (Figures 
6.51). 

Additional  consequences  of  the  interaction  of  the  two-dimensional  and  three- 
dimensional  disturbances  are  obvious  from  observing  the  behavior  of  the  wake 
half- width  b.  In  Figure  6.54,  the  wake  half-width  b  corresponding  to  three  different 
calculations:  Case-4  {A2  =  .01,  A3  =  0),  Case-7  {A2  =  O^Aj  =  .001),  and  Case- 
8  {A2  =  .01,  A3  =  .001)  is  displayed.  For  reference  the  wake  half-width  of  the 
base  flow  is  also  shown.  When  only  two-dimensional  disturbances  are  present 
(A2  =  .01,  A3  =  0),  the  wake  ais  che^acterized  by  its  half-width  b  becomes  much 
broader.  However,  when  forced  with  only  a  three-dimensionad  disturbance  (Aj  = 
0,A3  =  .001),  the  wake  width  does  not  differ  significantly  from  the  width  of 
the  bcise  flow.  The  strongest  effect  of  broadening  is  observed  when  the  wake  is 
excited  with  both  two-dimensional  and  three-dimensional  disturbances  (Case-8). 
For  Case-8  (A2  =  .01,  A3  =  .001),  the  wake  half- width  development  is  initially 
very  similar  to  that  of  Case-4  (A2  =  .01,  A3  =  0).  However,  while  for  Case-4  the 
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half-width  stops  increasing  when  the  two-dimensional  fundamental  disturbance 
saturates,  for  Case-8  the  half-width  continues  increasing  beyond  this  point.  This 
additional  increeise  in  the  wake  half-width  is  due  to  the  presence  of  the  three- 
dimensional  disturbance. 

In  Figure  6.55,  the  wake  half-width  for  both  Case-8  and  Case-9  is  shown. 
For  Case-9,  in  which  the  two-dimensional  disturbance  is  larger  than  in  Case-8,  the 
initial  increase  of  the  half- width  begins  further  upstream.  However,  the  variation  of 
the  wake  half- width  is  similar  for  both  cases  once  the  two-dimensional  disturbance 


has  saturated. 
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CHAPTER  7 
CONCLUSIONS 

A  numerical  method  has  been  developed  for  studying  the  evolution  of  two-  and 
three-dimensional  disturbances  in  high-deficit  wakes.  Comparison  of  the  results 
of  this  method  were  made  to  both  linear  stability  theory  and  experiments.  The 
numerical  method  was  found  to  be  capable  of  simulating  both  small  and  large 
aunplitude  disturbances. 

Simulations  of  two-dimensional  sinuous  mode  disturbances  in  a  high-deficit 
flat  plate  wake  were  undertaken.  At  small  amplitude  levels,  the  disturbances  grew 
exponentially  at  rates  predicted  by  linear  stability  theory.  At  higher  amplitude 
levels,  nonlinear  effects  became  important  and  the  disturbances  saturated.  The 
saturation  of  the  fundamental  disturbance  was  found  to  be  related  to  the  stabil¬ 
ity  characteristics  of  the  mean  flow.  The  main  influence  of  a  larger  excitation 
amplitude  on  the  resulting  disturbances  was  to  accelerate  saturation.  For  large 
disturbance  amplitudes,  the  wake  devdoped  a  Karmin  vortex  street  pattern.  Fur¬ 
thermore,  the  influence  of  the  outflow  boundary  conditions  on  the  results  of  the 
numerical  simulations  was  foimd  to  be  negligible.  The  subharmonic  component 
that  appeared  in  the  wake  was  found  to  be  due  to  a  nonperiodic  wake  disturbance. 
For  a  fixed  spatial  location,  the  subharmonic  component  decreased  with  time. 

Investigations  of  three-dimensionai  disturbances  were  also  undertaken.  Asso¬ 
ciated  with  three-dimensional  disturbances  were  pairs  of  counter-rotating  stream- 
wise  vortices  that  appeared  to  rotate  as  a  result  of  the  velocity  field  induced  by 
these  vortices.  When  both  two-  and  three-dimensional  disturbances  were  present, 
the  saturation  of  the  two-dimensional  disturbance  caused  the  three-dimensional 
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disturbance  to  saturate.  Shortly  after  saturation,  the  three-dimensionsJ  distur¬ 
bance  resumed  stronger  growth,  possible  because  of  a  secondary  instability  mech¬ 
anism.  Larger  two-dimensional  forcing  amplitudes  accelerated  the  saturation  of 
the  two-dimensional  and  three-dimensional  disturbances  amd  also  accelerated  the 
resumption  of  strong  three-dimensional  growth. 

The  interaction  of  two-dimensional  and  three-dimensional  disturbances  re¬ 
sulted  in  complicated  distributions  of  vortidty.  Instead  of  a  small  number  of  dis¬ 
crete  streamwise  vortices,  as  when  the  wake  was  excited  with  only  three-dimensional 
disturbances,  a  larger  number  of  vortices  spread  over  a  much  wider  transverse  re¬ 
gion  were  present.  Furthermore,  these  interactions  resulted  in  a  much  broader 
wake  distribution  than  that  observed  when  only  two-dimensional  or  three-dimen¬ 
sional  disturbances  were  excited. 

Future  simulations  of  three-dimensional  disturbances  in  wakes  will  undoubt¬ 
edly  require  larger  computational  grids  in  order  to  better  resolve  the  flow  Held. 
For  these  simulations  to  be  practical,  the  numerical  method  should  be  modified.  It 
is  felt  that  greater  efficiency  could  be  obtained  by  solving  the  Navier-Stokes  equa¬ 
tions  in  velocity-pressure  formulation.  This  would  result  in  both  reduced  memory 
requirements  because  a  smaller  number  of  variables  would  have  to  be  stored,  and 
in  reduced  computation  times  because  of  the  smaller  number  of  nonlinear  terms. 

Furthermore,  for  simulations  of  three-dimensional  disturbances,  the  stability 
of  the  current  numerical  method  was  found  to  be  quite  restraining.  This  had 
an  adverse  effect  on  the  computationid  efficiency.  In  the  future,  this  inefficiency 
should  be  avoided  by  employing  numerical  methods  with  more  favorable  stability 
characteristics. 

In  this  work,  interactions  of  two-dimensional  and  three-dimensional  distur¬ 
bances  appeared  to  be  an  important  factor  in  the  three-dimensional  breakdown 
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of  the  wake.  In  these  simulations,  only  the  initial  stage  of  these  interactions  was 
observed.  In  the  future,  more  complete  simulations  of  these  interactions  should 
be  attempted  as  these  interactions  most  likely  play  an  important  role  in  the  de¬ 
velopment  of  the  three-dimensionality  in  transitional  wakes.  Of  particidar  interest 
is  how  these  interactions  lead  to  the  formation  of  the  dominant  three-dimensional 
structures  that  are  observed  in  experiments. 


131 


FIGURES 


Figure  2.1  Spatial  domain. 


SINUOUS  MODE  EIGENFUNCTION  (OMPLITUDE) 
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Figure  6.2  Sinuous  mode  eigenfunctions  of  the  two-dimensional  Oi 
feld  equation  corresponding  to  =  .317  and  Re^  =  ^ 
a)  amplitudes  of  17®,  F®,  and  Oj;  b)  phases  of  £7®,  F® 
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Figure  6.4  Comparison  of  the  amplitude  distributions  of  the  fundamental  dis¬ 
turbance  {F  =  1)  to  linear  stability  theory.  Case-1:  Aj  =  .001, 
A3  =  0,  and  =  .317.  a)  streamwise  velocity,  b)  transverse 
velocity,  V”;  c)  spanwise  vorticity, 
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Figure  6.5  Comparison  of  the  phase  distributions  of  the  fundamental  distur¬ 
bance  (F  =  1)  to  linear  stability  theory.  Case-1:  A2  =  >001, 
i4j  =  0,  and  =  .317.  a)  streamwise  velocity,  17°;  b)  transverse 
velocity,  V°;  c)  spanwise  vorticity,  flj. 
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y  =  y/6 

Fioobe  13.  Non-duneniionat  m«an>v«locity  distribution  in  linear  region. 


Figure  6.7  Comparison  of  a  Gaussian  streamwise  velocity  distribution  to  ex¬ 
perimental  data  and  to  the  similarity  solution  of  Goldstein  (1929), 
reproduced  from  Sato  and  Kuriki  (1961). 


Figure  6.10  Orr>Sommerfeld  eigenvalues  for  the  sinuous  mode  (mode-1)  and 
the  varicose  mode  (mode-2),  corresponding  to  the  strcamwise  ve- 
lodtj  distribution  at  the  inflow  boundary  (Case-2)  and  Rt0  = 
iu^m.  —  594.  a)  streamwise  wavenumber,  b)  amplification 
rate,  Oj;  c)  phase  velocity,  c,. 
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Figure  6.10  Continued. 
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gure  6.11  Sinuous  mode  eigenfunctions  of  the  two-dimensional  O 
feld  equation  corresponding  to  /?  =  .51  and  Rt0  =  ^ 
a)  amplitudes  of  17®,  V®,  and  b)  p* - ^ 
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Frouu  9.  Diatributiona  of  fundamental  (930  Hz)  component.  Natural  transition. 

The  ordinate  scale  ia  arbitrary.  X  (mm):  V.  20;  T>  30;  •,  40;  60;  150. 

Figure  6.12  Comparison  of  streamwise  velocity  amplitude  distributions  of  the 
fundamental  disturbance  (F  =  1)  to  experimental  data,  a)  Case-2: 
Ai  —  .000667,  Ai  =  0,  and  ^  =  .51;  b)  experimental,  reproduced 
from  Sato  (1970). 


Figure  3.  Streamwise  variations  of  velocity  on  the  centre-line  17^,  O 
and  half- value-breadth  b,  #.  Xatural  transition. 

Figure  6.13  Comparison  of  the  streamwise  variation  of  the  mean  streamwise 
centerline  velocity  and  the  meam  wake  half-width  to  experi¬ 
mental  data,  a)  Ca8e-2:  At  s  .000667,  Aj  s  0,  and  =  .51;  b) 
experimental,  reproduced  from  Sato  (1970). 


CM 

n 


Figure  6.14  Contours  of  C/'®,  V®,  and  (1,  for  t  =  2176At.  Case-3:  j4j  = 
.001,  Aj  =  0,  and  =  .51.  Twenty  contour  intervals  between 
the  minimum  and  maximum  values.  Solid  lines  denote  positive 
vorticity  and  dashed  lines  denote  negative  vorticity.  a)  streamwise 
velocity,  C/"®;  b)  transverse  velocity,  V”®;  c)  spanwise  vorticity,  flj- 


Figure  6.14  Continued. 
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Figure  6.15  Contours  of  spanwise  vorticitj  0“  for  i  =  2176At.  Case-3:  Aj  = 
.001,  i43  =  0,  and  ;3  =  .51.  Thirty  contour  intervals  between  the 
minimum  and  maximum  values..  Solid  lines  denote  positive  vortic- 
ity  and  dashed  lines  denote  negative  vorticity.  a)  0  <  x/Ax  <  512, 
b)  0  <  x/Ax  <  340. 
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Figure  6.16  Contours  of  spanwise  vorticitj  for  i  =  2176At.  Case-3;  A2  = 
.001,  Az  =  0,  and  =  .51.  Thirty  contour  intervals  between  the 
minimum  and  maximum  values.  Solid  lines  denote  positive  vortic- 
ity  and  dashed  lines  denote  negative  vorticity.  a)  0  <  x/Ax  <  128, 
b)  128  <  x/Ax  <  256,  c)  256  <  x/Ax  <  384,  d)  384  <  i/Ai  < 
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Figure  6.17  Amplification  curvet  of  the  kinetic  energy  F)  obtained  from 
the  time  interval  ISTjr  <  t  <  15Tp.  Case-3:  A2  =  .001,  Ai  =  0, 
and  0  =  .51. 
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Figure  6.21  Streamwise  variation  of  and  Cp  for  the  fundamental  disturbance 
component  {F  =  1).  Case-3:  A2  =  .001,  A3  =  0,  and  —  .51.  a) 
streamwise  wavenumber,  ar‘,  b)  phase  velocity,  c-. 


Case-3 
F  =  2. 


Figure  6.26  Contours  of  spanwise  Torticity  0®  for  Case-S.  Aj  =  .001,  Aj  =  0, 
$  s  ,51.  Solid  lines  denote  positive  vortidty  and  dashed  lines 


denote  negative  vortidty.  a)  t  =  2176At  and  0  <  i/Ai  <  1024, 
twenty  contour  intervals;  b)  t  =  3200 At  and  0  <  x/Ax  <  1536, 
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Figure  6.30  Sinuous  mode  eigenfunctions  of  the  three>dimen8ionaJ  Orr-Som- 
merfeld  equation  corresponding  to  ^  =  .28,  7  =  .5,  and  = 
=  594.  a)  amplitudes  of  17^,  and  W^\  b)  amplitudes  of 
n^,  nj,  and  fli;  c)  phases  of  and  d)  phases  of 

Oy,  and  n^. 
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Figure  6.32  Comparison  of  the  amplitude  distributions  of  the  fundamental  dis¬ 
turbance  (F  =  1,  ib  =  1)  to  linear  stability  theory.  Case-6:  A2  =  0, 
^3  =  .001,  ^  s=  .28,  and  7  =  .5  a)  streamwise  velocity,  U^;  b) 
transverse  velocity,  V^;  c)  spanwise  velocity, 
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Figure  6.33  Comparison  of  the  phase  distributions  of  the  fundamental  distur¬ 
bance  (F  =  l,fc  =  1)  to  linear  stability  theory.  Case-6:  Aj  =  0, 
As  =  .001,  0  =  .28,  and  7  =  .5.  a)  streamwise  velocity,  b) 
transverse  velocity,  V*;  c)  spanwise  velocity,  . 


Figure  6.34  Sinuous  mode  eigenfunctions  of  the  three-dimensional  Orr-Som- 
m<^eld  equation  corresponding  to  s  .51,  7  s  .5,  and  Rt0  = 
s=  594.  a)  amplitudes  of  and  W^;  b)  amplitudes  of 

fl^,  flj,  and  nj;  c)  phases  of  and  W*;  d)  phases  of  fl^, 

nj,  and 
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Figure  6.35  Amplification  cxirves  of  the  kinetic  energy  il^{zyF).  Chse^lzAt  = 
0,  As  =  .001,  =  .51,  and  7  =  .5.  a)  STp  <  t  <  STjp,  b) 

STj?  <  f  <  lOTjp. 
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Figure  6.36  Contours  of  streamwise  vortidty  w.  in  the  yz-plane  for  <  ss  4096 At. 

Case-7:  Aj  =  0,  i4j  *=  .001,  0  s=  .51,  and  7  *  .5.  Twenty  con- 
totir  intervals  between  the  minimum  and  maximum  values.  Solid 
lines  denote  positive  vortidty  and  dashed  lines  denote  negative 
vortidty.  a)  z  =  .030,  b)  z  =  .146. 
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Figure  6.36  Continued, 
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of  conatant  atreamwiae  vortirity  are  plotted  for  a  reitical  (}’.Z)-piane  at  x  ■  4.6. 


Figure  6.37  Contours  of  streamwise  Tortidty  u>^  in  the  yz>plane,  reproduced 
from  Meiburg  and  Lasheras  (1988). 
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Figure  6.41  Temporal  ▼ariation  of  two-dimensional  and  three-dimensional  dis¬ 
turbance  kinetic  energy  in  a  free  shear  layer,  reproduced  from 
Metcalfe  et.  al  (1987). 


Figure  6.42  Contours  of  streamwise  Tortidty  in  the  yz-plane  for  t  s  4096 At. 

X  =  .146.  Twenty  contour  intemds  between  the  minimum  and 
maximum  values.  Solid  lines  denote  positive  vortidty  and  dashed 
lines  denote  negative  vortidty.  a)  Case-8,  A2  0.01  and  Ai  s 
0.001;  b)  Case-O,  A2  ~  0.05  and  ^4$  s  0.001. 


Figure  6.43  Contours  of  spanwise  Torticity  w,  m  the  yz>plane  for  t  s  4096 At. 

X  s  .146.  Twenty  contour  intervals  between  the  minimum  and 
maximum  values.  Solid  lines  denote  positive  vortidty  and  dashed 
lines  denote  negative  vortidty.  a)  Case*?,  At  —  0.00  and  A3  = 
0.001;  b)  Case-S,  A3  s  0.01  and  A3  »  0.001;  c)  Case-9,  A3  s  0.05 
and  A3  =  0.001. 
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Figure  6.44  Contours  of  streamwise  rortidtj in  the  xy-plane  for  i  s  4096 At. 

z  =  ir/2.  Twenty  contour  interrals  between  the  minimum  and 
maximum  values.  Solid  lines  denote  positive  vortidty  and  dashed 
lines  denote  negative  vortidty.  a)  Case- 7,  At  —  0.00  and  ^43  = 
0.001;  b)  Case-8,  At  —  0.01  and  As  s  0.001;  c)  Case-9,  At  ^  0.05 
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Figure  6.44  Continued 
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Figure  6.45  Contours  of  spanwise  Torticity  u,  in  the  sy>plane  for  t  =  4096At. 

z  =  0.  Twenty  contour  intervals  between  the  minimum  and  majci- 
m\im  values.  Solid  lines  denote  positive  vorticity  and  dashed  lines 
denote  negative  vortidty.  a)  Case-?,  s  0.00  and  —  0.001; 
b)  Case-8,  Ai  s  0.01  and  Ai  s  0.001;  c)  Case-9,  A3  s  0.05  and 
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Figure  6.45  Continued 
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Figure  6.46  Contours  of  spanwise  Torticity  u,  m  the  zy-plane  for  i  ss  4096At 
and  z  =  0.  Case-S:  A2  —  <01,  Aj  =  .001,  0  =  .51,  and  7  s  .5. 
Thirty  contour  intervals  between  the  minimum  and  maximum  val¬ 
ues.  Solid  lines  denote  positive  vortidty  and  dashed  lines  denote 
negative  vortidty. 


Figure  6.47  Contours  of  spanwise  Tortidiy  nf,  in  the  «z>plane  for  t  =  4096At. 

Case-S:  A2  =  .01,  Aj  =  .001,  ^  =  .51,  and  7  =  .5.  Thirty  con¬ 
tour  intervals  between  the  minimum  and  maximum  values.  Solid 
lines  denote  positive  vortidty  and  dashed  lines  denote  negative 
vortidty.  a)  y  s  2  (y/Ay  s=  72),  b)  y  =  -2  (y/Ay  =  56). 
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Figure  6.48  Contours  of  streamwise  Tortidty  in  the  xz-plane  for  t  s  4096At. 

Case-S:  A2  —  .01,  ^43.=  .001,  =  .51,  and  7  =  .5.  Thirty  con¬ 

tour  intervals  between  the  minimum  and  maximum  values.  Solid 
lines  denote  positive  vortidty  and  dashed  lines  denote  negative 
vorticity.  a)  y  =  2  (y/Ay  »  72),  b)  y  =  -2  (y/Ay  =  56). 
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Figure  6.49  Contours  of  spanwise  Tortidiy  w,  in  the  yz-plane  for  i  =  ^96Ai, 
Case-8:  A2  =  <01,  ^4$  s  .001,  =  .51,  and  7  =  .5.  Thirty  con¬ 

tour  intervals  between  the  minimum  and  majumum  values.  Solid 
lines  denote  positive  vortidty  and  dashed  lines  denote  negative 
vortidty.  a)  z  =  .172  (z/Az  =  110),b)  z  =  .160  (z/Az  =  101). 


Figure  6.49  Continued 
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Figure  6.50  Contours  of  streamwise  Tortidty  u;«  in  the  yz-plane  for  t  s  4096 At. 

Case-8:  Aj  =  .01,  j4j  =  .001,  /?  =  .51,  and  7  =  .5.  Thirty  con¬ 
tour  intervals  between  the  minimum  and  maximum  values.  Solid 
lines  denote  positive  vorticity  and  dashed  lines  denote  negative 
vortidty.  a)  *  *  .172  (x/Ai  =  110), b)  *  =  .160  (i/Ax  =  101). 
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Figure  6.51  Vorticity  amplitude  distributions  for  F  =  1.  Case- 7:  A2  —  0, 
Ai  =  .001,  ^  =  .51,  and  7  =  .5.  a)  spanwise  vorticity,  fl®;  b) 
spanwise  vorticity,  c)  streamwise  vorticity, 


Figure  6.52  Vorticity  amplitude  distributions  for  F  =  1.  Case-8:  Aj  =  .01, 
Az  =  .001,  ^  =  .51,  and  7  =  .5.  a)  spanwise  vorticity,  fl®;  b) 
spanwise  vorticity,  fii;  c)  streamwise  vorticity,  17^. 


Figure  6.55  Streamwise  Tariation  of  the  wake  half-width  h  for  Case-8,  A, 
0.01,  A3  =  0.001;  and  Case-9,  A2  —  0.05,  A3  =  0.001. 


Figure  B.l  Influence  of  the  grid  increment  on  the  amplification  of  the  kinetic 
energy  =  1)  for  Case-l.  a)  Af  =  64  with  TV  =  64,  iV  = 

128,  and  TV  =  256;  b)  TV  =  128  with  A/  =  32,  A/  =  64,  and 


Figure  B.2  Influence  of  the  grid  increment  on  the  amplitude  distributions  of 
the  fundamental  disturbance  (F  =  1)  of  the  streamwise  velocity, 
a)  Af  =  64  with  N  =  U,  N  =  128,  and  iV  =  256;  h)  N  =  128 


Figure  B.3  Influence  of  the  grid  increment  on  the  phase  distributions  of  the 
fundamental  disturbance  (F  =  1)  of  the  streamwise  velocity,  a) 
M  =  64  with  N  =  64,  N  =  128,  and  ^■  =  256;  b)  AT  =  128  with 
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Figure  6.4  Influence  of  the  streamwise  grid  increment  on  the  .unplification  of 
the  kinetic  energy  E^{x,F).  Case-2:  M  =  256  with  N  =  512, 
N  =  768,  and  N  =  1024.  a)  the  mean  component,  F  =  0;  b)  the 
fundamental  component,  F  =  1;  c)  the  second  harmonic,  F  =  2. 
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Figure  B.5  Influence  of  the  transvere  grid  increment  on  the  amplification  of 
the  kinetic  energy  E°{x,F).  Ca8e-2:  N  =  1024  with  M  =  128, 
M  =  192,  and  M  =  256.  a)  the  mean  component,  F  =  0;  b)  the 
fundamental  component,  F  =  1;  c)  the  second  harmonic,  F  =  2. 
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Figure  B.6  Influence  of  the  streamwise  grid  increment  on  the  amplitude  dis¬ 
tributions  of  the  streamwise  velocity.  Case-2:  M  =  256  with 
N  =  512,  N  =  768,  and  N  =  1024.  a)  the  mean  component, 
F  =  0;  b)  the  fundamental  component,  F  =  1;  c)  the  second 
harmonic,  F  =  2. 
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Figure  B.7  Influence  of  the  transverse  grid  increment  on  the  amplitude  dis* 
tributions  of  the  streamwise  velodty.  Case-2:  N  =  1024  with 
M  =  128,  M  =  192,  and  M  =  256.  a)  the  mean  component, 
F  =:  0;  b)  the  fundamental  component,  F  =  1;  c)  the  second 
hsirmonic,  F  =  2. 


Figure  B.8  Influence  of  the  streamwise  grid  increment  on  the  phase  distribu¬ 
tions  of  the  streamwise  velocity.  Case-2:  M  =  256  with  N  =  512, 
N  =  768,  and  N  =  1024.  a)  the  mean  component,  F  =  0;  b)  the 
fundamental  component,  F  =  1;  c)  the  second  harmonic,  F  =  2. 


Figure  B.9  Influence  of  the  transverse  grid  increment  on  the  phase  distribu 
tions  of  the  streamwise  velocity.  Ciise-2:  N  =  1024  with  M  =  128 
M  =  192,  and  M  =  256.  a)  the  mean  component,  F  =  0;  b)  th' 
fundamental  component,  F  =  1;  c)  the  second  harmonic,  F  =  2. 
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APPENDIX  A 

SOLUTION  OF  THE  THREE-DIMENSIONAL 
ORR-SOMMERFELD  EQUATION 

In  this  appendix,  the  finite-difference  method  used  to  obtain  the  eigenvalues 
and  eigenfunctions  of  the  three-dimensional  Orr-Sommerfeld  equation  for  the  spa¬ 
tial  stability  problem  is  discussed.  Both  the  Orr-Sommerfeld  equation  and  the 
related  Squire’s  equation  axe  solved  using  fourth-order  finite-differences. 

For  a  small  amplitude  disturbance  with  velocity  and  vorticity  components 


u(x,v,J,i)  =  Rxal(u(j,)e‘<“*+’-«‘>)  , 

(A. la) 

(A.16) 

«;(x,  y,  X,  ()  =  R<!al(ii(y)€‘<"+'"-«‘>)  , 

(A.lc) 

u.(x,V,z,()  =  Real(<x.(y)e«“'+-"-«)  , 

(A.ld) 

w,(x,y,z,()  =  , 

(A.le) 

“.(i.y.x.i)  =  Real(w.(s)e‘<“'+-"-»»)  , 

(A.l/) 

the  amplitude  v{y)  of  the  transverse  velocity  is  an  eigenfunction  of  the  Orr- 
Sommerfeld  equation 

{U-c){v  v-f— -2(a^-t-7^)v  -|-(a^-f-7*)^u)  =  0  .  (A.2) 

axie 

In  equation  (A.2),  U  =  U'(y)  is  the  streamwise  velocity  of  the  imdisturbed  flow 
and  is  assumed  to  be  independent  of  x  and  z.  Furthermore,  a  is  the  streamwise 
wavenumber  of  the  disturbance,  7  is  the  spanwise  wavenumber  of  the  disturbance, 
c  =  /3/q  is  the  phase  velocity  of  the  disturbance  where  is  the  temporal  frequency. 
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Re  is  the  Reynolds  number,  and  i  =  is  the  imaginary  number.  The  notation 
(  )'  denotes  differentiation  with  respect  to  y. 

Solutions  of  equation  (A.2)  are  sought  in  the  transverse  interval  ymin  ^  y  ^ 


Vniax-  For  wakes,  equation  (A.2)  is  subject  to  boundary  conditions  at  y  =  ymin 
and  y  =  ymax  of  the  form 

V  (ymin)  —  Ot  V^ymin)  >  (A. 3a) 

u'(ymin)  =  ot;(y„»in)  ,  (A. 36) 

v"(yma*)  =  Q!^v(yma*)  ,  (A.3c) 

and  v'iymaz)  =  -Q!v(ymo*)  .  (A.3d) 


Equation  (A.2),  coupled  with  the  boundary  conditions  (A. 3),  is  an  eigenvalue 
problem  of  the  form 

Liv  =  CL2V  (A. 4a) 

where 

and 

I2  =  ^  -  7^  .  (A.4c) 

dy* 

For  specified  parameters  a,  7,  Re,  and  the  base  flow  velocity  17,  equation  (A. 4a) 
can  be  solved  for  the  eigenvalue  c  and  the  eigenfunction  v. 

For  spatiadly  amplifying  disturbances,  a  is  obtained  as  a  function  of  /?,  Re, 
and  7  by  iteratively  solving  equation  (A.4a)  for  the  eigenvalue  c,  subject  to  the 
constraint, 

R(a)  =  —  —  c  =  0  .  (A.5) 

a 
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Subsequent  values  of  a  are  obtained  by  applying  the  secant  method  to  equation 
(A. 5)  to  get 

U.6, 

where  j  denotes  the  iteration  level.  The  iteration  scheme  is  repeated  until 
a{l3^Re,j)  is  obtained  for  the  desired  accuracy.  For  each  iteration  level  j,  the 
Orr-Sommerfeld  equation  (A.4a)  is  solved  using  a  finite-difference  method. 

The  eigenfunction  v{y)  of  equation  (A. 4a)  is  sought  at  the  discrete  points 

2/m  =  2/0  +  mAy  for  m  =  0, . . . ,  M  (A. 7) 

where  j/o  =  2/minj  2/M  =  2/max  >  Ay  is  the  spacing  between  grid  points,  and  M  +  1 
is  the  total  number  of  points  in  the  domain.  The  fimction  v{y)  and  its  derivatives 
can  be  approximated  to  fourth-order  accuracy  (Kurtz  and  Cr^lndall,  1962)  using 

,  (A.8a) 
(AM) 

,(A.8c) 
(A.Sd) 

where  =  v(ym)- 


1  7  41  7  1 

1  /  1  ^  2  3  2  1  \ 

^m  —  ^  I  12^"'“^  ^  3Sm-l  -  ) 

and  v^'  =  ^  (gm.2  -  4«7m-i  +  6^,^  -  Ag^-i  +  ^m-2) 


Using  equations  (A.8)  to  approximate  the  derivatives  in  the  Orr-Sommer- 
feld  equation  (A.4a)  and  the  boundary  conditions  (A.3),  the  following  matrix 
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eigenvalue  problem 


is  obtained.  The  coefficients  of  the  matrices  in  equation  (A. 9)  are 


Am  =  U {ym)  ■‘r  <i\U" {ym)  +  0’2  +  >  (A. 10a) 

B„  =  Cf(y„.)  +  03Cf"(j/m)  +  a4  +  &2  ,  (A.106) 

C^  =  C7(y^)  +  a5f^"(y^)  +  a«  +  6,  ,  (A.lOc) 

Dm  =  Bm  ,  (AAOd) 

Em  =  Am  ,  (^-lOe) 

and 


Am  =  ot  > 


(A.llo) 


Brn  =  0»  , 

Cm  =  0»  , 

■Dnv  =  as  , 
Brn  =  Ot  , 


^,  =  ^VS^  +  5  , 


'^  =  ^v'<»=  +  -r"-i  . 

j  3  41/^2/  /  2  I  2N 

*  =  -2--6r(“ 

The  constants  in  equations  (A. 10)  and  (A. 11)  are 


1  ^  /  2  2\ 


360  ’ 


03  = 


3Ay^  45 
7 

=  ’ 

^  /  2  .  2' 

= ‘iZ? ■  55 (“  ' 

41 


“•  =  "5o 


(A.116) 

(yl.llc) 

(A.lld) 

(A.lle) 

(A.12a) 

{A.l2b) 

(A.12c) 

(A.nd) 

(44.12e) 

(44.12/) 

(44.125) 

(A.13a) 

(A.136) 

(44,13c) 

(A.13(i) 

(A.13e) 


(A.13/) 

(A.I35) 
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7^2  2\  2 

“•  =  -«(“  +^)  +  5Z7  ’ 

41  /  2  2\  ^ 

■..  =  -55(0  *y)-^  . 


and 


*  i 

1 

“  aRe  ' 

=  -i-i 

aRe  ' 

\  Ay' 

(  ® 

“  aRe  * 

Uy^ 

3Ay 


(A.13/1) 

(A.13i) 

r) 

,  (A.14a) 

,  (A.146) 

)  • 

(A. 14c) 

The  eigenvalues  «md  eigenfunctions  of  equation  (A. 9)  are  obtained  using  the 
IMSL  routine  EIGZC.  As  equation  (A.9)  is  a  M  +  5  dimensional  matrix  eigenvalue 
problem,  M  +  5  eigenvalues  and  eigenfunctions  are  obtained.  However,  only  two 
eigenvalues  and  two  eigenfunctions  have  physical  meaning. 

Of  the  M  4-  5  values  of  the  complex  eigenvalue  c,  only  those  for  which  the  real 
part  of  c, 

c.  =  R.al(£)=^(^)  .  (A.15) 

satisfies 


Cp  <  1  for  ^  <  1  (A. 16) 

are  considered  to  be  physically  meaningful.  In  this  work,  0  is  always  less  than 
one.  Of  the  eigenvalues  that  satisfy  (A. 16),  one  corresponds  to  a  sinuous  mode 
disturbance  and  one  corresponds  to  a  varicose  mode  disturbance.  The  sinuous 
mode  eigenvalue  is  the  one  for  which 


Ci 


(A.17) 


attains  its  maximum  value  and  for  which  the  corresponding  centerline  transverse 
velocity  is  nonzero.  Alternatively,  the  varicose  mode  eigenvalue  is  the  one  for 


I 
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which  Cj  attains  its  maximum  value  and  for  which  the  corresponding  centerline 
transverse  velocity  is  zero. 

Once  a  is  obtained  for  the  given  parameters  /?,  7,  and  Re,  it  remains  to 
construct  the  velocity  and  vorticity  components  from  the  eigenfunction  Qm-  For 
each  transverse  location  ym,  the  transverse  velocity  Vm  is  obtained  from 


1  7  41  7  1 

~  360^"^“^  60^"*  ^  45^"*“^ 


(A.18) 


Procedures  to  obt<un  the  remaining  velocity  and  vorticity  components  depend  on 
the  whether  the  disturbance  is  two-dimensional  or  three-dimensional. 

For  two-dimensional  disturbances  (7  =  0),  the  streamwise  velocity  is  obtained 
from  the  continuity  equation 

a 

where  v'„^  is  given  by  equation  (A. 8b).  The  spanwise  vorticity  component  Ug  is 
obtained  from 

^zm  =  -  (C  -  (A.20) 

a 

where  is  given  by  equation  (A.8c).  The  remaining  velocity  and  vorticity  com¬ 
ponents  are 

Wm  =  Vxm  =  =  0  •  (A-21) 

For  three-dimensional  disturbances  (7  ^  0),  the  velocity  and  vorticity  compo¬ 
nents  are  obtained  in  the  following  manner.  The  transverse  vorticity  (Jy  is  obtained 
from  Squire’s  equation 


Jy"  —  (a^  4-  7^  +  iaRe{U  —  c))  Jy  =  —iReU'jv  .  (A. 22) 
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Sqmre’s  equation  is  solved  in  the  transverse  interval  ymin  <  y  <  l/ma*  subject  to 
the  Dirichlet  boundary  conditions 

Wj,(ymin)  =  0  ,  (>1.23a) 

and  ojyiymaz)  =  0  •  (A. 236) 

Fourth-order  finite-differences  are  used  to  discretize  the  derivatives  in  equation 
(A. 22).  The  resulting  system  of  equations, 


is  solved  for  the  transverse  vorticity  Uy.  The  coefficients  of  the  matrix  in  equation 
(A. 24)  are 


Gm  =  -12Ay’  (a^  +  7^  +  iaRe{U{ym)  -  c))  -  30 


(A.25) 
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and  the  right  hand  side  of  equation  (A.24)  is  composed  of  the  elements 

fn  =  . 


(A.26) 


The  remaining  velocity  and  vorticity  components  can  be  obtained  from  the 
transverse  velocity  v  and  the  transverse  vorticity  (Jy.  The  spanwise  vorticity  cJ,  is 
obtained  from  the  equation 

,  .  ((«^  +  7^)  Wm  +  -  cu'y.m)  f  .  o-s 


=  ~ta- 


+  7^ 


(A.27) 


Equation  (A.27)  is  derived  from  the  v  Poisson  equation  and  the  relationship 


iQ(Jx  +  +  i-yi^z  =  0 


(A.28) 


In  equation  (A.27),  the  derivative  v'^  is  obtained  from  equation  (A. 8c)  and  wj, 
is  obtained  using  the  fourth-order  finite-difference  formula 


“  1  2/i  '^ym+2) 


(A.29) 


The  streamwise  vorticity  is  obtained  from  equation  (A. 28),  rewritten  as 


_  +  <^y,m 


(A.30) 


where  cJy^  is  obtained  from  equation  (A.29).  The  streamwise  velocity  is  obtained 


from  the  relation 


=  (^)  +  7”) 


(A.31) 


which  is  derived  from  the  definition  of  the  transverse  vorticity.  Finally,  the  span- 


wise  velocity  is  obtained  from  the  continuity  equation 


VJm  =  — 


iaum  +  vL 


(A.32) 


In  both  equations  (A. 31)  and  (A.32),  the  velocity  derivative  v'^  is  obtained  from 
equation  (A.8b). 
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APPENDIX  B 

INFLUENCE  OF  THE  GRID  INCREMENT 
ON  THE  NUMERICAL  SOLUTION 


For  the  calculations  that  are  discussed  in  Chapter  6,  it  is  desired  to  determine 
the  influence  of  the  grid  increment  on  the  numerical  solutions.  This  is  done  by 
repeating  the  calculations  of  Case-1  and  C2ise-2  for  several  different  computational 
grids.  For  Case-1,  small  amplitude  disturbances  are  calculated  and  therefore,  this 
case  is  used  to  test  the  influence  of  the  grid  increment  when  nonlinear  effects  are 
unimportant.  For  Case-2,  larger  amplitude  disturbances,  as  compared  to  Case-1, 
are  calculated.  This  case  is  used  to  test  the  influence  of  the  grid  increment  when 
nonlinear  effects  are  significant. 

For  both  Case-1  and  Case-2,  the  streamwise  and  transverse  grid  increments 
are  varied  independently  so  that  the  variation  of  the  calculated  flow  field  with  the 
grid  increment  may  be  observed.  With  the  domzdn  size  held  constamt,  the  grid 
increments  are  varied  by  changing  the  number  of  grid  points.  For  the  temporal 
discretization,  numerical  stability  considerations  require  that  a  much  smaller  time 
step  be  used  than  is  required  based  solely  on  temporal  accuracy.  Therefore,  the 
time  step  is  not  varied  for  either  Case-1  or  Case-2. 

Case-1  is  recalculated  using  the  grids 


M  =  64  ,  N  =  64, 


(B.la) 

(S.16) 

(R.lc) 

(R.ld) 


t 


M  =  64  ,  N  =  256  , 
N  =  128,  M  =  32  , 
and  N  =  128  ,  M  =  128  . 


274 


In  addition,  the  results  of  Case-1  that  are  described  in  Chapter  6,  for  which 

iV  =  128  and  Af  =  64  ,  (5.2) 

axe  compared  with  the  calculations  described  here.  All  other  parameters  for  these 
calculations  are  identical  to  those  discussed  in  Chapter  6  with  regard  to  Case-1. 

Figures  B.l  show  amplification  curves  of  the  fundamental  disturbance  for 
various  grid  increments.  In  Figure  B.la,  the  amplification  curves  corresponding  to 
AT  =  128  and  N  =  256  are  relatively  close  together  and  exhibit  exponential  growth. 
However,  for  N  =  64  the  amplification  curve  is  significantly  different  from  those 
corresponding  to  N  ■=  128  and  N  =  256  and  does  not  exhibit  exponenticd  growth 
as  would  be  expected  for  a  small  sunplitude  disturbance.  In  Figure  B.lb,  the 
amplification  curves  corresponding  to  M  =  64  and  M  =  128  are  almost  identical 
while  the  amplification  curve  corresponding  to  M  =  32  exhibits  small  differences 
with  the  other  two  curves.  For  all  values  of  M,  the  fundamental  disturbance  grows 
exponentially. 

Figures  B.2  display  the  amplitude  distribution  of  the  fundamental  disturbance 
component  of  the  streamwise  velocity  for  various  grid  increments.  In  Figure  B.2a, 
the  amplitude  distributions  for  N  =  128  and  N  =  256  are  quite  close  and  compare 
well  to  linear  stability  theory.  For  N  =  64,  the  shape  of  the  amplitude  distribution 
is  similar  to  those  for  N  =  128  and  N  =  256,  but  its  3unplitude  level  is  significantly 
higher  than  for  the  other  curves.  In  Figure  B.2b,  the  eunplitude  distributions  for 
M  =  64  and  M  =  128  axe  quite  close.  For  M  =  32,  the  amplitude  distribution  is 
similar  in  shape  to  those  for  M  =64  and  M  =  128  but  clearly  requires  a  smaller 
transverse  grid  increment. 

Figures  B.3  display  the  phase  distribution  of  the  fundamental  disturbance 
component  of  the  streamwise  velocity  for  various  grid  increments.  In  Figxire  B.3a, 


the  phase  distributions  for  N  =  256  and  N  =  128  are  similar  in  shape  and  rea¬ 
sonably  close  together.  For  N  =  64,  a  significant  phase  difference  between  this 
curve  and  those  for  N  =  128  and  N  —  256  is  present.  In  Figure  B.3b,  the  phase 
distributions  for  zdl  values  of  M  are  in  reaisonable  agreement. 

From  these  results  it  is  apparent  that  the  size  of  the  grid  increment  has 
significant  influence  on  the  calculated  flow  field.  However,  it  appears  that  for  the 
grid  used  for  Case-1  in  Chapter  6,  for  which 

Af  =  64  and  N  =  128  ,  (B.3) 

the  calculated  flow  field  is  reasonable  well  resolved  and  is  not  significantly  altered 
when  the  finer  grids  {N  =  256,  A/  =  64  and  N  =  128,  Af  =  i28)  are  employed. 


Case-2  is  recalculated  using  the  grids 

M  =  256  ,  =  512  ,  (BAa) 

M  =  256  ,  N  =  768  ,  {BAb) 

N  =  1024  ,  A/  =  128  ,  (19.4c) 

and  N  =  1024  ,  Af  =  192  .  {BAd) 

In  addition,  the  results  of  Case-2  that  are  discussed  in  Chapter  6,  for  which 

M  =  256  and  N  =  1024  ,  {B.5) 


are  compared  with  the  calculations  discussed  here.  All  other  parameters  for  these 
Ccilculations  are  identical  to  those  discussed  in  Chapter  6  with  regard  to  Case-2. 

Since  nonlinear  effects  are  important  for  this  flow  field,  the  fundamental  dis¬ 
turbance  as  well  as  the  mean  and  second  heirmonics  components  are  displayed  for 
various  grid  increments.  Figures  B.4  and  B.5  display  the  amplification  curves  for 


276 


the  mean  component  F  =  0,  the  fundamental  disturbance  F  =  1,  and  the  second 
harmonic  F  —  2-  In  Figxire  B.4a  (F  =  0)  and  B.4b  (F  =  1),  the  amplification 
curves  corresponding  to  N  =  768  and  N  =  1024  are  virtually  identical.  Beyond 
the  point  of  saturation,  the  amplification  curves  for  N  =  512  deviate  from  those 
corresponding  to  N  =  768  and  N  =  1024.  In  Figure  B.4c,  the  second  harmonic 
(F  =  2)  appears  to  be  more  sensitive  to  the  streamwise  grid  increment  than  the 
other  harmonic  components  and  appears  to  require  greater  streamwise  resolution 
than  was  used  in  these  calculations.  In  Figures  B.5,  the  amplification  curves  for 
all  three  harmonic  components,  F  =  0,  F’  =  1,  and  F  =  2,  are  virtually  identical 
regardless  of  the  transverse  grid  increment. 

Figures  B.6  through  B.9  display  the  amplitude  and  phase  distributions  of  the 
streamwise  velocity  for  F  =  0,  F  =  1,  and  F  =  2.  As  was  observed  for  the 
amplification  curves,  the  amplitude  and  phase  distributions  vary  more  with  the 
streamwise  grid  increment  than  with  the  transverse  grid  increment.  In  particular, 
the  second  harmonic  component  displays  the  most  sensitivity  with  respect  to  the 
streamwise  grid  increment. 

From  the  results  presented  in  Figures  B.4  through  B.9,  it  would  appear  that 
for  M  =  256,  the  transverse  grid  increment  is  sufficiently  small  so  that  the  mean 
component,  the  fundamental  component,  and  the  second  harmonic  component  are 
adequately  represented.  However,  for  the  streamwise  grid  increment  corresponding 
to  N  =  1024,  the  results  are  not  quite  as  conclusive.  It  appears  that  the  mean 
component  and  the  fundamentsd  component  are  adequately  represented  with  this 
grid  increment.  However,  it  is  not  quite  eis  clear  whether  the  second  harmonic  is 
adequately  resolved  as  it  displays  a  significant  variation  with  the  streamwise  grid 
increment. 
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